Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,576,080	0	G	A	42.6%	88.7 / 23.3	61	L514F (CTT→TTT)	yddB	putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base G (30/5);  new base A (0/26);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

CGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAATGAC  >  NC_000913/1575931‑1576228                                                                                                                                                         |                                                                                                                                                     cgcAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGa                                                                                                                                                     <  2:485549/151‑1 (MQ=255)  gcAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGaa                                                                                                                                                    <  1:976192/151‑1 (MQ=255)      tAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATAtt                                                                                                                                                >  1:866646/1‑151 (MQ=255)       aaCGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATAtt                                                                                                                                                <  2:769999/150‑1 (MQ=255)             cTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGcc                                                                                                                                          >  2:880749/1‑150 (MQ=255)               tACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCAt                                                                                                                                        >  1:214140/1‑150 (MQ=255)                    aTTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAAc                                                                                                                                  >  2:859245/1‑151 (MQ=255)                           gctgcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTtat                                                                                                                            >  1:372598/1‑150 (MQ=255)                            ctgcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTtata                                                                                                                           >  2:99742/1‑150 (MQ=255)                            ctgcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTtata                                                                                                                           >  2:72388/1‑150 (MQ=255)                            ctgcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTtata                                                                                                                           >  2:990526/1‑150 (MQ=255)                            ctccAACCCCATTCCTATTTAATCTTTATAGGGTTTTTAAAAA‑CCCTGATAGCGCTCCAGATTTTTTTTCCCTGATCCGGATCCCTT‑CCATCTTTGGGGAATACCACTAAACCCAAT‑AAAAAAAATCTTCCCGCCATAGTAACGGTtata                                                                                                                           <  1:967697/147‑1 (MQ=255)                            ctccAACCCCATGGCCAGTTCATCGTTATAAGGTGTTTTAAAT‑TCCTGATAAGGCGTAAGTTTTTTATTACCTGAT‑CCCATTCCCTCCCACCTATGGGGAAAATCACGTAACCCAAT‑ACAAAAAATATTCCCGCCATAGTAACGGtt                                                                                                                              <  1:335739/144‑1 (MQ=255)                              gcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAAcc                                                                                                                        >  1:685282/1‑151 (MQ=255)                                aaTCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACctg                                                                                                                      >  2:621545/1‑151 (MQ=255)                                 aTCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACctgc                                                                                                                     >  1:430665/1‑151 (MQ=255)                                   cccATTGCCAGTTTATCTTTATAAGGTTTTTAAAAA‑TCCTAATAACGCGTCAGATTTTTTTTACCTGATACCGTTTCCGC‑CCAGCTTTTGCGGATATCACGTAACCCCAT‑ATCAAAAATTTTCCCGCCATAGTAACGGTTATAACctgctg                                                                                                                   <  2:604170/151‑1 (MQ=255)                                   cccATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACctgctg                                                                                                                   >  1:194120/1‑151 (MQ=255)                                       ttGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAAt                                                                                                               >  1:439894/1‑151 (MQ=255)                                        tGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATc                                                                                                              >  1:688692/1‑151 (MQ=255)                                          ccAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCAt                                                                                                            >  2:562882/1‑151 (MQ=255)                                           cAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCAtt                                                                                                           >  2:730548/1‑151 (MQ=255)                                              ttcttcTTTATAGGGCTTTTTAAAA‑TCTTGAAACCGCGTAAGATTTTTTTTCCCTGATCCCGTTCCCGT‑CCAGCTATGGGGGATATCCCTTAACCCCAA‑ATCAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGa                                                                                                         <  2:176878/150‑1 (MQ=255)                                              ttCCTCGTTATAAGGTGTTTACAAT‑CCCTGATAACGCGTAAAATTTTTTTTCCCTGAAACCGTTTCCGC‑CCGGCTTTCGCGAATACCACGTAATCCCAT‑ATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGa                                                                                                         <  2:459761/150‑1 (MQ=255)                                              ttCATTTTTATAGGGCTTTTACAAT‑TCTTAATAACGCGTCAGTTTTTTTTTACCTGATACAGATCCCCT‑CCAGCTATGGCGGATATCACGTACCCCAAT‑ATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGa                                                                                                         <  1:385758/150‑1 (MQ=255)                                              ttCATCTTTATAAGGCTTTTTAAAA‑TCTTGAAAACGCGTCAGATTTTTTTTCCCTGATCCCGATCCCTC‑CCATTTATGGCGGATACCACTTAACCCCAT‑ATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTg                                                                                                          <  2:452930/149‑1 (MQ=255)                                              ttCATCATTATCGGGC‑TTTTCAAACTCCTAAAAACGTATAAGATTTTTATTACCTGACCCAGATTCCTA‑CTAGTTATGGTGGTTAACAGGAAACCCAAC‑TTAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGa                                                                                                         <  2:994699/150‑1 (MQ=255)                                                 aTCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGtt                                                                                                     >  2:238751/1‑151 (MQ=255)                                                       aaagggCTTTTTCAAA‑TCGTGATACCGCGTCAGTTTTTTTTTCCCTGATCCCGATTCCTC‑CAAGTTATGGCGGATATCACGTAACCCAAT‑ATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTg                                                                                                   <  2:427065/143‑1 (MQ=255)                                                       aaagggCTTTTTCAAA‑TCCTGATACCGCGTCAGATTTTTATTCCCTGATCCCTATCCCGT‑CCAGCTTTTGCGAATATCACGTAACCCAAT‑ATCAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAt                                                                                                 <  2:248868/145‑1 (MQ=255)                                                         aaGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGc                                                                                             >  2:726059/1‑151 (MQ=255)                                                              gTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                                         >  2:172856/1‑150 (MQ=255)                                                               ttttCAAA‑GCCTAATAACGCGTCATTTTTTTATTCCTTGATACTTTTCCCTT‑CCAGCTTTTGCGGATATCACATAACCCAAT‑ATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                                         <  1:1025221/149‑1 (MQ=255)                                                               tttaaaaa‑CCCTGATAACGCGTAAGGTTTTTATTCCCTAATCCCATTCCCGT‑CCAGCTATTGGGAATACCACGTAACCCAAT‑AAAAAAAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                                         <  1:909989/144‑1 (MQ=255)                                                                  tacaattctttaTACCGCGTAAGATTTTTATTACCTAATCCCGATCCCGC‑CCATCTATTGGGAATTTATCTTAACCCAAT‑ATAAAAATTTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAata                                                                                      <  1:731169/139‑1 (MQ=255)                                                                   cAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATcc                                                                                   >  2:694590/1‑151 (MQ=255)                                                                     aa‑TCCTAATAACGCGTCAGATTTTTATTCCCTAATCCCGTTTCCGC‑CCAGCTATTGCGGATATCACTTAACCCCAT‑ATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCa                                                                                 <  2:440262/151‑1 (MQ=255)                                                                             aTAACGCGTAAGATTTTTATTCCCTGATACCGATCCCGT‑CCAGCTATGGCGAATATCACGTACCCCAAT‑ATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTc                                                                          <  2:445183/151‑1 (MQ=255)                                                                                 cgcgTCAGTTTTTTTTTCCCTGATCCCTATCCCGC‑CCAGCTATTGCGAATATCCCGAACTCCCAT‑AAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCAt                                                                        <  1:654411/149‑1 (MQ=255)                                                                                   cgTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATaa                                                                      >  1:110880/1‑149 (MQ=255)                                                                                           ttttttTCCCGGATCCCGATCCCGT‑CGAGTTTTGGCGGATACCGCTAACCCCCAT‑ACAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGcgcg                                                                <  2:93337/147‑1 (MQ=255)                                                                                           ttttATTCCCTGATCCCGATCCCGT‑CCAGCTTTTGCGAATATCACGTAACCCCAA‑ATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGcgcg                                                                <  2:787638/147‑1 (MQ=255)                                                                                            ttttttCCCTGAACCTGATTCCGT‑CCAGTTATTGCGAATACCACGTACCCCCAT‑AAAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGga                                                            <  1:650562/150‑1 (MQ=255)                                                                                             ttattaCCTAATCCCGATTCCGC‑CCAGTTATTGCGATTATCACGTAACCCCAT‑AAAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGgag                                                           <  1:274155/150‑1 (MQ=255)                                                                                               attaCCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGat                                                         >  2:964001/1‑150 (MQ=255)                                                                                                ttccgtGATCCCGTTTCCGC‑CCATCTATTGCGAATATCACGTAACCCAAT‑ACAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGat                                                         <  2:23980/144‑1 (MQ=255)                                                                                                 taCCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATAtt                                                      >  1:330642/1‑151 (MQ=255)                                                                                                     tGATGCCGATCCCGC‑CCAGCTATTGGGGTTAGCACGTAA‑CCCATAATAAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATtgtg                                                   <  2:624517/150‑1 (MQ=255)                                                                                                     tGATCCCGATTCCTC‑CCAGTTATTGGGAATATCACTTAATCCCAT‑AAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGata                                                        <  1:898302/145‑1 (MQ=255)                                                                                                          ccGATTCCGC‑CCACTTTTTGCGGATATCACGAAATCCCAT‑AAAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTgg                                                  <  1:292923/146‑1 (MQ=255)                                                                                                                ccGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGAt                                        >  1:75374/1‑150 (MQ=255)                                                                                                                ccGC‑CCAGCTTTTGCGGAATTCCCTTAACCCCAT‑AAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAg                                          <  2:132098/148‑1 (MQ=255)                                                                                                                 cGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGAt                                        >  1:717860/1‑149 (MQ=255)                                                                                                                  cgccaaattgttGGGGATATCAGGAATTCCTAT‑ATAAAGAATGTGCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATa                                       <  2:570754/140‑1 (MQ=255)                                                                                                                   t‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGt                                     >  1:221761/1‑150 (MQ=255)                                                                                                                     ccAGTTATTGCGGATATCACTTAACCCCAT‑ATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGt                                     <  1:2268/149‑1 (MQ=255)                                                                                                                             tGCGGATATCACGTACTCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCCTTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatagtc                           >  1:125773/1‑151 (MQ=255)                                                                                                                                       aCGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACg                 >  1:505671/1‑151 (MQ=255)                                                                                                                                       aCGTAATCCCAT‑ATCAAGAATATTACCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGTGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGCGATATTGTGGTTTGACAGCTAGTTGTCATAGTCATACCGCAc                  >  1:323130/1‑150 (MQ=255)                                                                                                                                                       aaGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTaa      <  2:254300/147‑1 (MQ=255)                                                                                                                                                       aaGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAATGac  >  2:991797/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                     CGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGT‑CCAGCTATTGCGGATATCACGTAATCCCAT‑ATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAATGAC  >  NC_000913/1575931‑1576228

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