Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,501,341 A→T 6.5% intergenic (‑236/+19) PP_3989 ← / ← PP_t60 DNA‑cytosine methyltransferase/tRNA‑Ser

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,501,3410AT 6.5% 232.9 / 6.7 77intergenic (‑236/+19)PP_3989/PP_t60DNA‑cytosine methyltransferase/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (39/33);  new base T (2/3);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.51e-01

AAAGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGCT  >  NC_002947/4501196‑4501489
                                                                                                                                                 |                                                                                                                                                    
aaaGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGGAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGaataat                                                                                                                                                  >  1:2093351/1‑150 (MQ=255)
    ctctAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAgg                                                                                                                                              >  1:1539686/1‑150 (MQ=255)
     tctAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAggg                                                                                                                                             >  2:806975/1‑150 (MQ=255)
     tctAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGAAACCAGATCTTAATAATCAggg                                                                                                                                             >  2:1023249/1‑150 (MQ=255)
       tAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCt                                                                                                                                           >  1:2435657/1‑150 (MQ=255)
        aGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGAAAGCAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCAGGGCtt                                                                                                                                          <  1:2126605/150‑1 (MQ=255)
         ggggAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCttt                                                                                                                                         >  1:1949410/1‑150 (MQ=255)
           ggAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCAGGGCttttt                                                                                                                                       <  2:1365662/150‑3 (MQ=255)
           ggAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGa                                                                                                                                       <  1:1734963/150‑1 (MQ=255)
               gcgcGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAatat                                                                                                                                   >  1:2219521/1‑150 (MQ=255)
                 gcgTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCAGGGCtttt                                                                                                                                        <  2:1966503/143‑2 (MQ=255)
                 gcgTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCAGGGCtttt                                                                                                                                        >  1:1966503/1‑142 (MQ=255)
                   gTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGaataat                                                                                                                                                  >  1:2172533/1‑131 (MQ=255)
                   gTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGaataat                                                                                                                                                  <  2:2172533/131‑1 (MQ=255)
                      aCGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGaa                                                                                                                            <  2:2219521/150‑1 (MQ=255)
                       cGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAg                                                                                                                           >  2:466972/1‑150 (MQ=255)
                        gCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGc                                                                                                                          <  2:2446961/150‑1 (MQ=255)
                          cAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGt                                                                                                                        >  2:2069175/1‑150 (MQ=255)
                             gCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGCAGCGTaga                                                                                                                     >  2:2351267/1‑150 (MQ=255)
                                gAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGAt                                                                                                                  >  1:2133503/1‑150 (MQ=255)
                                      gcTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAAc                                                                                                            >  2:1804196/1‑150 (MQ=255)
                                       cTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAAct                                                                                                           >  2:533521/1‑150 (MQ=255)
                                        tCNCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAANTGCGGANACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAActc                                                                                                          >  1:1229499/1‑150 (MQ=255)
                                        tCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGaa                                                                                                                            >  1:394018/1‑132 (MQ=255)
                                        tCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGaa                                                                                                                            <  2:394018/132‑1 (MQ=255)
                                          gCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTa                                                                                                        >  2:1442330/1‑150 (MQ=255)
                                          gCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTa                                                                                                        <  2:1949410/150‑1 (MQ=255)
                                            acacaGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAgg                                                                                                      >  2:285289/1‑150 (MQ=255)
                                               caGCATGGGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATa                                                                                                   <  1:199978/150‑1 (MQ=255)
                                                aGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGaa                                                                                                                                      >  2:1351000/1‑114 (MQ=255)
                                                aGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGaa                                                                                                                                      <  1:1351000/114‑1 (MQ=255)
                                                   aTGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTag                                                                                                                      >  1:376758/1‑127 (MQ=255)
                                                   aTGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTag                                                                                                                      <  2:376758/127‑1 (MQ=255)
                                                     gtgtCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGcc                                                                                             <  1:1657325/150‑1 (MQ=255)
                                                           gATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCGTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGa                                                                                       >  2:1299755/1‑150 (MQ=255)
                                                           gATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGGTAGTTGCCCATCGa                                                                                       <  1:466972/150‑1 (MQ=255)
                                                                       aTCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGa                                                                           <  2:2435657/150‑1 (MQ=255)
                                                                       aTCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGa                                                                           <  2:2133503/150‑1 (MQ=255)
                                                                        tCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCaa                                                                             <  1:135757/147‑1 (MQ=255)
                                                                        tCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCaa                                                                             >  2:135757/1‑147 (MQ=255)
                                                                         cAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGAcc                                                                         <  1:285289/150‑1 (MQ=255)
                                                                            agagTTGCGGAAAAGACCCTGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTATGATAGTTGCCCATCGACGGTTTTCAAGACCGtt                                                                      <  2:271774/150‑1 (MQ=255)
                                                                              aaTTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGc                                                                    >  1:641784/1‑150 (MQ=255)
                                                                                ttGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCt                                                                  >  2:2437157/1‑150 (MQ=255)
                                                                                    ggAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGaa                                                                                                                            <  2:2221845/88‑1 (MQ=255)
                                                                                    ggAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGaa                                                                                                                            >  1:2221845/1‑88 (MQ=255)
                                                                                    ggAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAACCAGGGCTTCGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTaaa                                                              <  1:1483057/150‑1 (MQ=255)
                                                                                        aCGACCCGGATGACAGAAGGGCAAAAACTAGAGCCTCAGAAACGCAAAAAGACCTGAATAATCCGGGCTCTGAATATGACGGAAACGTAGAGACTCGAACTCTAGGATAGTTGCCAATCGACGCTGTTCAAGAACTCGGCCTTAAACTAc                                                          >  1:225599/1‑150 (MQ=255)
                                                                                             ccGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGc                                                                                              >  1:1134408/1‑109 (MQ=255)
                                                                                             ccGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGc                                                                                              <  2:1134408/109‑1 (MQ=255)
                                                                                                gATGACAGAAGGCAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCAc                                                  <  1:2069175/150‑1 (MQ=255)
                                                                                                    aCAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCtt                                              >  2:20687/1‑150 (MQ=255)
                                                                                                     cAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTACCGCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTc                                             <  1:1023249/150‑1 (MQ=255)
                                                                                                      aGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAActc                                                                                                          <  2:1995167/88‑1 (MQ=255)
                                                                                                      aGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTcc                                            >  2:754979/1‑150 (MQ=255)
                                                                                                      aGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTcc                                            >  2:235529/1‑150 (MQ=255)
                                                                                                      aGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGATGTGAATATGGCGGAAGCGTAGAGATTCGAActc                                                                                                          >  1:1995167/1‑88 (MQ=255)
                                                                                                       gAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCa                                           <  1:20687/150‑1 (MQ=255)
                                                                                                         aGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGc                                         <  1:782089/150‑1 (MQ=255)
                                                                                                           ggAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTc                                       <  2:355953/150‑1 (MQ=255)
                                                                                                                  aCTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGc                                <  1:2351267/150‑1 (MQ=255)
                                                                                                                   cTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCg                               <  1:754979/150‑1 (MQ=255)
                                                                                                                     agagCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGgc                             >  1:1940461/1‑150 (MQ=255)
                                                                                                                         ccTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCAc                                                          >  2:2078303/1‑117 (MQ=255)
                                                                                                                         ccTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCAc                                                          <  1:2078303/117‑1 (MQ=255)
                                                                                                                         ccTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGgccgcc                         >  1:815353/1‑150 (MQ=255)
                                                                                                                           tCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCat                       >  2:687135/1‑150 (MQ=255)
                                                                                                                           tCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTCTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCat                       >  2:687100/1‑150 (MQ=255)
                                                                                                                            cAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATc                                                                                         <  2:2104918/83‑1 (MQ=255)
                                                                                                                            cAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATc                                                                                         >  1:2104918/1‑83 (MQ=255)
                                                                                                                                 aCGCAAGAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACACTAGGATAGTTGCCCATCAACGGTTTTCAAGACCGTTGCCTTAACCCCCTCGTCCACGCTTCCAGCTCGTTTTGCGGCCGCCACAATAcc                 <  2:1602010/150‑1 (MQ=37)
                                                                                                                                   gCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAAc                                                             >  2:2037154/1‑104 (MQ=255)
                                                                                                                                   gCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAAc                                                             <  1:2037154/104‑1 (MQ=255)
                                                                                                                                     aaaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAActc                                                                                                          >  1:1726036/1‑57 (MQ=255)
                                                                                                                                     aaaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAActc                                                                                                          <  2:1726036/57‑1 (MQ=255)
                                                                                                                                      aaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAAt            <  1:1804196/150‑1 (MQ=255)
                                                                                                                                       aaaGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATg           >  2:620271/1‑150 (MQ=255)
                                                                                                                                          gCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGACCTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGaaa        <  1:2366291/150‑1 (MQ=255)
                                                                                                                                            gcTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATATGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTAATTTAAACCACACAACCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACAGTAATGAAgca      >  2:562170/2‑150 (MQ=37)
                                                                                                                                               gAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGc   >  2:2011273/1‑150 (MQ=255)
                                                                                                                                                aataatCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGATCTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCCGCCACGCTTCCAGCTCGTTTTGCGGCCCCCATAATACCTTAATGAAACACGCt  <  1:562170/150‑1 (MQ=255)
                                                                                                                                                aataatCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCtt               >  2:1393733/1‑137 (MQ=255)
                                                                                                                                                aataatCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCtt               <  1:1393733/137‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AAAGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGCT  >  NC_002947/4501196‑4501489

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