Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,751,5370AC40.2% 126.0 / 20.7 82V111G (GTT→GGT) yiaTputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (35/14);  new base C (0/33);  total (35/47)
Fisher's exact test for biased strand distribution p-value = 7.07e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTC  >  NC_000913/3751388‑3751686
                                                                                                                                                         |                                                                                                                                                     
gaACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAAc                                                                                                                                                      >  1:509464/1‑151 (MQ=255)
 aaCGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAAcc                                                                                                                                                     >  2:1002822/1‑151 (MQ=255)
 aaCGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAAcc                                                                                                                                                     >  1:1261116/1‑151 (MQ=255)
  aCGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAAccc                                                                                                                                                    <  1:916050/151‑1 (MQ=255)
     cccAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAAccccc                                                                                                                                                  >  2:1098961/1‑150 (MQ=255)
     cccAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCg                                                                                                                                                 <  2:1099830/151‑1 (MQ=255)
               gtgtCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGt                                                                                                                                       >  2:398282/1‑151 (MQ=255)
                  tCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACt                                                                                                                                    >  1:1245778/1‑151 (MQ=255)
                   cAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACt                                                                                                                                    >  1:448768/1‑150 (MQ=255)
                    aGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTAt                                                                                                                                  >  1:959696/1‑151 (MQ=255)
                    aGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTAt                                                                                                                                  >  1:710224/1‑151 (MQ=255)
                     gCGACAGCCGCCCTATCTGCTTTTTGTGAAATACCGATACCCCCCCCCCCCAGCCTTT‑GCTGTTATCCAGAACGCCCGCAGCTGCAGAGGTTTT‑CCCGTTT‑CCCCACCG‑CTGGTGGTAATCCCAGGCACCCCCCGCCATAGCCGTACTAt                                                                                                                                  <  2:559396/150‑1 (MQ=255)
                      cGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATc                                                                                                                                 >  1:1023431/1‑150 (MQ=255)
                       gACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCAc                                                                                                                               >  1:873050/1‑151 (MQ=255)
                        aCAGACGACCTATCTGCATTTTGTGGAATCCCGATAGCTCCCCCCCCCAGCCGTT‑GCTGTTATCCAGACCGTCCGCAGCTGCAGAGGCTTT‑CACGTTT‑CCCCAACG‑CTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCAc                                                                                                                               <  1:914600/150‑1 (MQ=255)
                          aGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCt                                                                                                                            >  2:871656/1‑151 (MQ=255)
                          aGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCt                                                                                                                            >  2:584225/1‑151 (MQ=255)
                                ccTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCa                                                                                                                       >  2:90553/1‑150 (MQ=255)
                                ccTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCa                                                                                                                       >  1:921200/1‑150 (MQ=255)
                                catATCTCCATTTGGTGAAATACCGATAGCCCCCCCCCCCGCCCGTT‑GCTTTTACCCAGAACGCCCGCAGCTGCAGAGCCTTC‑CACGTTC‑CCCCAACG‑CCGGGGTTGACACAAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATc                                                                                                                        <  2:962266/147‑1 (MQ=255)
                                 cTATCTGATTTTGGTGGAATCCCGATAGCCCCCCCCCCCAGCCTTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGCCTTT‑CACTGTC‑CCCCACCC‑CTCGGGTTGATACCGGGCAACCCCCGCCATAGCCGTACTATCa                                                                                                                                <  2:892743/140‑1 (MQ=255)
                                    tcggctttttGTGGAATCCCGATAGCCCCCCCCCCCACCCTTT‑GCTGTTATCCAGAACGTCCGCACCGGCAGAGGTTTT‑CCCGTTT‑CCCCACCG‑CTCGTGGTCATCCCGGCCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCt                                                                                                                   <  1:1065382/144‑1 (MQ=255)
                                       gCATTTTGTGAAATCCCGATAGCCCCCCCCCCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGAGGAGGCTTT‑CCCGGTC‑CCCCACCG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCt                                                                                                                   <  1:414235/147‑1 (MQ=255)
                                          ttttGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCAt                                                                                                            >  1:944409/1‑151 (MQ=255)
                                           tttGTGGAATGCCGATACCCCCCCCCCCCAGCCGTT‑GCTGTAACCCAAATCCCCAGCAGTCGCGAAGTTTCC‑CCCGGTC‑CCCCACCG‑CCCGTGGGAACCCCGGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCAt                                                                                                            <  2:722118/150‑1 (MQ=255)
                                           tttGTGGAATACCGATAGCCCCCCCCCCCACCCGTT‑GCTTTTACCCAGAACGCCCGCAGCTGCAGAGGCTTT‑CCCGTTC‑CCCCAACC‑CTCGTGGTAATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTg                                                                                                               <  2:477042/147‑1 (MQ=255)
                                            ttGTGGAATACGGATATCTCCCCCACCCAGCCTTG‑GCTGTTTTCCAGACCGTCCCCAGCTGCAGATGCTTC‑CCCGTTC‑CCCAAACC‑CTGGTGGTATTCCTAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGc                                                                                                              <  2:1084594/147‑1 (MQ=255)
                                                   aTACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATTGCATAATc                                                                                                   >  2:638747/1‑151 (MQ=255)
                                                    tACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCg                                                                                                  >  2:824227/1‑151 (MQ=255)
                                                      ccGATAGCCCCCCCCCCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CCCGGTC‑CCCCACCG‑CTCGGGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTc                                                                                                <  1:682666/151‑1 (MQ=255)
                                                        gATAGCCCCCCCCCCCAGCCGTT‑GCTTTTACCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CCCGTGGGGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCAt                                                                                              <  1:1258555/151‑1 (MQ=255)
                                                          tacccccccccccCAGCCGTT‑GTTGTTATCCAGAACGCCCGCAGCTGCAGAGGTTTT‑CACGTTC‑CCCCACCG‑CTCGTGGTGATCCCAGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCAt                                                                                              <  1:629184/144‑1 (MQ=255)
                                                               cccccACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                       >  2:781722/1‑151 (MQ=255)
                                                                ccccccCAAGCCGTG‑GCTGTTATCAAAACGGCCCGCAGCGGGAGAGCCTTC‑CACGGTC‑CCCCAACG‑CTCGGGGTAATCCCGGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTa                                                                                      <  2:954032/151‑1 (MQ=255)
                                                                  cgccccaGACGTTGGCTTTGATCCAGACCCGCCGCGGCTG‑AGAGGGCTT‑CATGGTTCCCCCAACT‑CGCGTGGTGGTACCAGGAACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAg                                                                                     <  2:852093/148‑1 (MQ=255)
                                                                   cacccaGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAggg                                                                                   >  1:389983/1‑151 (MQ=255)
                                                                     cccaGCCTTT‑GCTTTTACCCAAACCGCCCGCAGCGGCAGAGTCTTT‑CACGGTT‑CCCCAACG‑CTCGGGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAggg                                                                                   <  2:880438/149‑1 (MQ=255)
                                                                       cacccttt‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGCCTTT‑CCCGTTT‑CCCCAACG‑CTCGTGGTAATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCaa                                                                               <  2:411740/145‑1 (MQ=255)
                                                                       caGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCa                                                                                >  2:82929/1‑150 (MQ=255)
                                                                             tt‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTcc                                                                          >  2:1864/1‑150 (MQ=255)
                                                                             tt‑GCTGTTATCCAGAACGCCCGCAGCTGAAGAGTCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGAACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTcc                                                                          <  1:412115/150‑1 (MQ=255)
                                                                                gCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGc                                                                       >  2:1150956/1‑151 (MQ=255)
                                                                                      aTCCAGAGCGGCCGCAGCTCGGGAGCCTTT‑CACGGTT‑CCCCATCG‑CTCGTGTGGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCg                                                                      <  2:412528/146‑1 (MQ=255)
                                                                                           gAACGCCCGCAGCTGCAGAGGCTTT‑CCCGGTT‑CCCCAACG‑CTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAgg                                                             <  1:133253/150‑1 (MQ=255)
                                                                                                 ccGCAGCTGCGGAGGCTTT‑CCCGTTC‑CCCCGACG‑CCCTGGGTAAAACCAGACAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGt                                                      <  1:689221/151‑1 (MQ=255)
                                                                                                 ccGCAGCTGCAGAGGTTTT‑CCCGTTT‑CCCCAACG‑TTCGGGGTGTTCCCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATg                                                          <  1:858390/147‑1 (MQ=255)
                                                                                                 ccGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGt                                                      >  1:1108314/1‑151 (MQ=255)
                                                                                                 ccGCAGCGCCAGGGCTTTT‑CCCGGTC‑CCCAACCG‑CTCGTGGTCAACCCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATg                                                          <  2:69716/147‑1 (MQ=255)
                                                                                                   gCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTc                                                     >  1:163786/1‑150 (MQ=255)
                                                                                                   gCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTc                                                     >  2:387063/1‑150 (MQ=255)
                                                                                                      gCTGCAGAGGCTTT‑CACGGTT‑CCCCACCG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAgg                                                  <  1:584242/150‑1 (MQ=255)
                                                                                                        tGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCt                                                >  2:623491/1‑150 (MQ=255)
                                                                                                         gCAAAGTCTTG‑CAGGGTC‑CCCCAACG‑CCCGTGGCAACACAAGCAAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGt                                                      <  1:1142218/143‑1 (MQ=255)
                                                                                                             gagCTTTCCACGTTT‑CCCCAACG‑CGCGGGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAg                                             <  2:285111/147‑1 (MQ=255)
                                                                                                                 ttt‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAttt                                      >  2:910610/1‑151 (MQ=255)
                                                                                                                   tccactttt‑CCCCAACG‑GTCGTGTTGTACCCGGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAt                                        <  1:6913/141‑1 (MQ=255)
                                                                                                                   t‑CACGGTT‑CCCCAACG‑CTGGGGGTGATACCGGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAtttt                                     <  2:1222746/150‑1 (MQ=255)
                                                                                                                     ccgggTT‑CCCCAACG‑CCGGGGGGGATCCCGGGAACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAtt                                       <  1:71832/144‑1 (MQ=255)
                                                                                                                     cACGGTT‑CCCCACCG‑CTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAttttt                                    <  1:678315/150‑1 (MQ=255)
                                                                                                                     cACGGTC‑CCCCAACG‑CTCGGGTTGACCCCGGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGTATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAtttttt                                   <  2:944362/151‑1 (MQ=255)
                                                                                                                       cGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTc                                 >  2:925064/1‑151 (MQ=255)
                                                                                                                        ggTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTc                                 >  1:210505/1‑150 (MQ=255)
                                                                                                                        ggTC‑CCCCAACG‑CTCGGGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAttttttt                                  <  2:853582/149‑1 (MQ=255)
                                                                                                                         gTT‑TCCCAACGCCTCGTGTTAACACCGGGCACCCACCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAttttttt                                  <  2:630997/149‑1 (MQ=255)
                                                                                                                          ttcccccACCC‑CTGGTGGTATTACCAGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAATCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAttttttt                                  <  2:957912/145‑1 (MQ=255)
                                                                                                                           t‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACtt                             >  1:39922/1‑151 (MQ=255)
                                                                                                                             ccccAGCG‑CTCGTGGTAATCCCAGGCACCCCCCGCCATAGCCGTACTACCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACttt                            <  1:1153339/151‑1 (MQ=255)
                                                                                                                             ccccAGCG‑CTCGTGGTAATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACtt                             <  2:1261050/150‑1 (MQ=255)
                                                                                                                             ccccAACG‑CTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACttt                            <  1:245869/151‑1 (MQ=255)
                                                                                                                             ccccAACG‑CCGGGGGTGATACCTGCCCCCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCa                                <  2:818098/147‑1 (MQ=255)
                                                                                                                                      ctggtggtgATACCCGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACa                      <  2:988896/146‑1 (MQ=255)
                                                                                                                                      cTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGAc                       >  2:173866/1‑148 (MQ=255)
                                                                                                                                         gtggtgACACCAGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATa                 <  2:1241946/151‑1 (MQ=255)
                                                                                                                                         gtggtaATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAg                     <  1:470300/147‑1 (MQ=255)
                                                                                                                                         gtggtaATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGa                    <  1:999936/148‑1 (MQ=255)
                                                                                                                                          ttgtgATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATa                 <  2:321099/148‑1 (MQ=255)
                                                                                                                                          tggggATACCAGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATa                 <  2:91466/150‑1 (MQ=255)
                                                                                                                                           ggtgATCCCAGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGaa                   <  2:921154/147‑1 (MQ=255)
                                                                                                                                              gATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAAc             <  1:1128792/150‑1 (MQ=255)
                                                                                                                                              gATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAAcc            >  1:650953/1‑151 (MQ=255)
                                                                                                                                                   cAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACg       >  1:1016568/1‑151 (MQ=255)
                                                                                                                                                    aGGCAACCCCCGCCATAGCCGTACTATCACAATTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGAAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACg       >  2:1219143/1‑150 (MQ=255)
                                                                                                                                                        aaCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTc  >  1:336555/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
GAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTT‑GCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTT‑CACGGTT‑CCCCAACG‑CTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTC  >  NC_000913/3751388‑3751686

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