Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,501,344 A→T 6.5% intergenic (‑239/+16) 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,3440AT 6.5% 226.3 / 6.7 77intergenic (‑239/+16)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 = 3.85e-01

AAAGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGCTGT  >  NC_002947/4501196‑4501491
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                                                    <  2:2172533/131‑1 (MQ=255)
                   gTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGaataat                                                                                                                                                    >  1:2172533/1‑131 (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                                                                                                                        <  2:376758/127‑1 (MQ=255)
                                                   aTGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTag                                                                                                                        >  1:376758/1‑127 (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:2133503/150‑1 (MQ=255)
                                                                       aTCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGa                                                                             <  2:2435657/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                                                                                                                              >  1:2221845/1‑88 (MQ=255)
                                                                                    ggAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGaa                                                                                                                              <  2:2221845/88‑1 (MQ=255)
                                                                                    ggAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAACCAGGGCTTCGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTaaa                                                                <  1:1483057/150‑1 (MQ=255)
                                                                                        aCGACCCGGATGACAGAAGGGCAAAAACTAGAGCCTCAGAAACGCAAAAAGACCTGAATAATCCGGGCTCTGAATATGACGGAAACGTAGAGACTCGAACTCTAGGATAGTTGCCAATCGACGCTGTTCAAGAACTCGGCCTTAAACTAc                                                            >  1:225599/1‑150 (MQ=255)
                                                                                             ccGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGc                                                                                                <  2:1134408/109‑1 (MQ=255)
                                                                                             ccGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGc                                                                                                >  1:1134408/1‑109 (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:235529/1‑150 (MQ=255)
                                                                                                      aGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTcc                                              >  2:754979/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                                                                                                            <  2:1726036/57‑1 (MQ=255)
                                                                                                                                     aaaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAActc                                                                                                            >  1:1726036/1‑57 (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                 <  1:1393733/137‑1 (MQ=255)
                                                                                                                                                aataatCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCtt                 >  2:1393733/1‑137 (MQ=255)
                                                                                                                                                  taatCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGCTGt  >  2:677695/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAAGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGCTGT  >  NC_002947/4501196‑4501491

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: