Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,484,002 A→G E317E (GAA→GAG)  unknown → Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,484,0020AG98.0% 137.8 / ‑5.2 50E317E (GAA→GAG) unknownType II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (23/26);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

TTTGGAACAGCCTCATGTTGTACGCCTGGAAACGCGTTTAGCCGGCGTTGAGCATACAGGCGAAGTTACTATGCAAGACTTAGTAATAAATGCATTGCGTATGCGTCCTGAGCGTATTATTGTGGGCGAATGCCGTGGTGGAGAAGCATTCCAAATGTTACAAGCAATGAATACCGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATAGCCCTCGCGATGCAACATCTCGCTTGGA‑A‑AGTATGGTTATGATGTCGAATGCATCCCTACCGTTAGAGGCAATTCGCCGTAATATCTCTTCAGCAGTAAACATAATAGTGCAGGCGTCCCGTTTAAATGATGGTTCAAGAAAAATTATGAATATCACTGAGGTTATGGGAATGGAAAACGGTCAAATCGTTTTACAGGATATATTTTCCTATAAGGCAAGTAA  >  1/1483762‑1484196
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                   
tttGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTa                                                                                                                                                                                          <  1:452672/251‑1 (MQ=255)
tttGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTa                                                                                                                                                                                          <  1:383791/251‑1 (MQ=255)
tttGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTa                                                                                                                                                                                          <  1:323908/251‑1 (MQ=255)
         accTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCg                                                                                                                                                                                  <  1:475616/249‑1 (MQ=255)
         accTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGa                                                                                                                                                                                 >  1:782473/2‑251 (MQ=255)
              aTGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGt                                                                                                                                                                                            <  1:655968/235‑1 (MQ=255)
                gttgtACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGcgt                                                                                                                                                                           >  1:407762/1‑248 (MQ=255)
                             aaaCACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGttggaa                                                                                                                                                             >  1:41889/1‑247 (MQ=255)
                                            gCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGt                                                                                                                                                                                    >  1:643934/1‑213 (MQ=255)
                                                tGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattg                                                                                                                                          <  1:117009/251‑3 (MQ=255)
                                                    cATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgca                                                                                                                                        <  1:292012/249‑5 (MQ=255)
                                                          ggCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGttgga                                                                                                                                                              <  1:15671/221‑4 (MQ=255)
                                                           gCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcag                                                                                                                                 <  1:538650/249‑12 (MQ=255)
                                                            cgaggTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACGCTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaa                                                                                                                              <  1:28129/247‑15 (MQ=255)
                                                                 ttACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaatatc                                                                                                                         <  1:689306/251‑20 (MQ=255)
                                                                 ttACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaatatc                                                                                                                         <  1:215324/251‑20 (MQ=255)
                                                                   aCTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGACGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTa                                                                                                                                                                                          <  1:754580/184‑1 (MQ=255)
                                                                                           gcactacGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatat                                                                                                                                            <  1:872275/200‑1 (MQ=255)
                                                                                              ctacGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaatatc                                                                                                                         <  1:275003/219‑20 (MQ=255)
                                                                                                acGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaat                                                                                                                            <  1:743762/216‑17 (MQ=255)
                                                                                                acGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatc                                                                                                                                      <  1:659428/206‑7 (MQ=255)
                                                                                                acGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatat                                                                                                                                            >  1:776148/2‑201 (MQ=255)
                                                                                                  gTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTa                                                                                                                                                                                          <  1:257180/153‑1 (MQ=255)
                                                                                                        gTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGcgtc                                                                                                                                                                          >  1:175010/1‑160 (MQ=255)
                                                                                                           cTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAAtc                                                                               <  1:281034/251‑2 (MQ=255)
                                                                                                            tGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAAtca                                                                              <  1:782966/251‑3 (MQ=255)
                                                                                                                     tatcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAatatta                                                                     >  1:405504/11‑249 (MQ=255)
                                                                                                                      atcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGAGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattg                                                                                                                                          >  1:807965/10‑181 (MQ=255)
                                                                                                                      atcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAatattac                                                                    <  1:563338/242‑4 (MQ=255)
                                                                                                                      atcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAatatta                                                                     >  1:629222/10‑248 (MQ=255)
                                                                                                                      atcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAatatta                                                                     >  1:418153/10‑248 (MQ=255)
                                                                                                                       tcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAGCACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTatg                                                                                                                                                                                        >  1:346730/9‑134 (MQ=255)
                                                                                                                       tcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTAcc                                                                                                                                                                    >  1:706188/9‑154 (MQ=255)
                                                                                                                                         gtggAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTatga                                                                                                                                                                                       >  1:330552/1‑117 (MQ=255)
                                                                                                                                          tggAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGaaa                                                <  1:462289/251‑1 (MQ=255)
                                                                                                                                          tggAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGaa                                                 <  1:676093/250‑1 (MQ=255)
                                                                                                                                              gaAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGaaaatgg                                            >  1:390491/1‑248 (MQ=255)
                                                                                                                                                    ttCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaat                                                                                                                            <  1:316392/165‑17 (MQ=255)
                                                                                                                                                    ttCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAatattac                                                                    >  1:846452/1‑218 (MQ=255)
                                                                                                                                                    ttCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAAt                                      >  1:145787/1‑251 (MQ=255)
                                                                                                                                                    ttCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAAt                                      >  1:775030/1‑251 (MQ=255)
                                                                                                                                                                   gcgaTGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatat                                                                                                                                            >  1:12367/4‑134 (MQ=255)
                                                                                                                                                                    cgaTGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcat                                                                                                                                       <  1:30257/136‑6 (MQ=255)
                                                                                                                                                                     gaTGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGaaataatgggtatgga                                                  >  1:737220/2‑207 (MQ=255)
                                                                                                                                                                         aacactggTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAATTGTTTTACAAGATATTTTTTc                 <  1:401941/245‑1 (MQ=255)
                                                                                                                                                                                    cGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAATTGTTTTACAAGATATTTTTTCTTATAAGGCaa      >  1:105586/1‑251 (MQ=255)
                                                                                                                                                                                    cGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAATTGTTTTACAAGATATTTTTTCTTATAAGGCaa      >  1:833597/1‑251 (MQ=255)
                                                                                                                                                                                        ggTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTAc                                                                                                                                                                     >  1:218801/1‑88 (MQ=255)
                                                                                                                                                                                        ggTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAGAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAATTGTTTTACAAGATATTTTTTCTTATAAGGCAAGTaa  >  1:438528/1‑251 (MQ=255)
                                                                                                                                                                                        ggTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAATTGTTTTACAAGATATTTTTTCTTATAAGGCAAGTaa  >  1:323313/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                   
TTTGGAACAGCCTCATGTTGTACGCCTGGAAACGCGTTTAGCCGGCGTTGAGCATACAGGCGAAGTTACTATGCAAGACTTAGTAATAAATGCATTGCGTATGCGTCCTGAGCGTATTATTGTGGGCGAATGCCGTGGTGGAGAAGCATTCCAAATGTTACAAGCAATGAATACCGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATAGCCCTCGCGATGCAACATCTCGCTTGGA‑A‑AGTATGGTTATGATGTCGAATGCATCCCTACCGTTAGAGGCAATTCGCCGTAATATCTCTTCAGCAGTAAACATAATAGTGCAGGCGTCCCGTTTAAATGATGGTTCAAGAAAAATTATGAATATCACTGAGGTTATGGGAATGGAAAACGGTCAAATCGTTTTACAGGATATATTTTCCTATAAGGCAAGTAA  >  1/1483762‑1484196

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑