Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,449,361 G→C 7.1% I107M (ATC→ATG)  yhaV ← toxin antitoxin system SohB(PrlF)‑YhaV toxin endoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,449,3610GC 7.1% 209.6 / 6.8 70I107M (ATC→ATG) yhaVtoxin antitoxin system SohB(PrlF)‑YhaV toxin endoribonuclease
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base C (3/2);  total (39/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 = 7.66e-01

AATCTGCGGCAGCCTCTTGCAGCAACTGGTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAAC  >  NC_002947/1449216‑1449508
                                                                                                                                                 |                                                                                                                                                    
aaTCTGCGGCAGCCGCTTGCAGCAACTGGTCCCAATCGTGCGGGGGGTTGCCGGTAAGCAACATCTTTTGAAAGAACTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGTTCGACTCGTCATTGATCCAGGCCAGCACGATGATc                                                                                                                                                   <  6:238943/149‑1 (MQ=255)
 aTCTGCGGCAGCCTCTTGCAGCAACTGGTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATcc                                                                                                                                                  >  8:100714/1‑149 (MQ=255)
  tCTGCGGCAGCATCTTGCAGCAACTGGTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGATAGACTCGTAATTGATCCAGGCCAGCCCGATGATCct                                                                                                                                                 >  6:360369/1‑149 (MQ=255)
      cGGCAGCCTCTTGCAGCAACTGGTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCctg                                                                                                                                                <  7:100714/146‑1 (MQ=255)
            ccTCTTGCAGCAACTGGTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGgtg                                                                                                                                       >  1:212766/1‑149 (MQ=255)
                  gcagcaACTGGTCCCACTCGTCCGGGTGGTTGCCGCTATGCAACTTCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTgg                                                                                                                                    <  2:395236/146‑1 (MQ=255)
                   cagcaACTGGTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATTGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGc                                                                                                                                <  4:125115/149‑1 (MQ=255)
                   cagcaACTGGTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGc                                                                                                                                <  1:214969/149‑1 (MQ=255)
                       aaCTGGTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGgaa                                                                                                                            <  4:133814/149‑1 (MQ=255)
                       aaCTGGTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGgaa                                                                                                                            <  1:8315/149‑1 (MQ=255)
                            gTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGaac                                                                                                                        >  6:32081/1‑148 (MQ=255)
                              cccAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGaacaac                                                                                                                     >  6:41577/1‑149 (MQ=255)
                                  aTCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCgg                                                                                                                  >  7:186314/1‑148 (MQ=255)
                                      tCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTActg                                                                                                             <  5:209501/149‑1 (MQ=255)
                                      tCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTActg                                                                                                             >  7:133314/1‑149 (MQ=255)
                                      tCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTActg                                                                                                             <  8:186314/149‑1 (MQ=255)
                                                  acGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCTCTACTTTCATATGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGgaatatct                                                                                                  <  5:360369/147‑6 (MQ=255)
                                                        tGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGctct                                                                                           <  4:135996/149‑1 (MQ=255)
                                                          cAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTg                                                                                          >  5:177184/1‑148 (MQ=255)
                                                           aaCATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGaa                                                                                        >  2:391463/1‑149 (MQ=255)
                                                                 ttttGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACCGCTGGAAGAACTTTGCTCTGAACCAGTg                                                                                  >  8:9874/1‑149 (MQ=255)
                                                                  tttGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGc                                                                                 >  6:201790/1‑149 (MQ=255)
                                                                      aaaGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCAGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTtgt                                                                             >  6:161587/1‑149 (MQ=255)
                                                                            cTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACCGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCAc                                                                       <  7:9874/149‑1 (MQ=255)
                                                                                tAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGcc                                                                   <  5:41577/149‑1 (MQ=255)
                                                                                tAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGc                                                                    >  7:26894/1‑148 (MQ=255)
                                                                                 aaGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGcc                                                                   >  1:221745/1‑148 (MQ=255)
                                                                                 aaGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGGGAGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGcc                                                                   >  7:249941/1‑148 (MQ=255)
                                                                                  aGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGcc                                                                   >  8:354923/1‑147 (MQ=255)
                                                                                     aTCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAg‑tg                                                             >  1:4110/1‑149 (MQ=255)
                                                                                        gtcgCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTg                                                          <  8:133314/149‑1 (MQ=255)
                                                                                         tcgCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGc                                                         <  5:32081/149‑1 (MQ=255)
                                                                                            cTACTTTCATAGGCCCGCGTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGcccc                                                      <  4:378030/149‑1 (MQ=255)
                                                                                                tttCATAGGCCCGCTTGCTCGCCTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGc                                                  >  6:43539/1‑149 (MQ=255)
                                                                                                tttCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACACCCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGc                                                  >  3:352362/1‑149 (MQ=255)
                                                                                                tttCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGCAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGc                                                  >  2:181493/1‑149 (MQ=255)
                                                                                                  tCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGcc                                                 <  1:181493/148‑1 (MQ=255)
                                                                                                    aTAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCgg                                               >  1:338379/1‑148 (MQ=255)
                                                                                                          ccGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCgg                                        >  1:18631/1‑149 (MQ=255)
                                                                                                                gCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTgg                                  >  3:282390/1‑149 (MQ=255)
                                                                                                                 cTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCTGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGATTGTCGTCACCGCCAAG‑TGTTTCCCATTGACGGTATTCGGGTTTgg                                  >  4:105515/1‑148 (MQ=255)
                                                                                                                    gACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCgg                              >  1:197246/1‑149 (MQ=255)
                                                                                                                    gACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCgg                              <  6:177184/149‑1 (MQ=255)
                                                                                                                    gACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGCATTCGGGTGTGGTc                                <  2:338379/147‑1 (MQ=255)
                                                                                                                         gTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTc                           <  7:354923/147‑1 (MQ=255)
                                                                                                                         gTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCtt                         >  2:107084/1‑149 (MQ=255)
                                                                                                                          tCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCtt                         >  3:26383/1‑148 (MQ=255)
                                                                                                                            aTTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCgg                                                                                                                              <  3:82062/46‑1 (MQ=255)
                                                                                                                            aTTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCgg                                                                                                                              >  4:82062/1‑46 (MQ=255)
                                                                                                                            aTTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGa                       <  4:26383/148‑1 (MQ=255)
                                                                                                                            aTTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAg                      <  8:26894/149‑1 (MQ=255)
                                                                                                                            aTTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAg                      <  2:197246/149‑1 (MQ=255)
                                                                                                                            aTTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAg                      <  2:4110/149‑1 (MQ=255)
                                                                                                                            aTTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAg                      >  7:260073/1‑149 (MQ=255)
                                                                                                                            aTTGATCCAGGCCAGCACGATGATACTGCTGGGGGCGTGGTATCGGAAGAACAACCGGTACTGCTGGACGAACTTTCCGCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAg                      <  8:266798/149‑1 (MQ=255)
                                                                                                                             ttGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAgg                     >  7:19369/1‑149 (MQ=255)
                                                                                                                              tGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAgg                     >  4:11743/1‑148 (MQ=255)
                                                                                                                                aTCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGAt                  >  7:301556/1‑149 (MQ=255)
                                                                                                                                 tCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGAtt                 <  6:378893/149‑1 (MQ=255)
                                                                                                                                 tCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACAAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGAt                  >  6:18429/1‑148 (MQ=255)
                                                                                                                                 tCCAGGCCAGAACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTATGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCACTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGAt                  >  8:282166/1‑148 (MQ=255)
                                                                                                                                 acacccccAGCAGGATCATCGTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCgg                                        >  4:351579/7‑126 (MQ=255)
                                                                                                                                 acacccccAGCAGGATCATCGTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCgg                                        <  3:351579/120‑1 (MQ=255)
                                                                                                                                  ccAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTa                >  4:294246/1‑149 (MQ=255)
                                                                                                                                  ccAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTa                >  5:371506/1‑149 (MQ=255)
                                                                                                                                     ccccAGCAGGATCATCGTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACAt             >  7:332943/3‑149 (MQ=255)
                                                                                                                                      cccAGCAGGATCATCGTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCgg                              >  2:37094/2‑131 (MQ=255)
                                                                                                                                      cccAGCAGGATCATCGTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCgg                              <  1:37094/130‑1 (MQ=255)
                                                                                                                                          gCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAAc        <  1:391463/149‑1 (MQ=255)
                                                                                                                                             cGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGCTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGc      <  1:107084/149‑1 (MQ=255)
                                                                                                                                               atgatCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCaa   <  8:19369/149‑1 (MQ=255)
                                                                                                                                                tgatCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAAc  >  2:70228/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AATCTGCGGCAGCCTCTTGCAGCAACTGGTCCCAATCGTCCGGTGGGTTGCCGCTATGCAACATCTTTTGAAAGACCTTGTAAGCATCGTCGCTACTTTCATAGGCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAG‑TGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAAC  >  NC_002947/1449216‑1449508

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