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.0% 177.3 / NA 50L156* (TTG→TAG) yigADUF484 domain‑containing protein YigA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (19/30);  minor base C (1/0);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTTGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTG  >  NC_000913/3995799‑3996279
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      
gcgcTGTTGATGGAACAGGCTATCGCCAACGACGGCCTGTTTTATCGCCTACTCTACCTGCCGCGCAGTATCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCGGTCAGTCTTTCGAACCGCTGAGAATTCAGCGttt                                                                                                                                                                                                                                        >  2:175753/1‑251 (MQ=255)
             aaCAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAg                                                                                                                                                                                                                           <  1:51754/251‑1 (MQ=255)
                 ggCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACt                                                                                                                                                                                                                       <  2:276745/251‑1 (MQ=255)
                               aaGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTCGGGGCAGGAACAGCACTATCTTGGGCCGCtt                                                                                                                                                                                                         >  2:240272/1‑251 (MQ=255)
                               aaGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCtt                                                                                                                                                                                                         >  1:152608/1‑251 (MQ=255)
                               aaGGCCTGTGTTATCGCCTACTCTTCCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCAGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCtt                                                                                                                                                                                                         <  2:255418/251‑1 (MQ=255)
                                 ggCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTaa                                                                                                                                                                                                       <  2:128714/251‑1 (MQ=255)
                                  gCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAAc                                                                                                                                                                                                      <  1:306318/251‑1 (MQ=255)
                                      gTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCCCCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTTGGCGCGCGATCTCGGCCTGGCAGGTGCTAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCGTTCGAACCGCTGCTTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGa                                                                                                                                                                                                   <  1:210284/250‑1 (MQ=255)
                                      gTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGa                                                                                                                                                                                                   <  2:55774/250‑1 (MQ=255)
                                           aTCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCaga                                                                                                                                                                                              <  1:191070/250‑1 (MQ=255)
                                                cTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAgctgct                                                                                                                                                                                        >  1:278095/1‑251 (MQ=255)
                                                  aCTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGctgctg                                                                                                                                                                                       <  1:48603/250‑1 (MQ=255)
                                                   ctctACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGCGCGCGCGCTCTCGGCCTGGCAGGTGCGCGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGCGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCtggg                                                                                                                                                                                     >  2:180011/1‑250 (MQ=255)
                                                       cccTGCAGCGCAGTCTCCCCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCtggtgg                                                                                                                                                                                   <  2:278095/248‑1 (MQ=255)
                                                           gCAGCGCAGTCTCCCCGCCGCCAGCAGTCGCGACGATATGCTGATGCGCTTTCTCCGCTGGGCGCGCGATCTCGGCCTGGCAGGGGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTGCGAGCCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCGTTTGCCGCTTAACGCACCAGATCTGCCGGTGGTGCTa                                                                                                                                                                              <  2:158073/250‑1 (MQ=255)
                                                           gCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTAc                                                                                                                                                                             >  1:77923/1‑251 (MQ=255)
                                                              gcgcCGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCgg                                                                                                                                                                          <  2:107631/251‑1 (MQ=255)
                                                                gcAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGaa                                                                                                                                                                        <  1:359255/251‑1 (MQ=255)
                                                                   gTCTCCCCGCCGCCAGCAGTCTCGACGTTATGCTGATGCGCTTTCACCGCTGGGCGCTCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTGAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCg                                                                                                                                                                     <  2:265386/251‑1 (MQ=255)
                                                                      tcACCGCCGCCAGCAGTCTCAACGATATGCTGCTGCGCTTTCACCGCTGGGCCCGCGATCTCGGCCTGGCAGGTGCGAGTCTGAGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTCAGCCGTCAGTATTTCGAACCGCTGCGTCTTCAGCGTTAGGGGCAGGAACCGCACACTCTTGGGCCGCTTCACGGACCAGAGCTGCTGGTGGTGGTACCGGACGCGaa                                                                                                                                                                   >  2:218987/1‑250 (MQ=255)
                                                                             cgccAGCAGTCTCGACGATCTGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTgg                                                                                                                                                           >  2:234602/1‑251 (MQ=255)
                                                                                     gTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATAAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTg                                                                                                                                                   <  2:257968/251‑1 (MQ=255)
                                                                                      tctcGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTgg                                                                                                                                                  >  1:93403/1‑251 (MQ=255)
                                                                                              atatGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCg                                                                                                                                          <  2:197216/251‑1 (MQ=255)
                                                                                                         cgcTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGCAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGa                                                                                                                                 >  1:202261/1‑249 (MQ=255)
                                                                                                             ttCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCg                                                                                                                             <  1:240272/249‑1 (MQ=255)
                                                                                                              tCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCGTTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACGATCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGAt                                                                                                                           <  2:347206/250‑1 (MQ=255)
                                                                                                               cACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATg                                                                                                                          <  2:264363/250‑1 (MQ=255)
                                                                                                               cACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTAGCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGc                                                                                                                         <  2:348447/251‑1 (MQ=255)
                                                                                                                            gcgATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTc                                                                                                            <  1:364006/251‑1 (MQ=255)
                                                                                                                                tctcGGCCTGGCAGGGGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGAACGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGGGTCGTGc                                                                                                        >  2:271700/1‑251 (MQ=255)
                                                                                                                                   cGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGt                                                                                                     <  1:150799/251‑1 (MQ=255)
                                                                                                                                    ggCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGtt                                                                                                    <  1:28587/251‑1 (MQ=255)
                                                                                                                                            aGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCGGATTTGGGTGTCGTGCTGTTTACCAGTc                                                                                            >  2:19419/1‑251 (MQ=255)
                                                                                                                                              gTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGCTTTGGGGGTCGTGCTGTTTACCAGTcgc                                                                                          >  2:70594/1‑251 (MQ=255)
                                                                                                                                                  gAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATg                                                                                      <  1:145340/251‑1 (MQ=255)
                                                                                                                                                         cgcCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCa                                                                               <  1:250971/251‑1 (MQ=255)
                                                                                                                                                                 tCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGGAGCGATGCTGATTTGGGTGGCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATcaac                                                                       >  2:364558/1‑251 (MQ=255)
                                                                                                                                                                            ggCGCGTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGGTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGa                                                            <  2:166739/251‑1 (MQ=255)
                                                                                                                                                                                      tGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTCCCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTAc                                                  >  1:133231/1‑251 (MQ=255)
                                                                                                                                                                                           cGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGACCAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCAt                                             >  1:212740/1‑251 (MQ=255)
                                                                                                                                                                                            gTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATg                                            >  2:363090/1‑251 (MQ=255)
                                                                                                                                                                                                        tCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGGGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGCTTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTg                                 >  2:202945/1‑250 (MQ=255)
                                                                                                                                                                                                         cATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGAt                               <  1:187371/251‑1 (MQ=255)
                                                                                                                                                                                                         cATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGAt                               <  1:104330/251‑1 (MQ=255)
                                                                                                                                                                                                            cTGGCATTAAGCCGTCAGTCTTTCGACCCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATcaa                                                                        >  1:289719/1‑207 (MQ=255)
                                                                                                                                                                                                            cTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATcaa                                                                        <  2:289719/207‑1 (MQ=255)
                                                                                                                                                                                                              ggCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATg                              <  1:285311/247‑1 (MQ=255)
                                                                                                                                                                                                                       ccGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTc                 >  1:115807/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ccGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTg  >  1:372143/1‑250 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      
GCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTTGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTG  >  NC_000913/3995799‑3996279

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