Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,501,341 A→T 6.7% 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.7% 233.3 / 6.4 74intergenic (‑236/+19)PP_3989/PP_t60DNA‑cytosine methyltransferase/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (27/42);  new base T (2/3);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GAAAAGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGC  >  NC_002947/4501194‑4501488
                                                                                                                                                   |                                                                                                                                                   
gAAAAGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGaata                                                                                                                                                   >  1:1342619/1‑150 (MQ=255)
    aGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCa                                                                                                                                               <  2:2982608/150‑1 (MQ=255)
       tctAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGGACACCGACCGTCTGGTTCTTGTAAACTATCAGAGAGGTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAggg                                                                                                                                            >  1:835681/1‑150 (MQ=255)
       tctAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAggg                                                                                                                                            <  2:2631103/150‑1 (MQ=255)
       tctAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAggg                                                                                                                                            <  1:680853/150‑1 (MQ=255)
         tAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCGAAAAGCCCTGAATAATCAGGGCt                                                                                                                                          >  1:8273/1‑150 (MQ=255)
           ggggAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCttt                                                                                                                                        <  1:945227/150‑1 (MQ=255)
              gAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGaa                                                                                                                                     <  2:398383/150‑1 (MQ=255)
              gAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGaa                                                                                                                                     >  1:77349/1‑150 (MQ=255)
                aGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAata                                                                                                                                   <  1:553191/150‑1 (MQ=255)
                            cAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGt                                                                                                                       >  1:1533075/1‑150 (MQ=255)
                             aaGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTa                                                                                                                      >  1:894012/1‑150 (MQ=255)
                             aaGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTa                                                                                                                      <  1:2937032/150‑1 (MQ=255)
                             aaGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGAGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTa                                                                                                                      >  2:1991857/1‑150 (MQ=255)
                                   aaGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGAtt                                                                                                                <  1:768611/150‑1 (MQ=255)
                                   aaGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGAtt                                                                                                                <  1:1001052/150‑1 (MQ=255)
                                            gCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTa                                                                                                       >  1:1866126/1‑150 (MQ=255)
                                             cacacaGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAg                                                                                                      <  1:1170665/150‑1 (MQ=255)
                                             cacacaGCATGTGACGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCGGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAg                                                                                                      <  2:1118028/150‑1 (MQ=255)
                                                acaGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGGCGAAGCGGAGAGATTCGAACTCTAGGAt                                                                                                   >  2:2689653/1‑150 (MQ=255)
                                                 caGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATa                                                                                                  <  1:1828143/150‑1 (MQ=255)
                                                   gCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGt                                                                                                >  2:2620753/1‑150 (MQ=255)
                                                    cATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGtt                                                                                               <  2:1069946/150‑1 (MQ=255)
                                                     aTGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCttt                                                                                                                                        <  1:1576413/108‑1 (MQ=255)
                                                     aTGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCttt                                                                                                                                        >  2:1576413/1‑108 (MQ=255)
                                                     aTGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTg                                                                                              <  1:563960/150‑1 (MQ=255)
                                                      tgtgTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGc                                                                                             <  2:835681/150‑1 (MQ=255)
                                                        tgtCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGc                                                                                             <  2:2947690/148‑1 (MQ=255)
                                                        tgtCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGc                                                                                             >  1:2947690/1‑148 (MQ=255)
                                                         gtCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGAAAAAATCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCa                                                                                          <  1:1458639/150‑1 (MQ=255)
                                                           cGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATc                                                                                        <  1:448285/150‑1 (MQ=255)
                                                            ggATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCAGGGCtttt                                                                                                                                       >  2:1726701/1‑101 (MQ=255)
                                                            ggATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCAGGGCtttt                                                                                                                                       <  1:1726701/102‑2 (MQ=255)
                                                                 cGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGt                                                                                  >  1:2029224/1‑150 (MQ=255)
                                                                  ggAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGtt                                                                                 >  2:464373/1‑150 (MQ=255)
                                                                   gAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCAggg                                                                                                                                            >  2:744250/1‑90 (MQ=255)
                                                                   gAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCAggg                                                                                                                                            <  1:744250/90‑1 (MQ=255)
                                                                    aaaCGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGtttt                                                                               >  2:1789133/1‑150 (MQ=255)
                                                                     aaCGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTc                                                                              <  2:2029224/150‑1 (MQ=255)
                                                                       cGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCaa                                                                            <  2:1342619/150‑1 (MQ=255)
                                                                       cGATCAGAGAATTAGGGAAACGACCCGGATTACAGAAGGGAAAAAACTAGAGCCTCATAAACGCAAAAAGCTCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTCGAGATTCGAACTCTAGGATAGTTGCCCATCGATGGTGTTCaa                                                                            <  2:128221/150‑1 (MQ=37)
                                                                             gagaATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGt                                                                      >  2:611316/1‑150 (MQ=255)
                                                                                aaTTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGCCCGTTGc                                                                   <  1:1789133/150‑1 (MQ=255)
                                                                                aaTTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGc                                                                   <  1:2620753/150‑1 (MQ=255)
                                                                                 aTTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGcc                                                                  <  2:22722/150‑1 (MQ=255)
                                                                                  ttGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCttt                                                                                                                                        >  2:267184/1‑79 (MQ=255)
                                                                                  ttGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCttt                                                                                                                                        <  1:267184/79‑1 (MQ=255)
                                                                                    gCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCGCATCGACGGTTTTCAAGACCGTTGCCTTa                                                               >  2:2786058/1‑150 (MQ=255)
                                                                                                       cAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTc                                            <  1:464373/150‑1 (MQ=255)
                                                                                                        aGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGaa                                                                                                            >  2:1935345/1‑85 (MQ=255)
                                                                                                        aGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGaa                                                                                                            <  1:1935345/85‑1 (MQ=255)
                                                                                                          aaGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTg                                                                    <  2:327494/123‑1 (MQ=255)
                                                                                                          aaGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTg                                                                    >  1:327494/1‑123 (MQ=255)
                                                                                                              gAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCg                                     >  1:2640788/1‑150 (MQ=255)
                                                                                                              gAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCg                                     >  2:2691911/1‑150 (MQ=255)
                                                                                                               aaaaaaCTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGaa                                                                                                                                     <  1:1906746/53‑1 (MQ=255)
                                                                                                               aaaaaaCTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGaa                                                                                                                                     >  2:1906746/1‑53 (MQ=255)
                                                                                                                 aaaaCTAGAGTCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGttt                                  <  1:604229/150‑1 (MQ=255)
                                                                                                                          gCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGgccgc                         <  2:2640788/150‑1 (MQ=255)
                                                                                                                           ccTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGgccgc                         >  1:2277759/1‑149 (MQ=255)
                                                                                                                           ccTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGgccgc                         <  2:2277759/149‑1 (MQ=255)
                                                                                                                             tCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCat                      <  1:2691911/150‑1 (MQ=255)
                                                                                                                               aGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCataa                    >  1:1599305/1‑150 (MQ=255)
                                                                                                                                  aaCGCGAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATAc                 <  2:8273/150‑1 (MQ=255)
                                                                                                                                  aaCGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTATAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATAc                 <  1:389183/150‑1 (MQ=255)
                                                                                                                                       tcaaaGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAGCCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTaa            <  2:1448066/148‑1 (MQ=255)
                                                                                                                                       aaaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTGGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTaa            <  1:1991857/150‑1 (MQ=255)
                                                                                                                                       aaaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGc                                                                   >  1:2247899/1‑95 (MQ=255)
                                                                                                                                       aaaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGc                                                                   <  2:2247899/95‑1 (MQ=255)
                                                                                                                                          aaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGa         <  2:2525899/150‑1 (MQ=255)
                                                                                                                                             cccTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGTCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGgccgc                         <  1:1173980/131‑1 (MQ=255)
                                                                                                                                             cccTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGTCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGgccgc                         >  2:1173980/1‑131 (MQ=255)
                                                                                                                                              ccTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAaca     >  1:1457707/1‑150 (MQ=255)
                                                                                                                                                tGAATAATCAGGGCTTCGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTACGGCCGCCATAATACCTTAATGAAACACg   <  2:77349/150‑1 (MQ=255)
                                                                                                                                                 gAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGc  <  2:1457707/150‑1 (MQ=255)
                                                                                                                                                  aataatCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTc                                                       >  1:183701/1‑96 (MQ=255)
                                                                                                                                                  aataatCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTc                                                       <  2:183701/96‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAAAAGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGC  >  NC_002947/4501194‑4501488

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