Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,644 C→T G95E (GGG→GAG)  thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,6440CT98.9% 322.5 / ‑4.7 93G95E (GGG→GAG) thiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (42/50);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

TAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGG  >  NC_000913/72508‑72773
                                                                                                                                        |                                                                                                                                 
tAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCa                                                                                                                                >  1:1274269/1‑140 (MQ=255)
tAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCa                                                                                                                                <  1:749853/140‑1 (MQ=255)
  aCGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGc                                                                                                                              <  2:1274269/140‑1 (MQ=255)
     ccACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCg                                                                                                                           <  1:1393770/140‑1 (MQ=255)
        cTCGCTGTCGCTGACCGCCGGAAAGCTCGCACGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCGGATg                                                                                                                        <  2:269261/140‑1 (MQ=255)
          cgctgtcgctgACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGtt                                                                                                                      <  2:1133691/140‑1 (MQ=255)
            ctgtcgctgACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCt                                                                                                                    >  2:851465/1‑140 (MQ=255)
            ccgtcgctgACCACCGGAAAGCGCCCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCt                                                                                                                    >  2:278921/3‑140 (MQ=255)
                cgctgACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGc                                                                                                                >  2:480322/1‑140 (MQ=255)
                 gctgACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCg                                                                                                               <  1:480322/140‑1 (MQ=255)
                  ctgACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGa                                                                                                              >  1:571010/1‑140 (MQ=255)
                  ctgACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGa                                                                                                              >  2:42208/1‑140 (MQ=255)
                      ccgccgGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGt                                                                                                          <  1:319187/140‑1 (MQ=255)
                        gccgGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCa                                                                                                        <  1:42208/140‑1 (MQ=255)
                         ccgGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCg                                                                                                           <  2:598606/136‑1 (MQ=255)
                         ccgGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCg                                                                                                           >  1:598606/1‑136 (MQ=255)
                         ccgGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAg                                                                                                       >  1:1386821/1‑140 (MQ=255)
                         ccgGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAg                                                                                                       >  2:135085/1‑140 (MQ=255)
                         ccgGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAg                                                                                                       >  1:762001/1‑140 (MQ=255)
                                 cTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGa                                                                                               >  2:1389359/1‑140 (MQ=255)
                                   cGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAAc                                                                                             <  1:638830/140‑1 (MQ=255)
                                    gCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACa                                                                                            >  1:713808/1‑140 (MQ=255)
                                          tAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGgttgt                                                                                      <  2:400597/140‑1 (MQ=255)
                                          tAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGgttgt                                                                                      <  1:774684/140‑1 (MQ=255)
                                          tAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGgttgt                                                                                      <  2:1386821/140‑1 (MQ=255)
                                          tAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGgttgt                                                                                      <  1:513891/140‑1 (MQ=255)
                                           aaCCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGgttgtt                                                                                     <  2:713808/140‑1 (MQ=255)
                                                cgcCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCt                                                                                >  2:1199098/1‑140 (MQ=255)
                                                cgcCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCt                                                                                >  2:992961/1‑140 (MQ=255)
                                                    aTTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGaaa                                                                            >  1:4223/1‑140 (MQ=255)
                                                     ttAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAAc                                                                           <  1:278921/140‑1 (MQ=255)
                                                        aaTTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGc                                                                        >  2:372372/1‑140 (MQ=255)
                                                          ttATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCAt                                                                      <  1:1202174/140‑1 (MQ=255)
                                                          ttATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCAt                                                                      <  2:762001/140‑1 (MQ=255)
                                                             tCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGa                                                                   >  1:530856/1‑140 (MQ=255)
                                                               aaTCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTAAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGaca                                                                 <  1:1059310/140‑1 (MQ=255)
                                                                     cATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGAc                                                           >  2:719322/1‑140 (MQ=255)
                                                                     cATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGAc                                                           >  2:747796/1‑140 (MQ=255)
                                                                           gcgggcgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTg                                                     >  2:1453161/1‑140 (MQ=255)
                                                                           gcgggcgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTg                                                     >  1:394473/1‑140 (MQ=255)
                                                                           gcgggcgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTg                                                     >  2:685363/1‑140 (MQ=255)
                                                                           gcgggcgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTg                                                     >  1:1191310/1‑140 (MQ=255)
                                                                             gggcgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGAc                                                   >  2:881417/1‑140 (MQ=255)
                                                                              ggcgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGAcg                                                  >  2:486762/1‑140 (MQ=255)
                                                                               gcgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGAcgg                                                 <  2:1412098/140‑1 (MQ=255)
                                                                               gcgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGAcgg                                                 <  2:530856/140‑1 (MQ=255)
                                                                               gcgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGAcgg                                                 <  1:965284/140‑1 (MQ=255)
                                                                               gcgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGAcgg                                                 <  1:1199098/140‑1 (MQ=255)
                                                                                cgATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGAcggc                                                >  2:1010041/1‑140 (MQ=255)
                                                                                 gATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGAcggcg                                               <  1:991644/140‑1 (MQ=255)
                                                                                    aGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCa                                            <  1:372372/140‑1 (MQ=255)
                                                                                    aGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCa                                            <  1:486762/140‑1 (MQ=255)
                                                                                    aGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCa                                            <  1:719322/140‑1 (MQ=255)
                                                                                    aGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCa                                            <  1:1389359/140‑1 (MQ=255)
                                                                                    aGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCa                                            <  1:135085/140‑1 (MQ=255)
                                                                                    aGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCa                                            <  1:747796/140‑1 (MQ=255)
                                                                                     gCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCa                                            <  2:593425/139‑1 (MQ=255)
                                                                                     gCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCa                                            >  1:593425/1‑139 (MQ=255)
                                                                                     gCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCAt                                           <  1:685363/140‑1 (MQ=255)
                                                                                      cGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATa                                          <  2:1474350/140‑1 (MQ=255)
                                                                                         gCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGACCAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGAc                                                   >  1:1128367/1‑128 (MQ=255)
                                                                                         gCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGAc                                                   <  2:1128367/128‑1 (MQ=255)
                                                                                            ttttCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTtgtg                                    <  1:992961/140‑1 (MQ=255)
                                                                                            ttttCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTtgtg                                    <  2:1191310/140‑1 (MQ=255)
                                                                                            ttttCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTtgtg                                    <  2:571010/140‑1 (MQ=255)
                                                                                              ttCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTtgtgtg                                  <  1:1453161/140‑1 (MQ=255)
                                                                                                ccccTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGAt                                <  2:4223/140‑1 (MQ=255)
                                                                                                 cccTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATc                               <  1:27974/140‑1 (MQ=255)
                                                                                                  ccTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTACCTCCAGCCCGATGTTCTGTGCGACCGTCAGGGGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGTCACCGGACGGCGTGACGGCGGCATGGTtgtgtcaggt                              <  2:1470996/140‑6 (MQ=255)
                                                                                                     gctgTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACg                           <  2:1248287/140‑1 (MQ=255)
                                                                                                      ctgTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCa                                            <  1:404943/122‑1 (MQ=255)
                                                                                                      ctgTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCa                                            >  2:404943/1‑122 (MQ=255)
                                                                                                      ctgTACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGc                          <  1:881417/140‑1 (MQ=255)
                                                                                                         tACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTtgt                                     <  2:264110/126‑1 (MQ=255)
                                                                                                         tACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTtgt                                     <  2:1291746/126‑1 (MQ=255)
                                                                                                         tACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTtgt                                     >  1:264110/1‑126 (MQ=255)
                                                                                                         tACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCAt                       >  1:1076001/1‑140 (MQ=255)
                                                                                                         tACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCAt                       >  1:191819/1‑140 (MQ=255)
                                                                                                         tACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCAt                       <  1:851465/140‑1 (MQ=255)
                                                                                                         tACCGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCGGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTtgt                                     >  1:1291746/1‑126 (MQ=255)
                                                                                                           ccGCGTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCg                     >  2:915444/1‑140 (MQ=255)
                                                                                                              cgTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTAc                            >  2:712270/1‑130 (MQ=255)
                                                                                                              cgTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTAc                            <  1:712270/130‑1 (MQ=255)
                                                                                                              cgTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATa                  >  2:13474/1‑140 (MQ=255)
                                                                                                              cgTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATa                  >  1:1200576/1‑140 (MQ=255)
                                                                                                              cgTTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATa                  >  2:1265047/1‑140 (MQ=255)
                                                                                                                ttCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGt                >  1:1326816/1‑140 (MQ=255)
                                                                                                                ttCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGt                >  1:827050/1‑140 (MQ=255)
                                                                                                                ttCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGt                >  2:1289471/1‑140 (MQ=255)
                                                                                                                       tCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGAt                   <  2:456049/130‑1 (MQ=255)
                                                                                                                       tCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGAt                   >  1:456049/1‑130 (MQ=255)
                                                                                                                          aTCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAg             >  2:1382041/1‑133 (MQ=255)
                                                                                                                          aTCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAg             <  1:1382041/133‑1 (MQ=255)
                                                                                                                          aTCCCGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCg      >  1:255238/1‑140 (MQ=255)
                                                                                                                              cGGATTTAGCTCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCgctgg  <  2:191819/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                 
TAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGG  >  NC_000913/72508‑72773

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: