Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,4940TG42.6% 93.5 / 21.2 62intergenic (‑512/+91)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (6/29);  new base G (26/0);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.37e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGC  >  NC_000913/987347‑987642
                                                                                                                                                    |                                                                                                                                                    
cttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGt                                                                                                                                                      >  2:235807/1‑148 (MQ=255)
  ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTg                                                                                                                                                    <  2:727174/148‑1 (MQ=255)
   tCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGc                                                                                                                                                   <  2:183425/148‑1 (MQ=255)
     aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATa                                                                                                                                              <  2:63407/151‑1 (MQ=255)
         aTCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGAtt                                                                                                                                          <  1:1136483/151‑1 (MQ=255)
         aTCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGaaa                                                                                                                                                >  1:525626/1‑145 (MQ=255)
         aTCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGCTGTaaaaaaaa                                                                                                                                           >  1:1089237/1‑145 (MQ=255)
           ccAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTg                                                                                                                                         <  2:276547/150‑1 (MQ=255)
             aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATGCtt                                                                                                                                      >  1:668418/1‑151 (MQ=255)
                  gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAGATTGtttttttt                                                                                                                                 >  1:610063/1‑146 (MQ=255)
                   ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTGAAAAAATTGCtttgattt                                                                                                                                >  1:420723/1‑145 (MQ=255)
                         aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaac                                                                                                                          <  1:843145/151‑1 (MQ=255)
                         aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCTAAAAAAATTCTTGACTTaaaaaa                                                                                                                           >  2:744704/1‑150 (MQ=255)
                            aGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGTGGTGAAAAAACACTGTTTGaactcaaaaa                                                                                                                           >  2:372462/1‑138 (MQ=255)
                                tGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAAATGCTTTTCTTattaaaaaca                                                                                                                       >  1:802722/1‑138 (MQ=255)
                                  aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACtt                                                                                                                 <  1:527640/151‑1 (MQ=255)
                                  aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACtt                                                                                                                 <  1:267215/151‑1 (MQ=255)
                                   aaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTg                                                                                                                <  2:31424/151‑1 (MQ=255)
                                   aaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGTAATAAATTGCTTGAATTTAAAAACAAATaaa                                                                                                                   >  2:8722/1‑147 (MQ=255)
                                       tgtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGATTAGTTTACTTAAAAAACAAAAAAAATGACa                                                                                                             >  1:506780/1‑150 (MQ=255)
                                         tgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGCAAAAAAATTCTTGACTTAACAAAAAGATAAATTGACAGa                                                                                                           >  2:1002628/1‑150 (MQ=255)
                                           cTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGACAAAAGATTGCATTTCTTAACAAACAAAAAAATTTACAGaaa                                                                                                         >  2:650630/1‑150 (MQ=255)
                                             aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAg                                                                                                       <  2:583844/150‑1 (MQ=255)
                                             aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGa                                                                                                      <  2:1127955/151‑1 (MQ=255)
                                              ttAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGGAAAAAAATTTGTTTAATTAAAAAAAAAAAAAAATTACAGaaaa                                                                                                        >  1:743584/1‑148 (MQ=255)
                                               tAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAt                                                                                                     <  1:509471/150‑1 (MQ=255)
                                               tAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAAAAGATTTCTTGACTTTAAAAAAAAATAAATTTACAGaaaa                                                                                                        >  1:18853/1‑147 (MQ=255)
                                                 gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaa                                                                                                  <  1:760588/151‑1 (MQ=255)
                                                 gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaa                                                                                                  <  1:593263/151‑1 (MQ=255)
                                                  cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTGCTTGACTTTACAAAAAGAAAAATTGACaaaaaaaaaaaa                                                                                                  >  2:88222/1‑150 (MQ=255)
                                                         aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGACGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc                                                                                          <  1:1079279/151‑1 (MQ=255)
                                                          aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc                                                                                         <  2:1091201/151‑1 (MQ=255)
                                                           aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCa                                                                                        <  2:670760/151‑1 (MQ=255)
                                                            aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAg                                                                                       <  1:565483/151‑1 (MQ=255)
                                                             cTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTCTTTACTTTACAAAAAAATAAAATTTAAGaaaaaaaaaaaaaaaa                                                                                            >  1:129664/1‑143 (MQ=255)
                                                                 tCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTGCTTTACTTTACAAAAAAAAAAATTTACAGAAAAAATAAAAAAAAAACCAGCggg                                                                                   >  2:427630/1‑150 (MQ=255)
                                                                  ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg                                                                                 <  2:940953/151‑1 (MQ=255)
                                                                  ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc                                                                                  <  2:523358/150‑1 (MQ=255)
                                                                        tAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAAATTCTTAGATTAAAAAACAAAAAAACTGACAAAAAAAATAAAAAAAGAAACAGCGGgcgggg                                                                              >  1:2601/1‑145 (MQ=255)
                                                                           gAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCCAAAAAATTGCTTTACTTTACAAAAAAAAAAAATTAAAGaaaaaaaaaaaaaaaaacccacggg                                                                                   >  2:275617/1‑129 (MQ=255)
                                                                             aTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATGCTTTACTTAAAAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAAACAGCGGGGGCCGGCTCtt                                                                        >  2:1068567/1‑149 (MQ=255)
                                                                               gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAAAAAATTGCTTTACTTAAAAAAAAAATAAATTTACAGAAAAGATAAAAAAAAAAACAGCGGGCGCgggggctt                                                                        >  2:1031108/1‑142 (MQ=255)
                                                                                 ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAc                                                                  <  1:889838/151‑1 (MQ=255)
                                                                                 ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTGGGGAAAAAATAGGTTTTCTTTACAAAAAAATAACATGAACGAAAAAAAAAAAAAAAAGACAACGGGGGCggggactt                                                                        >  1:452681/1‑140 (MQ=255)
                                                                                          tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAtt                                                          <  2:1140411/150‑1 (MQ=255)
                                                                                                      tGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTGCTTTACTTAACAAAAAAATAAATTGACAGAAAAAAAAAAAAAAAAGCCAGCGGGCGCTGGCGCTTTAAGAAGCCAAAATAAAAAATGGCggt                                             >  2:394306/1‑149 (MQ=255)
                                                                                                       gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATGATTTGATTAAAAAAAAAAAAAATTTATAGAAAAAAAAAAAAAAAAACAAGCGGGCGGTGGCTCTTTAAGAAGCAAAATTATAAGATgg                                                 >  2:977468/1‑146 (MQ=255)
                                                                                                          gCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgc                                         <  2:356520/151‑1 (MQ=255)
                                                                                                          gCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGCAATAAAATTGTTTACTTAAAAAAAAGATAAATTTACAGAAAAAATAAAAAAAAAAACAAGGGGGGCTGGCTCTTGAAGACGCTAAAATAAAAGATGGGGGTAc                                           >  2:64677/1‑149 (MQ=255)
                                                                                                           cAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgcg                                        <  2:922372/151‑1 (MQ=255)
                                                                                                                 aaTATTCCGTAACAGGATGATCGTTATGGCGGGGTTTAGAATATAATTCTTTACTTAAAAAAAACAGAACTTGATAAAAAAGAAAAAAAAAAAAACCACGGGGGCGCGGTTTTCGAGAAACTAAATTAGAAGCATGCGGTGcgcgcgaa                                    >  2:1040676/1‑145 (MQ=255)
                                                                                                                    atTCCGTAACAGGATGATCGTTATGGCGGGGTTGCTAAAAAATTGGTTTTATTTAAAAACAAAAAAATTTGAAGAAAAAATAAAAAAAAAACAAACGGGCGCTGGGTCTTTAAAACGCCAAATTAAAAAATGGGGGTACGCggggga                                   >  1:352647/1‑143 (MQ=255)
                                                                                                                    atTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTGCTTTTCTTTACAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCTGGCTCTTTAAGAAGCTAAAATAAAAAATGGGGTTTAGCGGAGGACGTg                               >  2:966934/1‑151 (MQ=255)
                                                                                                                     ttCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTTGTTGTCTTTACAAAAAAATAACATGGCAGAAAAAATAAAAAAAAAAACAGCGGGGGCTGGCTCTTTAAAACGCTAAAATATAAAGTGGGGGTACGCGGAGGACGTgg                              >  1:476707/1‑151 (MQ=255)
                                                                                                                      tCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                              <  2:43679/150‑1 (MQ=255)
                                                                                                                      tCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                              <  1:340677/150‑1 (MQ=255)
                                                                                                                         gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                          <  1:461541/151‑1 (MQ=255)
                                                                                                                          tAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTTTTTTACTTTAAAAACAAAAAAATTGAAAGAAAAAAAAAAAAAAAAAACAGCGGGGGGCGGGTCTTGTAAAAGCCAAAATAAAAAATGGGGGTTCGCGGAggaagcggg                             >  1:1061826/1‑139 (MQ=255)
                                                                                                                              aGGATGATCGTTATGGCGGGGGTGGAAATAAATTGGTTTACTTTAAAAAAAAAAAAATTGACAGAAAAAAAAAAAAAAAAACCCGCGGGGGCTGGCGCTTTAAGAAGCTAAAATAAAAGCTGGCGTTACGCCGAGGAAGTGGggaaggaa                      >  2:522403/1‑150 (MQ=255)
                                                                                                                              aGGATGATCGTTATGGCGGGGGTGCAAAAAAAATTCTTGACTTAACAAAAAAATAACTTGACAGAAAAAAAAAAAAAAAAACCAGCGGGGGCGGGCTCTTTTAGAAGCTAAATTAGAAGCTGGCGGTACGCGGAGGACGTggggaag                         >  2:933368/1‑143 (MQ=255)
                                                                                                                                       gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacg            <  1:961285/151‑1 (MQ=255)
                                                                                                                                                  ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc  <  2:666926/151‑1 (MQ=255)
                                                                                                                                                  ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc  >  2:434010/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGC  >  NC_000913/987347‑987642

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