Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,147,196 C→T P135S (CCG→TCG)  yceD → DUF177 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,1960CT96.8% 1441.4 / 22.3 434P135S (CCG→TCG) yceDDUF177 family protein
Reads supporting (aligned to +/- strand):  ref base C (8/6);  new base T (208/211);  total (216/218)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 434 total aligned reads displayed.

CACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTCCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAAC  >  NC_000913/1147048‑1147345
                                                                                                                                                    |                                                                                                                                                     
cacaACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTCCg                                                                                                                                                     <  8:1536580/151‑1 (MQ=255)
  caACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTCCGGt                                                                                                                                                   <  7:1102805/151‑1 (MQ=255)
  caACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTAGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGt                                                                                                                                                   <  8:828723/151‑1 (MQ=255)
    aCGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGc                                                                                                                                                 >  6:1377447/1‑151 (MQ=255)
     cGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCa                                                                                                                                                <  8:462482/151‑1 (MQ=255)
      gTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTNCCGGAAGCGTATGAACCGNNNGNNNTNNNNGNNTTCNGNNAAATCGNTCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCNGTGCAt                                                                                                                                               >  2:319830/1‑151 (MQ=255)
             tttAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCt                                                                                                                                         <  6:197305/150‑1 (MQ=255)
             tttAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCt                                                                                                                                         <  3:1297859/150‑1 (MQ=255)
              ttAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTg                                                                                                                                        >  8:666112/1‑150 (MQ=255)
                 gTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAac                                                                                                                                     >  7:907234/1‑150 (MQ=255)
                    cTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACtg                                                                                                                                 >  6:394578/1‑151 (MQ=255)
                         cGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGa                                                                                                                              >  8:129991/1‑149 (MQ=255)
                          gTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAgt                                                                                                                           >  6:876720/1‑151 (MQ=255)
                            tCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAgtgt                                                                                                                         >  6:654630/1‑151 (MQ=255)
                                aCGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCg                                                                                                                      >  7:1352062/1‑150 (MQ=255)
                                aCGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCg                                                                                                                      <  1:1520966/150‑1 (MQ=255)
                                aCGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGa                                                                                                                     <  1:1441980/151‑1 (MQ=255)
                                    aCAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCg                                                                                                                 <  2:772120/151‑1 (MQ=255)
                                     cAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTCCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCgg                                                                                                                >  3:1638724/1‑151 (MQ=255)
                                      aGGCTGAAGCACTGACGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTCCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGa                                                                                                               <  5:1657801/151‑1 (MQ=255)
                                          tGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATg                                                                                                           >  6:1489135/1‑151 (MQ=255)
                                            aaGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGt                                                                                                         <  5:1456780/151‑1 (MQ=255)
                                             aGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTCAAGTGTCCGAAGCGGACATGGTc                                                                                                        >  4:1570322/1‑151 (MQ=255)
                                              gCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCt                                                                                                       >  2:917593/1‑151 (MQ=255)
                                              gCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCt                                                                                                       >  5:1618494/1‑151 (MQ=255)
                                               cACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCtt                                                                                                      <  3:462855/151‑1 (MQ=255)
                                                  tGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTgg                                                                                                   >  2:1410651/1‑151 (MQ=255)
                                                   gCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTgg                                                                                                   <  7:1287589/150‑1 (MQ=255)
                                                    ccGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTg                                                                                                 >  3:394643/1‑151 (MQ=255)
                                                     cGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTg                                                                                                 >  6:132017/1‑150 (MQ=255)
                                                     cGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGa                                                                                                <  3:1197913/151‑1 (MQ=255)
                                                         aGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTg                                                                                            >  4:1434000/1‑151 (MQ=255)
                                                           cGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGcc                                                                                          >  5:707413/1‑151 (MQ=255)
                                                             tATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTg                                                                                        >  2:1643254/1‑151 (MQ=255)
                                                              aTGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTga                                                                                       <  2:593039/151‑1 (MQ=255)
                                                              aTGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTg                                                                                        <  4:8569/150‑1 (MQ=255)
                                                              aTGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCCCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTga                                                                                       <  1:386875/151‑1 (MQ=255)
                                                              aTGAACCGATTGAGGTTAACAAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTg                                                                                        <  8:38128/150‑1 (MQ=255)
                                                                gAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTgaag                                                                                     <  2:1281669/151‑1 (MQ=255)
                                                                   ccGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTgaaga                                                                                    >  6:973799/1‑149 (MQ=255)
                                                                   ccGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGaagaag                                                                                  <  3:762345/151‑1 (MQ=255)
                                                                   ccGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGaagaag                                                                                  <  7:789508/151‑1 (MQ=255)
                                                                    cGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAgc                                                                                 >  2:1220653/1‑151 (MQ=255)
                                                                    cGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGAAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAgc                                                                                 >  7:1424282/1‑151 (MQ=255)
                                                                      aTTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAgcgc                                                                               <  2:151484/151‑1 (MQ=255)
                                                                       ttGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGACTTGCAGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGAGGACATGGTCTTTGGTGAACTGCCTGaagaagaga                                                                               >  7:1498167/1‑147 (MQ=255)
                                                                         gAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCaa                                                                             >  4:1271241/1‑150 (MQ=255)
                                                                         gAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCaa                                                                             >  5:1543910/1‑150 (MQ=255)
                                                                          aGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCaaaa                                                                           <  8:246559/151‑1 (MQ=255)
                                                                          aGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCaaaa                                                                           <  5:1226328/151‑1 (MQ=255)
                                                                            gTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGc                                                                         <  4:200650/151‑1 (MQ=255)
                                                                              tAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGcc                                                                        <  8:1096748/150‑1 (MQ=255)
                                                                               aaCGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCTGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCa                                                                       <  4:633473/150‑1 (MQ=255)
                                                                               aaCGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCa                                                                       >  3:224256/1‑150 (MQ=255)
                                                                                aCGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCaa                                                                      >  2:999414/1‑150 (MQ=255)
                                                                                 cGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCaacc                                                                    >  1:1269370/1‑149 (MQ=255)
                                                                                   aaTTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAAccc                                                                  <  7:646737/151‑1 (MQ=255)
                                                                                    aTTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCa                                                                 <  8:1375178/151‑1 (MQ=255)
                                                                                       cGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCAttt                                                              <  8:1600734/151‑1 (MQ=255)
                                                                                        ggTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTg                                                             <  5:570754/151‑1 (MQ=255)
                                                                                        ggTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTg                                                             <  4:742371/151‑1 (MQ=255)
                                                                                           gAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCg                                                          >  3:1418056/1‑151 (MQ=255)
                                                                                            aaaTCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGt                                                         <  4:498515/151‑1 (MQ=255)
                                                                                             aaTCGATCTGCTTGCAATCGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTCAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGt                                                         >  4:1185276/1‑150 (MQ=255)
                                                                                              aTCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTCCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTAt                                                       <  1:1634732/151‑1 (MQ=255)
                                                                                                 gATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTAtt                                                      >  4:459416/1‑149 (MQ=255)
                                                                                                      gCTTGCAATGGTTGAAGATGTAATCATCGTCGCCTTGGCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAg                                                <  3:510778/150‑1 (MQ=255)
                                                                                                      gCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAg                                                <  5:362046/150‑1 (MQ=255)
                                                                                                      gCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGc                                               <  2:213822/151‑1 (MQ=255)
                                                                                                         tGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCtt                                             <  6:288595/150‑1 (MQ=255)
                                                                                                         tGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCtt                                             >  7:179793/1‑150 (MQ=255)
                                                                                                           cAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTaaa                                          <  7:31149/151‑1 (MQ=255)
                                                                                                            aaTGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCCTTTGCCGTATTAGCCAGCTTaa                                           >  3:1495607/1‑149 (MQ=255)
                                                                                                              tGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCg                                       >  5:1133066/1‑151 (MQ=255)
                                                                                                                gTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCgta                                     <  1:485854/151‑1 (MQ=255)
                                                                                                                gTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCgt                                      >  5:922132/1‑150 (MQ=255)
                                                                                                                gTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTCCCGAATTAGcaccattaaagctt                                      >  7:1617961/1‑137 (MQ=255)
                                                                                                                   gAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCCCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCgtaagt                                  <  3:979162/151‑1 (MQ=255)
                                                                                                                    aaGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCgtaagt                                  >  7:508713/1‑150 (MQ=255)
                                                                                                                     aGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCgtaagta                                 >  3:300473/1‑150 (MQ=255)
                                                                                                                          aaaTCATGCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGt                           <  8:1013756/151‑1 (MQ=255)
                                                                                                                           aaTCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAATCCGATTTGCCGTATTAGCCAGCTTAACGCGTAAGTAATTGGTg                          >  6:507153/1‑151 (MQ=255)
                                                                                                                            atcatcCTCGCCTTGCCGGTAGTTTCGGTGCATGCTTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGc                         <  7:1210541/151‑1 (MQ=255)
                                                                                                                              catcCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTc                       >  2:827956/1‑151 (MQ=255)
                                                                                                                                    cGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGt                  <  5:1364867/150‑1 (MQ=255)
                                                                                                                                    cGCCTTGCCGGTAGTTTCGGGGGATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGTAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCGCCCCGtt                 <  5:1027122/151‑1 (MQ=255)
                                                                                                                                     gCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTg                >  6:721045/1‑151 (MQ=255)
                                                                                                                                     gCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTg                <  6:1585033/151‑1 (MQ=255)
                                                                                                                                     gCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGATg                >  1:265541/1‑151 (MQ=255)
                                                                                                                                      ccTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTgg               >  8:894489/1‑151 (MQ=255)
                                                                                                                                      ccTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTgg               <  2:522915/151‑1 (MQ=255)
                                                                                                                                       cTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTTTACGAAGCGGACATGGTCAATGGTGAACTGCCAGAAGAAGCGCAAAAGCCAAAACCATTTGCAGTATTAGCCAGCTTAAAGCCTAAGTAATTGGTGCTCCCCGTTgg               >  6:1189584/1‑150 (MQ=255)
                                                                                                                                       cTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTgg               >  7:428703/1‑150 (MQ=255)
                                                                                                                                       cTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTgg               >  5:853350/1‑150 (MQ=255)
                                                                                                                                         tGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGAt             <  1:799447/150‑1 (MQ=255)
                                                                                                                                         tGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAAGTGGTGCTCCCCGTGGGagg            >  1:1184290/1‑149 (MQ=255)
                                                                                                                                            cGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCggg         <  1:709153/151‑1 (MQ=255)
                                                                                                                                              gTAGTTTCGGTGCATGATTCTGCACACTGTGACGTGTCCGAAGCGGACATGGGCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCACCCCGTTGGATCGGGGa       <  6:1116006/151‑1 (MQ=255)
                                                                                                                                               tAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGAt      >  1:406516/1‑151 (MQ=255)
                                                                                                                                                   tCCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAAc  >  1:1118272/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTCCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAAC  >  NC_000913/1147048‑1147345

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