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

TCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGT  >  NC_000913/987350‑987637
                                                                                                                                                 |                                                                                                                                               
tCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCaa                                                                                                                                            <  1:25513/150‑1 (MQ=255)
tCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTaa                                                                                                                                            >  2:110560/1‑150 (MQ=255)
 cAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGaaa                                                                                                                                           >  2:65738/1‑150 (MQ=255)
  aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATa                                                                                                                                         <  1:321586/151‑1 (MQ=255)
  aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTggta                                                                                                                                            >  1:46529/1‑145 (MQ=255)
          aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTGCtt                                                                                                                                 >  2:878695/1‑151 (MQ=255)
           gAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTg                                                                                                                                <  1:421525/151‑1 (MQ=255)
               gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACtt                                                                                                                            <  1:822053/151‑1 (MQ=255)
               gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACt                                                                                                                             <  1:368285/150‑1 (MQ=255)
                  aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaac                                                                                                                         <  2:372837/151‑1 (MQ=255)
                    aaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaa                                                                                                                       <  2:773206/151‑1 (MQ=255)
                     aaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAATTTCTTTACTTtacaa                                                                                                                       >  1:449880/1‑150 (MQ=255)
                      aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTTTTTGACTTtaaaaa                                                                                                                      >  2:400196/1‑150 (MQ=255)
                            aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATa                                                                                                                <  2:592915/150‑1 (MQ=255)
                             ggAAAAGGGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAAc                                                                                                              <  1:56706/150‑1 (MQ=255)
                               aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTTTTTTAATTTAAAAACAAATAAAtt                                                                                                            >  1:458091/1‑151 (MQ=255)
                               aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTGCTTGACTTaaaaaaaaaaaaaatt                                                                                                            >  2:346536/1‑141 (MQ=255)
                                       gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaaa                                                                                                    <  1:465429/151‑1 (MQ=255)
                                       gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaaa                                                                                                    <  1:864491/151‑1 (MQ=255)
                                       gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATGTTTTTCTTAACAAAAAAATAAAATGACAGaa                                                                                                     >  1:48937/1‑150 (MQ=255)
                                        cTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCCAAAAGATTGCTTGACTTTAAAAAAAGATAAATTTACAGaaa                                                                                                    >  2:60333/1‑150 (MQ=255)
                                          aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGa                                                                                                 <  1:677271/151‑1 (MQ=255)
                                          aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTGTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATGGATTGCTGGACTTAACAAACAGAGGTCTTGACAGAAAAg                                                                                                  <  1:270760/150‑1 (MQ=255)
                                                 aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaa                                                                                           <  1:985120/150‑1 (MQ=255)
                                                 aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaa                                                                                            <  1:715132/149‑1 (MQ=255)
                                                 aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaa                                                                                            <  2:41612/149‑1 (MQ=255)
                                                     aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAagag                                                                                      <  2:46535/151‑1 (MQ=255)
                                                      aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc                                                                                     <  1:747142/151‑1 (MQ=255)
                                                      aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc                                                                                     <  1:40374/151‑1 (MQ=255)
                                                       aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc                                                                                    <  1:925073/151‑1 (MQ=255)
                                                       aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc                                                                                    <  2:235685/151‑1 (MQ=255)
                                                       aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGGTCATAGATTGCTTTAATTTAAAAAAAAAAAAATTCGAAGAAAAAAGAAAAAAAAAGcc                                                                                    >  2:658005/1‑151 (MQ=255)
                                                        aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCa                                                                                   <  1:268566/151‑1 (MQ=255)
                                                        aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTGAAAAAAATTGCTTGACTTTAAAAAAAAAAAAATTTACAGAAAAAATaaaaaaaaaaac                                                                                    >  2:561979/1‑145 (MQ=255)
                                                            aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                               <  1:461665/151‑1 (MQ=255)
                                                             aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAATTTCTTGTCTTTACAAAAAAAAAAATTTAAAGaaaaaaaaaaaaaaaaagcagaggg                                                                              >  1:30563/1‑140 (MQ=255)
                                                                 gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                          <  2:421267/151‑1 (MQ=255)
                                                                  cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGATTGCTTTACTTAAAAAACAAATAAATTGAAAGAAAAAATAAAAAAAAAACCAGCGGgcgccg                                                                         >  2:981811/1‑149 (MQ=255)
                                                                       aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTGAAAAAAAATGTTTTACTTAAAAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAGACAGCGGGGGCCGGctc                                                                     >  1:923280/1‑150 (MQ=255)
                                                                         aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTGTTTTCTTTAAAAAAAAAAAAATTGTCAGAAAAAAAAAAAAAAAAAACAGAGGGCGCTGGctct                                                                    >  1:841309/1‑149 (MQ=255)
                                                                            gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTGCTTTACTTAACAAAAAGATAAATTGACAGAAAAGAAAAAAAAAAAACCAGCGGGGGCCGGCTCtt                                                                   >  2:691460/1‑147 (MQ=255)
                                                                            gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAAATTCTTTACTTTACAAAAAAAAAAAATGACAGAAAAAAAAAAAAAAAAGACAGCGGGGGCCGGCGCtt                                                                   >  1:483298/1‑147 (MQ=255)
                                                                              ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAGATTGTTTTGCTTTACAAAAAAATAAATTGACAGAAAAAATAAAAAAAAAGCCAGCGGGGGCGGGGTCTTGAAGa                                                              >  2:814043/1‑150 (MQ=255)
                                                                                  tAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCt                                                          <  2:510734/150‑1 (MQ=255)
                                                                                    gTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAATTTTTTTACTTTAAAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCTGGGGCTTGTAGACGCCaa                                                        >  1:931394/1‑150 (MQ=255)
                                                                                      ttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGTCAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAt                                                      <  2:844119/150‑1 (MQ=255)
                                                                                      ttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAACAGAGTGGTTTTCTTCAAAAACAAAAAACATTTCAGAAAAGATAAAAAAAAAAACAGAGGGCGGCGGGCCTTGTAGACGCCaaa                                                       >  2:933294/1‑149 (MQ=255)
                                                                                      ttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGAATGTTTTTATTTAAAAAAAGATAAAATTTCAGGAAAAAAAAAAAAAAAGGCAACGGGGGCCGGCGCCTTAAGaagacaaaat                                                     >  1:349064/1‑142 (MQ=255)
                                                                                          gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGaa                                                 <  1:878595/151‑1 (MQ=255)
                                                                                           cgcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAGAATTCTTTACTTAACAAAAAAATAAATTGACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCCGGGTCCTGAAGACGACAAATTa                                                    >  2:275837/1‑147 (MQ=255)
                                                                                            gcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGc                                               <  2:209407/151‑1 (MQ=255)
                                                                                             cgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTGTTTTTATTAAAAAAAAAATAAATTGATAGAAAAAATAAAAAAAAAAACAGCGGGCGGCGGCTCTTGAAGAAGCtaaaataaaaaat                                              >  1:797961/1‑148 (MQ=255)
                                                                                              gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCCAAATTAGAAGCTg                                             <  1:732957/151‑1 (MQ=255)
                                                                                              gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTTTTTAATTTAAAAAAAAAAAAATTTAAAGAAAAAAAAAAAAAAAAGACAACGGGCGCCGGGTCTTTGAGAAGCCAAAATaaaa                                                 >  1:555547/1‑147 (MQ=255)
                                                                                                       gCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTGCTTGAATTAAAAAAAAGATAAATTTACAGAAAAAAAAAAAAAAAAACCAACGGGGGCTGGCTCTTGAAGAAGCTAAATTAGaaaatgggggttcg                                     >  1:614/1‑138 (MQ=255)
                                                                                                         aaTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                  <  2:339817/151‑1 (MQ=255)
                                                                                                         aaTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGATTTCTTGTCTTAAAAAAAAAATAACTTTACAGAAAAGATAAAAAAAAAACCAACGGGGGCTGGCTCTTGAAGACGCTAAATTAAAAGCTGGCGTTAcgc                                    >  2:671685/1‑149 (MQ=255)
                                                                                                              aaTATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTAc                             <  2:220375/151‑1 (MQ=255)
                                                                                                                tatTCCGTAACAGGATGATCGTTATGGCGGGGGTGGCAAAAAAATTGTTTTCTTTAAAAAAAAATAAATTGAAAGAAAAAAAAAAAAAAAAACCAACGGGGGCCGGCCCCTTTAGACGCTAAAATAAAAACTGGCGGTGCGCCGAgga                              >  1:600485/1‑146 (MQ=255)
                                                                                                                    ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                        <  1:368356/150‑1 (MQ=255)
                                                                                                                     cGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAGAATTCTTTGATTAACAAAAAAATAAATTGACAGAAAAGAAAAAAAAAAAGCCAGCGGGGGCGGGCCCTTGTAGACGCCAAAATAAAAACTGGCGATACGAGGAGTACGgggg                        >  1:585360/1‑149 (MQ=255)
                                                                                                                     cGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTGGTTTACTTTACAAAAAGATAAATTGACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCTGGCTCTTTAAAAAGCTAAATTAAAAAATGGCGGTAAGCGGGGGACGgggg                        >  1:98931/1‑149 (MQ=255)
                                                                                                                       tAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATTATTGTCTTTACAAAAAAAAAAACTTACAAAAAAAAAAAAAAAAAAAACAGCGGGCGCTGGGTCTTTTAGAGGCCAAATTATAAACAGGCGtttaccggaggac                             >  2:444816/1‑131 (MQ=255)
                                                                                                                          cAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTCTTGGCTTTACAAAAAAAAAAAATTACAGAAAAAAAAAAAAAAAAAGCAACGGGGGCCGGCTCTTTAAGAAGCCAAAATAAAAGGTGGGGGTACGCGGAGGACGgggggaagg                   >  2:237315/1‑144 (MQ=255)
                                                                                                                                      tATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgt      <  2:978284/150‑1 (MQ=255)
                                                                                                                                         ggCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacgt  >  2:910383/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGT  >  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: