Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,501,344 A→T 6.7% 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.7% 232.2 / 5.7 74intergenic (‑239/+16)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.98e-01

AGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGC  >  NC_002947/4501198‑4501488
                                                                                                                                                  |                                                                                                                                                
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)
                                                                                                                                                  |                                                                                                                                                
AGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACACGC  >  NC_002947/4501198‑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: