Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,996,049 T→A L156* (TTG→TAG)  yigA → DUF484 domain‑containing protein YigA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,996,0490TA98.4% 190.2 / NA 63L156* (TTG→TAG) yigADUF484 domain‑containing protein YigA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base A (30/31);  minor base C (1/0);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTTGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGTATGACCGAT  >  NC_000913/3995803‑3996295
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
tggtgatgaAACAGGGTATCGCCAACGAAGGCCTGTTTTATCGCCTACGCTCCCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAggg                                                                                                                                                                                                                                                     <  2:163551/247‑1 (MQ=255)
 gTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGc                                                                                                                                                                                                                                                   <  1:201164/251‑1 (MQ=255)
 gTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTACCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGGCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAgggg                                                                                                                                                                                                                                                    <  1:302644/250‑1 (MQ=255)
   tGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGc                                                                                                                                                                                                                                                   >  1:398989/1‑249 (MQ=255)
   ggATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCa                                                                                                                                                                                                                                                  <  1:213324/249‑1 (MQ=255)
       ggAACAGGTTCTCGCCAACGAAGGCCTTTTTTATCGCCTCCTCTACCTGCACCGCATTCTCCCCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGaa                                                                                                                                                                                                                                              <  2:392204/250‑1 (MQ=255)
       ggAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGaa                                                                                                                                                                                                                                              <  1:68438/250‑1 (MQ=255)
          aCAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGc                                                                                                                                                                                                                                          >  1:150849/1‑251 (MQ=255)
           cAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGc                                                                                                                                                                                                                                          >  1:100064/1‑250 (MQ=255)
                tATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATc                                                                                                                                                                                                                                    >  1:448177/1‑251 (MQ=255)
                          gAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCAACGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCAACCTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCGGCGTATTCAGCGTTTGGGGCAGGaaacgcaccagctttggccgat                                                                                                                                                                                                                          >  2:303826/1‑231 (MQ=255)
                           aaGGCCTGTTTTATCGACTACTCCACCTGCAGCGGAGTCCCACCGCCGCCAGCAGTCCCGAAGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCtt                                                                                                                                                                                                                         <  2:14315/251‑1 (MQ=255)
                              gCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTCTCTTGGGCCGCTTAAc                                                                                                                                                                                                                      >  2:268092/1‑251 (MQ=255)
                              gCCTGTTTTATCGCCTACTCTACATGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGTTATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCAGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGTGGCCGGAACAACACTAACGTGGGcagattaa                                                                                                                                                                                                                       >  1:457939/1‑243 (MQ=255)
                                    tttATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGAcc                                                                                                                                                                                                                 <  1:116837/250‑1 (MQ=255)
                                          gCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAgctg                                                                                                                                                                                                          >  1:179488/1‑251 (MQ=255)
                                           cctcctctACCTGCAGCGCAGTCCCACAGCCGCCAGCAGTCTCGACGAGATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAgc                                                                                                                                                                                                            <  2:399852/244‑1 (MQ=255)
                                              aCTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCtgg                                                                                                                                                                                                      >  1:232016/1‑251 (MQ=255)
                                                 cacccTGCGGCGCAGTCTCCCCGCCGCCAGCAGGCTCGACGATATGCTGATGCGCTTTCACCGCCGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCtggtgg                                                                                                                                                                                                   <  2:108743/248‑1 (MQ=255)
                                                             cAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGCCCAGAGCTGCTGGTGGGGCTACCGGa                                                                                                                                                                                         >  2:139697/1‑249 (MQ=255)
                                                                        ccgccAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTg                                                                                                                                                                            >  1:202096/1‑251 (MQ=255)
                                                                                      gacgaTATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCCGCGCTGCTGGGGGGGCTACCGTAAGCGAAAGCGGGGGGAGCGGTTGCggg                                                                                                                                                              >  2:398397/1‑249 (MQ=255)
                                                                                      gacgaTATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGGGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAAAAGCACTATCTTGGGCCGCTTAACGGCCAAGAGCTGCTGGTGGGGCTACCGGAAGCGagagggggggggacggggggggg                                                                                                                                                              >  2:190886/1‑229 (MQ=255)
                                                                                         gaTATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGTTCGGTGGCggtg                                                                                                                                                             >  1:267285/1‑246 (MQ=255)
                                                                                          atatGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCAATATCTTTGGCCGCTTAACGGACCAGCGCTGCTGGTGGTGCTACCGGCAGCGAAATGGGGGGGACCGGTGGGGATGTCg                                                                                                                                                          >  2:126397/1‑251 (MQ=255)
                                                                                            atGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGAt                                                                                                                                                        <  1:438006/251‑1 (MQ=255)
                                                                                                   tGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGAGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGGGGGATCGGTGGCGATGTCGATGCTggg                                                                                                                                                  >  1:147130/1‑250 (MQ=255)
                                                                                                    gcgcTTTCACCGCTGGGCGCGCGAGCTCGGCCTGGCAGGTGCTAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCATTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGaa                                                                                                                                                <  2:135327/251‑1 (MQ=255)
                                                                                                          tCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTTTGCGCCTGTGTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCCTCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGCACAGCACTTTCTTGGGCCGCTTAACGGACAAGGGCTGCTGGTGGGGCTACCGGAAGCGAAAGCGGGGGGATTGGTGGCGAGGTCGATGCTggggagtgagg                                                                                                                                          >  2:54398/1‑243 (MQ=255)
                                                                                                          tCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGAt                                                                                                                                           <  1:139697/250‑1 (MQ=255)
                                                                                                              cGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAGGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTg                                                                                                                                       >  1:396883/1‑250 (MQ=255)
                                                                                                                  gggCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTg                                                                                                                                  <  1:188309/251‑1 (MQ=255)
                                                                                                                         cgATCGCGGCCTGGTAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACTCTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCg                                                                                                                           <  1:303826/251‑1 (MQ=255)
                                                                                                                                ggCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGtt                                                                                                                    >  1:167095/1‑251 (MQ=255)
                                                                                                                                        aGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTc                                                                                                            >  1:412842/1‑251 (MQ=255)
                                                                                                                                         ggTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTcg                                                                                                           <  1:259481/251‑1 (MQ=255)
                                                                                                                                          gTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGGGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTAAGTCTTTCGAACCGCTGCGTATTCAGCGATAGGGGCAGGAACAGCGCTATCTTGGGACGCTTAACGGAGAAGAGCTGCTGGGGGAGGTACCGGAGGCGAACGCGGTGGGCTCGGTGGCGATGATGTTGCTGGTAAGCGATGGTGGTTTgggggcggcgttgggtgacagtggg                                                                                                          >  2:32224/1‑229 (MQ=255)
                                                                                                                                             cGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGGGGGATCGGTGGCGATGTCGGTGCTGGGGAGCGATGCTGATTTGGGGGTCGTGCTGTTTACCAGGcgcg                                                                                                         >  2:228874/1‑249 (MQ=255)
                                                                                                                                                    gcgcCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTc                                                                                                <  1:138586/251‑1 (MQ=255)
                                                                                                                                                        ccgTTTCCGGATCGCTGGCCCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGCATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTa                                                                                            <  1:162220/249‑1 (MQ=255)
                                                                                                                                                                     gCTGGCTCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTGTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCaa                                                                               <  2:449069/251‑1 (MQ=255)
                                                                                                                                                                        ggCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAgg                                                                             <  2:179332/250‑1 (MQ=255)
                                                                                                                                                                                ggggcgcCGGCGAACCACACTCATCCGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGAAGCTGGGAAGCGATGCTGAATTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGtt                                                                    <  2:232016/248‑1 (MQ=255)
                                                                                                                                                                                             aCCACACTCATCTGGCATTAAGCCGTCGGTCTTTCGAACCGCTGCTGATTCAGCGGTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGGAGTTACTTCATGAAAtt                                                       <  2:54103/251‑1 (MQ=255)
                                                                                                                                                                                                  acTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCt                                                  >  1:435522/1‑251 (MQ=255)
                                                                                                                                                                                                  acTCATCTGGCATTAAGCCGACAGGCTTTCGAACCGCTGCGTCTTCAGCGTTAGGGGCAGGCACAGCACTATCTTGGTCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGACAGCGGTGGGAGCGGTGGCGACGTCGGTGCTGGGCAGCGATTCTGCTTTGGGTGTCGTGCTGTTTACCCGTCGCGATCCCAGTCACCATCAACAAGGGCAAGGAACGCCGTTACTTCATTAAATTgcact                                                  >  2:284550/1‑248 (MQ=255)
                                                                                                                                                                                                     cATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGAt                                               <  2:412842/251‑1 (MQ=255)
                                                                                                                                                                                                          ggCATTAAGCCGTCAGGTGTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTGGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGc                                          <  1:316589/251‑1 (MQ=255)
                                                                                                                                                                                                            cATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGCACAGCACTATCTTGGGCCGCTTAACGGCCCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGGTCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGGCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTAATTCATGAAATTGCGCTGATGTTgacg                                        >  2:285594/1‑248 (MQ=255)
                                                                                                                                                                                                              ttAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGCTGCTGGGAAGCGGTGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCCACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCggg                                      >  2:14224/1‑250 (MQ=255)
                                                                                                                                                                                                                aaGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGc                                    >  1:428196/1‑251 (MQ=255)
                                                                                                                                                                                                                    cGTCAGTCTGTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCt                                <  1:59139/251‑1 (MQ=255)
                                                                                                                                                                                                                     gTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCCGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGTTTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGAAAGGAACGCAGTTACTTCATGAAATTGCGCTGAAGTTGCCGGAGCTTcgg                               >  2:117762/1‑249 (MQ=255)
                                                                                                                                                                                                                             ttCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGtt                        <  2:448177/250‑1 (MQ=255)
                                                                                                                                                                                                                             ttCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTGGCCGGAGCTTCTGGAGCg                          <  1:126397/248‑1 (MQ=255)
                                                                                                                                                                                                                                 aaCCGCAGCGTATTCAGCGTTCGTGGCAGGTACCGCACTATCCTGGGCCGCTTAACGGACCAGAGCAGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGCTGTCGTTGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCCGTCGCGATGCCAGTCAATATCAACAAGGGCAAGGAACGCAGTTAATTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTTGAtt                   >  2:60195/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cGCTGCGTATTCAGCGTTAAGGGAAGGCACAACACTATCTTGGGCCGCTTAACGGACCAGGGATGCTGGTGGTGCTACCGGAAGCGAAAGCGGAGGGATCGGTGGCGCTGTCGATGCTGGGAAGCGCTGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACGATCAACAAGGGCAAGGAACGCAGTTACTTCATAGAATTGCGCTGATGTTGACGGatcgtttggataagtggcttga                 >  2:393299/1‑229 (MQ=255)
                                                                                                                                                                                                                                     gCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGCGCTGCTGGTGGTGATACCGGAAGCGAAAGCGGTGGGATCGGTGGCGAAGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCCTGAAATTGCGCTGATGTTGCCGGAGATTCTGGAGCGTTGGATTGAAc               >  1:237687/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCGTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGAGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGatcg              <  2:93570/250‑4 (MQ=255)
                                                                                                                                                                                                                                        gCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAAcgc             <  1:292161/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gTATTCAGGGTGAGGGGCAGGAACAGCACTATCGTGTGCCGCTTAACGGTCCTGAGCTGCTGGTGGTGCTACCTGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCGGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTTATTAAATTGCGCTGATGTGGCCGGAGCTTCTGGAGCGTTGGATTGAATGCGTa          <  2:431270/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGTTAGGGGCAGGCACAGCACTATCTTGGGCCGCTTAACGGACCAGAGATGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGGGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGGTGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGTGGATGACCGa   >  2:409531/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  cGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGTATGACCGAt  <  1:14224/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
TGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTTGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGTATGACCGAT  >  NC_000913/3995803‑3996295

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