Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,4940TG48.5% 100.0 / 19.0 66intergenic (‑512/+91)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (2/32);  new base G (31/1);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 4.98e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGT  >  NC_000913/987350‑987637
                                                                                                                                                 |                                                                                                                                               
tCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCaa                                                                                                                                            <  1:25513/150‑1 (MQ=255)
tCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGTaa                                                                                                                                            >  2:110560/1‑150 (MQ=255)
 cAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGaaa                                                                                                                                           >  2:65738/1‑150 (MQ=255)
  aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATa                                                                                                                                         <  1:321586/151‑1 (MQ=255)
  aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTggta                                                                                                                                            >  1:46529/1‑145 (MQ=255)
          aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAATAAATTGCtt                                                                                                                                 >  2:878695/1‑151 (MQ=255)
           gAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTg                                                                                                                                <  1:421525/151‑1 (MQ=255)
               gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACtt                                                                                                                            <  1:822053/151‑1 (MQ=255)
               gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACt                                                                                                                             <  1:368285/150‑1 (MQ=255)
                  aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTaac                                                                                                                         <  2:372837/151‑1 (MQ=255)
                    aaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTaacaa                                                                                                                       <  2:773206/151‑1 (MQ=255)
                     aaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGTAAAAAAATTTCTTTACTTtacaa                                                                                                                       >  1:449880/1‑150 (MQ=255)
                      aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAATTTTTTGACTTtaaaaa                                                                                                                      >  2:400196/1‑150 (MQ=255)
                            aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATa                                                                                                                <  2:592915/150‑1 (MQ=255)
                             ggAAAAGGGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAAc                                                                                                              <  1:56706/150‑1 (MQ=255)
                               aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAATTTTTTTAATTTAAAAACAAATAAAtt                                                                                                            >  1:458091/1‑151 (MQ=255)
                               aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAATTGCTTGACTTaaaaaaaaaaaaaatt                                                                                                            >  2:346536/1‑141 (MQ=255)
                                       gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaaa                                                                                                    <  1:465429/151‑1 (MQ=255)
                                       gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaaa                                                                                                    <  1:864491/151‑1 (MQ=255)
                                       gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAAAATGTTTTTCTTAACAAAAAAATAAAATGACAGaa                                                                                                     >  1:48937/1‑150 (MQ=255)
                                        cTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCCAAAAGATTGCTTGACTTTAAAAAAAGATAAATTTACAGaaa                                                                                                    >  2:60333/1‑150 (MQ=255)
                                          aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGa                                                                                                 <  1:677271/151‑1 (MQ=255)
                                          aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTGTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAATGGATTGCTGGACTTAACAAACAGAGGTCTTGACAGAAAAg                                                                                                  <  1:270760/150‑1 (MQ=255)
                                                 aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaa                                                                                           <  1:985120/150‑1 (MQ=255)
                                                 aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaa                                                                                            <  1:715132/149‑1 (MQ=255)
                                                 aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaa                                                                                            <  2:41612/149‑1 (MQ=255)
                                                     aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAagag                                                                                      <  2:46535/151‑1 (MQ=255)
                                                      aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc                                                                                     <  1:747142/151‑1 (MQ=255)
                                                      aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc                                                                                     <  1:40374/151‑1 (MQ=255)
                                                       aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc                                                                                    <  1:925073/151‑1 (MQ=255)
                                                       aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc                                                                                    <  2:235685/151‑1 (MQ=255)
                                                       aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGGGGTCATAGATTGCTTTAATTTAAAAAAAAAAAAATTCGAAGAAAAAAGAAAAAAAAAGcc                                                                                    >  2:658005/1‑151 (MQ=255)
                                                        aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCa                                                                                   <  1:268566/151‑1 (MQ=255)
                                                        aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTTGAAAAAAATTGCTTGACTTTAAAAAAAAAAAAATTTACAGAAAAAATaaaaaaaaaaac                                                                                    >  2:561979/1‑145 (MQ=255)
                                                            aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                               <  1:461665/151‑1 (MQ=255)
                                                             aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGTAAAAAAATTTCTTGTCTTTACAAAAAAAAAAATTTAAAGaaaaaaaaaaaaaaaaagcagaggg                                                                              >  1:30563/1‑140 (MQ=255)
                                                                 gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                          <  2:421267/151‑1 (MQ=255)
                                                                  cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAGATTGCTTTACTTAAAAAACAAATAAATTGAAAGAAAAAATAAAAAAAAAACCAGCGGgcgccg                                                                         >  2:981811/1‑149 (MQ=255)
                                                                       aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTTGAAAAAAAATGTTTTACTTAAAAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAGACAGCGGGGGCCGGctc                                                                     >  1:923280/1‑150 (MQ=255)
                                                                         aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAAATTGTTTTCTTTAAAAAAAAAAAAATTGTCAGAAAAAAAAAAAAAAAAAACAGAGGGCGCTGGctct                                                                    >  1:841309/1‑149 (MQ=255)
                                                                            gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAATAAATTGCTTTACTTAACAAAAAGATAAATTGACAGAAAAGAAAAAAAAAAAACCAGCGGGGGCCGGCTCtt                                                                   >  2:691460/1‑147 (MQ=255)
                                                                            gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAATAAAATTCTTTACTTTACAAAAAAAAAAAATGACAGAAAAAAAAAAAAAAAAGACAGCGGGGGCCGGCGCtt                                                                   >  1:483298/1‑147 (MQ=255)
                                                                              ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAATAGATTGTTTTGCTTTACAAAAAAATAAATTGACAGAAAAAATAAAAAAAAAGCCAGCGGGGGCGGGGTCTTGAAGa                                                              >  2:814043/1‑150 (MQ=255)
                                                                                  tAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCt                                                          <  2:510734/150‑1 (MQ=255)
                                                                                    gTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGTAAAAAAATTTTTTTACTTTAAAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCTGGGGCTTGTAGACGCCaa                                                        >  1:931394/1‑150 (MQ=255)
                                                                                      ttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGTCAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAt                                                      <  2:844119/150‑1 (MQ=255)
                                                                                      ttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAACAGAGTGGTTTTCTTCAAAAACAAAAAACATTTCAGAAAAGATAAAAAAAAAAACAGAGGGCGGCGGGCCTTGTAGACGCCaaa                                                       >  2:933294/1‑149 (MQ=255)
                                                                                      ttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAGAATGTTTTTATTTAAAAAAAGATAAAATTTCAGGAAAAAAAAAAAAAAAGGCAACGGGGGCCGGCGCCTTAAGaagacaaaat                                                     >  1:349064/1‑142 (MQ=255)
                                                                                          gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGaa                                                 <  1:878595/151‑1 (MQ=255)
                                                                                           cgcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAATAGAATTCTTTACTTAACAAAAAAATAAATTGACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCCGGGTCCTGAAGACGACAAATTa                                                    >  2:275837/1‑147 (MQ=255)
                                                                                            gcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGc                                               <  2:209407/151‑1 (MQ=255)
                                                                                             cgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAATTGTTTTTATTAAAAAAAAAATAAATTGATAGAAAAAATAAAAAAAAAAACAGCGGGCGGCGGCTCTTGAAGAAGCtaaaataaaaaat                                              >  1:797961/1‑148 (MQ=255)
                                                                                              gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCCAAATTAGAAGCTg                                             <  1:732957/151‑1 (MQ=255)
                                                                                              gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAAATTTTTTTAATTTAAAAAAAAAAAAATTTAAAGAAAAAAAAAAAAAAAAGACAACGGGCGCCGGGTCTTTGAGAAGCCAAAATaaaa                                                 >  1:555547/1‑147 (MQ=255)
                                                                                                       gCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAAATTGCTTGAATTAAAAAAAAGATAAATTTACAGAAAAAAAAAAAAAAAAACCAACGGGGGCTGGCTCTTGAAGAAGCTAAATTAGaaaatgggggttcg                                     >  1:614/1‑138 (MQ=255)
                                                                                                         aaTGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                  <  2:339817/151‑1 (MQ=255)
                                                                                                         aaTGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAGATTTCTTGTCTTAAAAAAAAAATAACTTTACAGAAAAGATAAAAAAAAAACCAACGGGGGCTGGCTCTTGAAGACGCTAAATTAAAAGCTGGCGTTAcgc                                    >  2:671685/1‑149 (MQ=255)
                                                                                                              aaTATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTAc                             <  2:220375/151‑1 (MQ=255)
                                                                                                                tatTCCGTAACAGGATGATCGTTATGGCGGGGGTGGCAAAAAAATTG‑TTTTCTTTAAAAAAAAATAAATTGAAAGAAAAAAAAAAAAAAAAACCAACGGGGGCCGGCCCCTTTAGACGCTAAAATAAAAACTGGCGGTGCGCCGAgga                              >  1:600485/1‑146 (MQ=255)
                                                                                                                    ccGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                        <  1:368356/150‑1 (MQ=255)
                                                                                                                     cGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAATAGAATTCTTTGATTAACAAAAAAATAAATTGACAGAAAAGAAAAAAAAAAAGCCAGCGGGGGCGGGCCCTTGTAGACGCCAAAATAAAAACTGGCGATACGAGGAGTACGgggg                        >  1:585360/1‑149 (MQ=255)
                                                                                                                     cGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAAATTGGTTTACTTTACAAAAAGATAAATTGACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCTGGCTCTTTAAAAAGCTAAATTAAAAAATGGCGGTAAGCGGGGGACGgggg                        >  1:98931/1‑149 (MQ=255)
                                                                                                                       tAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAAAATTATTGTCTTTACAAAAAAAAAAACTTACAAAAAAAAAAAAAAAAAAAACAGCGGGCGCTGGGTCTTTTAGAGGCCAAATTATAAACAGGCGtttaccggaggac                             >  2:444816/1‑131 (MQ=255)
                                                                                                                          cAGGATGATCGTTATGGC‑GGGGGTGCAAAAAAATTTCTTGGCTTTACAAAAAAAAAAAATTACAGAAAAAAAAAAAAAAAAAGCAACGGGGGCCGGCTCTTTAAGAAGCCAAAATAAAAGGTGGGGGTACGCGGAGGACGgggggaagg                   >  2:237315/1‑144 (MQ=255)
                                                                                                                                      tATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgt      <  2:978284/150‑1 (MQ=255)
                                                                                                                                         ggC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacgt  >  2:910383/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGT  >  NC_000913/987350‑987637

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