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

TCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAC  >  NC_000913/987346‑987644
                                                                                                                                                    |                                                                                                                                                        
tcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGtt                                                                                                                                                         <  2:139051/150‑1 (MQ=255)
  tttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGt                                                                                                                                                          >  1:9124/1‑147 (MQ=255)
    tCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCaa                                                                                                                                                     <  1:1055097/150‑1 (MQ=255)
     cAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCaaa                                                                                                                                                    <  1:763953/150‑1 (MQ=255)
       gACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAg                                                                                                                                                 <  2:294644/151‑1 (MQ=255)
          aTCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGAtt                                                                                                                                              <  2:162060/151‑1 (MQ=255)
            ccAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTg                                                                                                                                             <  2:226688/150‑1 (MQ=255)
             cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAGATTGCt                                                                                                                                           >  2:441527/1‑151 (MQ=255)
              aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGaaaaaaaatgttt                                                                                                                                          >  1:463205/1‑141 (MQ=255)
              aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAATAAATTACtt                                                                                                                                          >  1:525157/1‑151 (MQ=255)
               gAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTg                                                                                                                                         <  1:1026491/151‑1 (MQ=255)
                 aTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAATAAATTGCTTgt                                                                                                                                        >  2:977735/1‑149 (MQ=255)
                    ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTa                                                                                                                                   <  1:949708/151‑1 (MQ=255)
                          aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaac                                                                                                                             <  1:367296/151‑1 (MQ=255)
                           aTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaaca                                                                                                                            <  1:436709/151‑1 (MQ=255)
                                aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATa                                                                                                                        <  2:6796/150‑1 (MQ=255)
                                   aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACtt                                                                                                                    <  1:715377/151‑1 (MQ=255)
                                    aaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTg                                                                                                                   <  1:616333/151‑1 (MQ=255)
                                           gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaaa                                                                                                            <  1:864674/151‑1 (MQ=255)
                                             tATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGAATGGTTTATTTTAAAAAAAAAAAACATGCCAGa                                                                                                              >  1:646029/1‑147 (MQ=255)
                                               ttAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGGAAATAGATTTTTTTTCTTTACAAAACAAAAAATTTACAGAAAAGa                                                                                                         >  1:944125/1‑150 (MQ=255)
                                                  gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTGGTTTATTTTACAAAAAAAAAAATTGACAGaaaaaaaaa                                                                                                      >  2:207778/1‑145 (MQ=255)
                                                   cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaa                                                                                                     <  2:716919/150‑1 (MQ=255)
                                                      aGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaaaa                                                                                                 <  1:163141/151‑1 (MQ=255)
                                                      aGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaaaa                                                                                                 <  2:962698/151‑1 (MQ=255)
                                                      aGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGAGTCTTGAAGAAAATGGTTTACTTTACAAAAAGAAAAATAGACCaaaaaaaaaaaaaa                                                                                                  >  1:707701/1‑150 (MQ=255)
                                                        aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTGATTTAGTTAACAAAAAAATAAATTTACAGAAAAGATaaaaaa                                                                                                  >  2:329945/1‑148 (MQ=255)
                                                            aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCa                                                                                          <  1:480607/151‑1 (MQ=255)
                                                            aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTTGTTTACTTTAAAAAAAAAAAAATTTACAGaaaaaaaaaaaaaaaa                                                                                               >  1:522612/1‑145 (MQ=255)
                                                            aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGCAAATAAAATTTTTTTCTTTAAAAAAAAATAAAATTACAGaaaaaaaaaaaaaaaa                                                                                               >  1:725436/1‑145 (MQ=255)
                                                                   ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc                                                                                    <  1:235253/150‑1 (MQ=255)
                                                                     gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGCGGTGTAAATAGATTGGTTGTACTAAAAAAAAGATAAATTCGGAGGAAAAATAAAAAAAAAGCCAGCggg                                                                                     >  2:566857/1‑147 (MQ=255)
                                                                         tAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGATTGGTTGACTTTACAAAAAGATAAATTTACAGAAAAAAAAAAAAAAAAGACCAGggggggggg                                                                                >  2:726083/1‑144 (MQ=255)
                                                                          aaGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTGTTTGGCTTTACAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAACCAACGGGGGCCgg                                                                               >  2:918658/1‑148 (MQ=255)
                                                                           aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGTAAAAAAAATTCTTGTATTTAAAAAAAAAAAAATTTAAAGAAAAAAAAAAAAAAAAAACAGCggggggcgg                                                                               >  1:379277/1‑141 (MQ=255)
                                                                            gAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTGCTTGACTTTACAAAAAAAAAAATTTACAAAAAAAAAAAAAAAAAAACCAGCGGGCGCCGGGTCtt                                                                          >  2:152817/1‑151 (MQ=255)
                                                                                gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAg                                                                      <  2:32202/151‑1 (MQ=255)
                                                                                  ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAAATGTTTTTATTAAAAAAAAAAAAAAATTTAAGAAAAGAAAAAAAAAAAAACAGCGGGCGCTGGCCCtt                                                                          >  1:375455/1‑145 (MQ=255)
                                                                                    aTTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGc                                                                  <  2:693124/151‑1 (MQ=255)
                                                                                     ttAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTTTTTGACTTAACAAAAAGAAAAATTTACAGAAAAAAAAAAAAAAAAACCAGCGGGGGGTGGGTCtt                                                                          >  2:1104265/1‑142 (MQ=255)
                                                                                       aGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAGATTTTTTTTCTTTAAAAAAAAAAAAAATGACAAAAAAAATAAAAAAAAAAACAACGGGGGCCGGCGCGTTAAGACGcca                                                                >  2:61508/1‑148 (MQ=255)
                                                                                         tttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAAATTCTTTGCTTTAAAAAAAGATAAAATTACAGAAAAAAAAAAAAAAAAACCAACGGGGGCTGGgtccttgagacagcaaa                                                              >  2:65201/1‑136 (MQ=255)
                                                                                           tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAtt                                                            <  1:730430/150‑1 (MQ=255)
                                                                                           tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAtt                                                            <  1:927817/150‑1 (MQ=255)
                                                                                           tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGGTAAAAAATTTGTTGAATTTAAAAAAAAATAAAATTAAAGAAAAAAAAAAAAAAAAAACACCGGGGGCGGGGTCTTGAAGAGGcccaa                                                              >  1:905147/1‑144 (MQ=255)
                                                                                            ttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTGCTTTGCTTAAAAAAAAAAAAAATTGACAGAAAAAAAAAAAAAAAAACCAGCGGGGGCTGGCTCTTGAAAAAGGTAAATTa                                                           >  2:1056723/1‑150 (MQ=255)
                                                                                            ttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGCTGGAAAAAGAATATTTTGCTTAAAAAAAAAAAAAATTTAACGAAAAAAAAAAAAAAAAAGCAGCGGGCGCCGGGGCTTTAAGAAGCCAAAATa                                                           >  1:992321/1‑150 (MQ=255)
                                                                                              gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGa                                                         <  2:673575/150‑1 (MQ=255)
                                                                                              gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGa                                                         <  2:1092586/150‑1 (MQ=255)
                                                                                              gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAAATGATTGTTTGTTTTATCAAAAAAAAAAATTTTCAGAAAAAAAAAAAAAAAAAACAACGGGGGCCGGGTCTTTAAGAAGCTAAAATAGa                                                         >  2:840110/1‑150 (MQ=255)
                                                                                              gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATGTTTTAGTTAACAAAAAAAAAAAATTAGAGAAAAAAAAAAAAAAAAAAAAACGGGCGGCGGGTCTTGAAGAAGCCaaaaaaaa                                                         >  1:84252/1‑145 (MQ=255)
                                                                                              gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGGTATGACGGGGTTGCAAATAGAATGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGa                                                         <  1:450726/150‑1 (MQ=255)
                                                                                               cgcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAGAATTCTTGACTTAACAAAAAAATAACTTTTCAGAAAAAATAAAAAAAAAGCCAGCGGGGGCTGGCTCTTTTAGAAGCTaaa                                                              >  2:417743/1‑144 (MQ=255)
                                                                                                 cgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTGAAATAAAATGCTTGAATTTAAAAAAAAAAAAAATTGGAGAAAAAAAAAAAAAAAAACCAGCGGGCGGCGGCCCCTGAAAAAGCGAAATTaaaaac                                                      >  2:888311/1‑148 (MQ=255)
                                                                                                   tAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAATAAATTGCTTGGCTTTACAAACAGAAAAATTGAAAGAAAAAAAAAAAAAAAAACCACCGGGCGCCGGGTCTTGGAGACGCCCAATTAGaa                                                        >  2:653427/1‑146 (MQ=255)
                                                                                                   tAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTTCTTTACTTTACAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAAACAACGGGCGCTGGGTCTTTTAAAAGCTAAATTaaaa                                                        >  2:1332/1‑146 (MQ=255)
                                                                                                    aaaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                  <  1:251059/151‑1 (MQ=255)
                                                                                                        gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTa                                              <  2:756677/151‑1 (MQ=255)
                                                                                                          tGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAc                                             <  2:45383/150‑1 (MQ=255)
                                                                                                               tGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGATTTCTTGAATTAACAAAAAAAAAAACTTACAGAAAAAATAAAAAAAAAGGCCACGGGCGCTGGCCCTTGTAGACGCCAAAATAAAAAATGGCGGTTcgcg                                          >  2:764826/1‑148 (MQ=255)
                                                                                                                gTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTCTTGTATTAAAAAAAAAATAAATTTGCAGAAAAAAAAAAAAAAAAACCAGCGGGGGCCGGGTCTTTGAGACGCCAAAATAAAAAATGGCGGTAcg                                            >  1:357989/1‑145 (MQ=255)
                                                                                                                 tAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTTCTTTTCTTTACAAAAAAAAAAATTTTCAGAAAAAAAAAAAAAAAAAACAGCGGGGGCCGGGGCTTTTAGACGCTAAAATAAAAACTGGCGGTAAGCgg                                         >  2:904853/1‑147 (MQ=255)
                                                                                                                     atTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGt                                  <  1:21714/150‑1 (MQ=255)
                                                                                                                       tCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                <  1:284052/150‑1 (MQ=255)
                                                                                                                        ccGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTGCTTGACTTTACAAAAAAAAAAAATTAAAGAAAAAATAAAAAAAAAACCAGCGGGCGCTGGCTCTTGAAGAAGCTAAAATAAAAGCTGGCGGTACGCGGAGGACGgggg                               >  2:285843/1‑150 (MQ=255)
                                                                                                                        ccGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGATTTCTTTACTTAAAAAAAAAAAAACTTTACAGAAAAAAAAAAAAAAAAAACAGCGGGGGCCGGCCCTTTAAAGAGCTAAAATAAAAAGTGGCGGTTAGCGGAggaaggggg                               >  2:489042/1‑142 (MQ=255)
                                                                                                                          gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                            <  1:660808/151‑1 (MQ=255)
                                                                                                                              cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                        <  2:257837/151‑1 (MQ=255)
                                                                                                                                           aTGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac           <  1:441518/151‑1 (MQ=255)
                                                                                                                                               cGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt       >  1:428225/1‑151 (MQ=255)
                                                                                                                                                  ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc    <  2:1097943/151‑1 (MQ=255)
                                                                                                                                                   gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc    >  2:848958/1‑150 (MQ=255)
                                                                                                                                                    ttGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc  >  1:697601/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                        
TCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAC  >  NC_000913/987346‑987644

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: