Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,161,966 C→A L66M (CTG→ATG)  yijD → DUF1422 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,9660CA98.8% 320.6 / NA 83L66M (CTG→ATG) yijDDUF1422 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (48/34);  total (49/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGC  >  NC_000913/4161819‑4162113
                                                                                                                                                   |                                                                                                                                                   
gCTGGCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATgg                                                                                                                                                  >  1:1196350/1‑151 (MQ=255)
 cTGGCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGc                                                                                                                                                 >  2:974227/1‑151 (MQ=255)
  tGGCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGc                                                                                                                                                 >  2:495496/1‑150 (MQ=255)
  tGGCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGc                                                                                                                                                 >  1:204915/1‑150 (MQ=255)
   ggCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGc                                                                                                                                                 <  1:270899/149‑1 (MQ=255)
   ggCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAg                                                                                                                                               <  1:263563/151‑1 (MQ=255)
    gCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGc                                                                                                                                              >  2:642292/1‑151 (MQ=255)
       tATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGc                                                                                                                                           >  1:1064378/1‑151 (MQ=255)
           gATTAATGGTACTTCCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATGATCGTTATCTTAATCGCACTATGCCGTTAGGCTTGCCGGGTATGGCTGCTGCCTg                                                                                                                                        >  1:33532/1‑150 (MQ=255)
           gATTAATGGCACTTTCGCAGCGCTGTTTAGCACCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGCGCTGACGGTTTACTGCCTGCCTCAACGTTATTTTAATCGCACTCTGCCGGTAGGATTGCCGGGGCTGGCCGCTGCCt                                                                                                                                         >  1:698393/1‑149 (MQ=255)
              tAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGttt                                                                                                                                     <  2:1036003/150‑1 (MQ=255)
              tAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGCCTTGCCGGGTATGGCAGCTGCCTGtttt                                                                                                                                    >  1:394940/1‑151 (MQ=255)
               aaTGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGtttt                                                                                                                                    >  1:613474/1‑150 (MQ=255)
                 tGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTAtt                                                                                                                                 >  2:468613/1‑151 (MQ=255)
                   gTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATtc                                                                                                                                >  1:1310160/1‑150 (MQ=255)
                    tACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATtctc                                                                                                                              >  2:1096137/1‑151 (MQ=255)
                       tttCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGc                                                                                                                           >  2:1099009/1‑151 (MQ=255)
                             aGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTActg                                                                                                                     <  2:1058664/151‑1 (MQ=255)
                              gcgcTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTActgc                                                                                                                    >  2:823066/1‑151 (MQ=255)
                               cgcTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTActgct                                                                                                                   >  2:286705/1‑151 (MQ=255)
                               cgcTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTActgct                                                                                                                   <  1:647687/151‑1 (MQ=255)
                                gcTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTActgc                                                                                                                    >  1:589000/1‑149 (MQ=255)
                                  tGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGt                                                                                                                 >  2:162094/1‑150 (MQ=255)
                                           ccATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGt                                                                                                       >  2:247514/1‑151 (MQ=255)
                                           ccATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGt                                                                                                       >  2:1097258/1‑151 (MQ=255)
                                             aTTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAg                                                                                                     <  1:185228/151‑1 (MQ=255)
                                              ttGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGt                                                                                                    >  2:1087714/1‑151 (MQ=255)
                                               tgtgCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGt                                                                                                    <  1:264877/150‑1 (MQ=255)
                                                gtgCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTc                                                                                                  >  2:672845/1‑151 (MQ=255)
                                                  gCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGt                                                                                                <  2:927587/151‑1 (MQ=255)
                                                  gCCATTTTCTGAATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTTCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCg                                                                                                 >  2:815737/1‑150 (MQ=255)
                                                        ttCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCATGCTTTATTCTCGGCGTACTGCTTTACACCACGGTAGTTCGTGCGGa                                                                                           >  1:4488/1‑150 (MQ=255)
                                                            gTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAata                                                                                        <  2:455412/149‑1 (MQ=255)
                                                            gTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAata                                                                                        <  2:1407/149‑1 (MQ=255)
                                                              aTTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATcc                                                                                     <  2:133324/150‑1 (MQ=255)
                                                              aTTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCg                                                                                    >  2:822625/1‑151 (MQ=255)
                                                                tCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCgg                                                                                   <  2:174987/150‑1 (MQ=255)
                                                                        attTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGct                                                                          >  1:660799/1‑151 (MQ=255)
                                                                         tttCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGctc                                                                         >  1:1082868/1‑151 (MQ=255)
                                                                         tttCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGctc                                                                         >  2:229515/1‑151 (MQ=255)
                                                                          ttCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGctc                                                                         <  1:487125/150‑1 (MQ=255)
                                                                            cccTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTaa                                                                      <  2:591782/151‑1 (MQ=255)
                                                                             ccTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAAc                                                                     >  1:326044/1‑151 (MQ=255)
                                                                             ccTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGCACTGCTGTACACCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCGaa                                                                      >  1:1019609/1‑150 (MQ=255)
                                                                                ggTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAActtc                                                                  <  1:1239397/151‑1 (MQ=255)
                                                                                 gTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAActtct                                                                 <  2:819753/151‑1 (MQ=255)
                                                                                   gCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAActtctt                                                                <  2:1305394/150‑1 (MQ=255)
                                                                                   gCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACttcttc                                                               <  1:653611/151‑1 (MQ=255)
                                                                                     tGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTccc                                                             >  1:546883/1‑151 (MQ=255)
                                                                                      gACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCg                                                            <  1:472228/151‑1 (MQ=255)
                                                                                        cGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGc                                                          >  2:1003124/1‑151 (MQ=255)
                                                                                         ggTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGc                                                          >  2:1002399/1‑150 (MQ=255)
                                                                                              aCTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTAct                                                    <  2:1009656/151‑1 (MQ=255)
                                                                                                tgcctgcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTTTTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTActc                                                   >  1:46383/1‑150 (MQ=255)
                                                                                                  cctgcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTcc                                                <  2:310847/151‑1 (MQ=255)
                                                                                                  cctgcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTcc                                                <  1:1163394/151‑1 (MQ=255)
                                                                                                    tgcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCg                                               >  1:87995/1‑150 (MQ=255)
                                                                                                      cATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTc                                             >  2:758240/1‑150 (MQ=255)
                                                                                                      cATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCa                                            <  2:24494/151‑1 (MQ=255)
                                                                                                        tCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCAtt                                          >  2:670084/1‑151 (MQ=255)
                                                                                                        tCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACCGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCAtt                                          >  2:679016/1‑151 (MQ=255)
                                                                                                           aCGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATg                                       <  1:91064/151‑1 (MQ=255)
                                                                                                           aCGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATg                                       >  1:1242950/1‑151 (MQ=255)
                                                                                                              ttATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGgtg                                    >  2:1333166/1‑151 (MQ=255)
                                                                                                                    tAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCt                                <  1:1182847/149‑1 (MQ=255)
                                                                                                                     aaTCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTgg                              <  1:971577/150‑1 (MQ=255)
                                                                                                                     aaTCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTgg                              >  2:674583/1‑150 (MQ=255)
                                                                                                                      aTCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGAt                            <  2:1158431/151‑1 (MQ=255)
                                                                                                                           aCTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGgcg                       >  2:1065482/1‑151 (MQ=255)
                                                                                                                            cTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGgcg                       >  1:1035033/1‑150 (MQ=255)
                                                                                                                              aTGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGgcgc                      <  1:1241366/149‑1 (MQ=255)
                                                                                                                               tGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGGGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGa                    >  1:2657/1‑150 (MQ=255)
                                                                                                                                 ccGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGa                    >  2:1242563/1‑148 (MQ=255)
                                                                                                                                  cGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGAt                <  2:522568/151‑1 (MQ=255)
                                                                                                                                    gTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGc              <  2:657285/151‑1 (MQ=255)
                                                                                                                                    gTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGc              <  1:722274/151‑1 (MQ=255)
                                                                                                                                      aGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGt            >  2:1284649/1‑151 (MQ=255)
                                                                                                                                        gCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTa           >  2:514/1‑150 (MQ=255)
                                                                                                                                           tGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAAc         <  1:1237742/149‑1 (MQ=255)
                                                                                                                                            gCCGGGTATGGCAGCTGCCTGTTTTATTCCCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCg       >  1:632594/1‑150 (MQ=255)
                                                                                                                                             ccGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTa     <  1:468557/151‑1 (MQ=255)
                                                                                                                                              cGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTa     >  2:190129/1‑150 (MQ=255)
                                                                                                                                                ggTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGc  <  2:464729/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTGGCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGC  >  NC_000913/4161819‑4162113

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