Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,887,352 T→A intergenic (‑133/+226) ygcB ← / ← cysH Cascade complex anti‑viral R‑loop helicase‑annealase Cas3/phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,887,3520TA100.0% 415.6 / NA 102intergenic (‑133/+226)ygcB/cysHCascade complex anti‑viral R‑loop helicase‑annealase Cas3/phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (68/34);  total (68/34)

Only 100 of 102 total aligned reads displayed.

ATTTAAAAGGTTCCATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCA  >  NC_000913/2887204‑2887495
                                                                                                                                                    |                                                                                                                                               
atTTAAAAGGTTCCATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGt                                                                                                                                               >  1:1376611/1‑151 (MQ=255)
 tttAAAAGGTTCCATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGtt                                                                                                                                              >  1:1470334/1‑151 (MQ=255)
  ttAAAAGGTTCCATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGttt                                                                                                                                             <  1:1153938/151‑1 (MQ=255)
        ggTTCCATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGtt                                                                                                                                        >  2:72937/1‑150 (MQ=255)
         gTTCCATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCg                                                                                                                                      >  1:1686317/1‑151 (MQ=255)
         gTTCCATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCg                                                                                                                                      >  2:1577396/1‑151 (MQ=255)
         gTTCCATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCATGCAAAGAAGAGACTAACATTCGAATAGTTTGGGtt                                                                                                                                        <  2:1374575/149‑1 (MQ=255)
           tCCATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGaa                                                                                                                                    >  2:1468254/1‑151 (MQ=255)
             cATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACg                                                                                                                                  >  1:424371/1‑151 (MQ=255)
             cATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACg                                                                                                                                  >  2:313939/1‑151 (MQ=255)
              aTTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGc                                                                                                                                 <  1:1546431/151‑1 (MQ=255)
               ttAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACg                                                                                                                                  >  1:1076163/1‑149 (MQ=255)
                        tCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAgg                                                                                                                       <  1:1061548/151‑1 (MQ=255)
                         cccTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGt                                                                                                                      <  2:1456753/151‑1 (MQ=255)
                             gTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                                                                  >  2:1500459/1‑151 (MQ=255)
                              tttttttAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATAT‑‑ACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg                                                                                                                <  2:132956/150‑1 (MQ=255)
                               ttttttAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg                                                                                                                <  1:1251203/151‑1 (MQ=255)
                               ttttttAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg                                                                                                                <  1:348467/151‑1 (MQ=255)
                                 ttttAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaa                                                                                                              >  1:1121268/1‑151 (MQ=255)
                                             cGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTggg                                                                                                  >  1:271534/1‑151 (MQ=255)
                                             cGATGAACTCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTggg                                                                                                  >  1:748/1‑151 (MQ=255)
                                                 gatATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCtt                                                                                              >  1:1074197/1‑151 (MQ=255)
                                                         cGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                                                                         <  1:29219/148‑1 (MQ=255)
                                                         cGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                                                                         <  1:1714529/148‑1 (MQ=255)
                                                           ggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                                                                     <  2:401470/150‑1 (MQ=255)
                                                                  tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGGAATATCGCCTGGGGCTTTTGTCCATCTGGAAc                                                                               >  1:1176911/1‑149 (MQ=255)
                                                                  tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt                                                                             <  1:585660/151‑1 (MQ=255)
                                                                          ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt                                                                      <  1:496327/150‑1 (MQ=255)
                                                                            gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                                                                     >  1:39991/1‑149 (MQ=255)
                                                                            gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                                                                     >  1:1523839/1‑149 (MQ=255)
                                                                            gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                                                                     >  1:997772/1‑149 (MQ=255)
                                                                            gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                                                                     >  2:559929/1‑149 (MQ=255)
                                                                            gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                                                                     >  1:1739702/1‑149 (MQ=255)
                                                                               cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTaa                                                                >  2:1176163/1‑151 (MQ=255)
                                                                                  aaaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc                                                             <  2:278838/151‑1 (MQ=255)
                                                                                  aaaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc                                                             <  2:422781/151‑1 (MQ=255)
                                                                                       tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGa                                                         >  2:1688376/1‑150 (MQ=255)
                                                                                       tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc                                                        <  2:1339440/151‑1 (MQ=255)
                                                                                       tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc                                                        <  2:106184/151‑1 (MQ=255)
                                                                                       tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc                                                        <  1:902422/151‑1 (MQ=255)
                                                                                       tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc                                                        <  1:840880/151‑1 (MQ=255)
                                                                                          aaTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAg                                                      >  2:1425611/1‑150 (MQ=255)
                                                                                           aTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc                                                     >  1:891625/1‑150 (MQ=255)
                                                                                            ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc                                                     >  2:1639339/1‑149 (MQ=255)
                                                                                                 gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAg                                               >  2:1520145/1‑150 (MQ=255)
                                                                                                 gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAg                                               >  1:238558/1‑150 (MQ=255)
                                                                                                 gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAg                                               >  2:1727370/1‑150 (MQ=255)
                                                                                                 gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAg                                               >  2:414453/1‑150 (MQ=255)
                                                                                                 gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAg                                               >  1:500612/1‑150 (MQ=255)
                                                                                                 gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAg                                               >  2:753960/1‑150 (MQ=255)
                                                                                                   gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCa                                             >  2:1633868/1‑150 (MQ=255)
                                                                                                   gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCa                                             >  1:759104/1‑150 (MQ=255)
                                                                                                   gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCa                                             >  1:580498/1‑150 (MQ=255)
                                                                                                   gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCa                                             >  1:967020/1‑150 (MQ=255)
                                                                                                   gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCa                                             >  1:476187/1‑150 (MQ=255)
                                                                                                   gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCa                                             >  2:1189134/1‑150 (MQ=255)
                                                                                                      ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCg                                         >  1:1507733/1‑151 (MQ=255)
                                                                                                      ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCg                                         >  1:1383814/1‑151 (MQ=255)
                                                                                                      ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCg                                         >  2:1188367/1‑151 (MQ=255)
                                                                                                        acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGa                                        >  1:444678/1‑150 (MQ=255)
                                                                                                        acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGa                                        >  2:1143107/1‑150 (MQ=255)
                                                                                                         cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACg                                      >  1:729618/1‑151 (MQ=255)
                                                                                                         cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACg                                      >  2:558782/1‑151 (MQ=255)
                                                                                                           caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGc                                     >  2:1410031/1‑150 (MQ=255)
                                                                                                           caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGc                                     >  2:623397/1‑150 (MQ=255)
                                                                                                           caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGc                                     >  1:107376/1‑150 (MQ=255)
                                                                                                             tGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCAt                                  >  1:561950/1‑151 (MQ=255)
                                                                                                               tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTc                                <  2:424518/151‑1 (MQ=255)
                                                                                                               tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTc                                <  1:1175613/151‑1 (MQ=255)
                                                                                                               tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTc                                <  2:291061/151‑1 (MQ=255)
                                                                                                               tCTTGCCTCCCTGCAAAGAAGAGACTAACATTAGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCAtt                                 >  2:655167/1‑150 (MQ=255)
                                                                                                                  tGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCtat                             <  1:514115/151‑1 (MQ=255)
                                                                                                                    ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATAc                           <  1:1275178/151‑1 (MQ=255)
                                                                                                                    ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATAc                           <  2:1130459/151‑1 (MQ=255)
                                                                                                                    ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATAc                           <  2:739727/151‑1 (MQ=255)
                                                                                                                    ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATAc                           <  1:701253/151‑1 (MQ=255)
                                                                                                                     cTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATAcc                          <  2:1218550/151‑1 (MQ=255)
                                                                                                                     cTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATAcc                          <  2:664176/151‑1 (MQ=255)
                                                                                                                         cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGAt                      >  1:244539/1‑151 (MQ=255)
                                                                                                                         cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGAt                      >  1:1679838/1‑151 (MQ=255)
                                                                                                                           gCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATaa                    >  1:628939/1‑151 (MQ=255)
                                                                                                                           gCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATaa                    >  1:487785/1‑151 (MQ=255)
                                                                                                                           gCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATaa                    >  2:1347987/1‑151 (MQ=255)
                                                                                                                           gCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAACCACCCGACGCCATTCTATACCTGATaa                    >  2:29511/1‑151 (MQ=255)
                                                                                                                            cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAAt                   <  1:313843/151‑1 (MQ=255)
                                                                                                                              aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAAttc                 <  2:1261219/151‑1 (MQ=255)
                                                                                                                              aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCCATACCTGATAAttc                 <  2:144606/151‑1 (MQ=255)
                                                                                                                               agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAgg                                                                                                                       >  2:1396172/1‑48 (MQ=255)
                                                                                                                               agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAgg                                                                                                                       <  1:1395550/48‑1 (MQ=255)
                                                                                                                                gaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAAttctt               >  2:1156885/1‑151 (MQ=255)
                                                                                                                                gaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAAttctt               >  2:234953/1‑151 (MQ=255)
                                                                                                                                   gagaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTcc             >  1:379651/1‑150 (MQ=255)
                                                                                                                                     gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCcg           >  1:362321/1‑150 (MQ=255)
                                                                                                                                     gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCcg           >  1:380731/1‑150 (MQ=255)
                                                                                                                                     gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCcg           >  1:92513/1‑150 (MQ=255)
                                                                                                                                      aCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCcgcg         >  2:1223511/1‑151 (MQ=255)
                                                                                                                                           cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGtt    >  1:1393560/1‑151 (MQ=255)
                                                                                                                                             ttCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCa  >  1:1224338/1‑151 (MQ=255)
                                                                                                                                             ttCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCa  >  2:692661/1‑151 (MQ=255)
                                                                                                                                             ttCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCa  >  1:715871/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ATTTAAAAGGTTCCATTAATAGCCTCCCTGTTTTTTTAGTATTATCGATGATATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCA  >  NC_000913/2887204‑2887495

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