Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,804,021 C→G 8.2% intergenic (‑105/‑39) PP_3361 ← / → PP_3362 aminotransferase/AsnC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,804,0210CG 8.2% 181.1 / 7.2 63intergenic (‑105/‑39)PP_3361/PP_3362aminotransferase/AsnC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (21/35);  new base G (2/3);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-02

GCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAGTGC  >  NC_002947/3803879‑3804166
                                                                                                                                               |                                                                                                                                                 
gCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCTGACCgg                                                                                                                                             >  5:323787/1‑149 (MQ=255)
gCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCg                                                                                                                                              <  3:327009/148‑1 (MQ=255)
 ccGAAAGCCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCTGAGCGCGCGGATTGTTTTTCAATAGCTTGTCTTGCGCCGCCGAAAGGGACCCCCTGCATAGGGTCGGGCGCAACAGCCATTATGCCGGTTCCGCCCgg                                                                                                                                             >  1:290435/1‑149 (MQ=255)
 ccGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCgg                                                                                                                                             >  3:117405/1‑148 (MQ=255)
 ccGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGc                                                                                                                                            <  1:149772/149‑1 (MQ=255)
   gAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTCGGAACCGGCAt                                                                                                                                          <  2:388556/149‑1 (MQ=255)
     aaaCCGAGGTGCCGTGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGGTTGTCATGCGCCGCCGAAACGGACACCCTGCATTCGGTCGGGCGCAACAGGCATTATGCCGGTCGGAACCGGCATa                                                                                                                                         <  6:6642/148‑1 (MQ=255)
       aCCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTCGGAACCGGCATAATg                                                                                                                                      <  4:159488/149‑1 (MQ=255)
         cGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGAc                                                                                                                                    >  2:159147/1‑149 (MQ=255)
          gAGGTGCCGGGTGGGAAAGGGCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACg                                                                                                                                   <  4:127467/149‑1 (MQ=255)
               gCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAg                                                                                                                              <  6:106146/149‑1 (MQ=255)
               gCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAg                                                                                                                              <  2:9406/149‑1 (MQ=255)
                          aaGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAgg                                                                                                                    <  8:195347/148‑1 (MQ=255)
                                    atatCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGc                                                                                                         >  1:103317/1‑149 (MQ=255)
                                       tcatcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCTGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGt                                                                                                       <  6:323787/148‑1 (MQ=255)
                                       tcatcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGt                                                                                                       <  2:220240/148‑1 (MQ=255)
                                        catcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGt                                                                                                       >  5:21437/1‑147 (MQ=255)
                                        catcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGt                                                                                                       <  6:21437/147‑1 (MQ=255)
                                         atcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGa                                                                                                     <  6:34738/148‑1 (MQ=255)
                                         atcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGa                                                                                                     >  5:34738/1‑148 (MQ=255)
                                          tcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTACAATGCGTGAcc                                                                                                   >  7:132067/1‑149 (MQ=255)
                                            aGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTg                                                                                                 >  2:20291/1‑149 (MQ=255)
                                            aGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTg                                                                                                 <  1:171625/149‑1 (MQ=255)
                                            aGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTg                                                                                                 <  1:165893/149‑1 (MQ=255)
                                              gTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGa                                                                                               <  1:370649/149‑1 (MQ=255)
                                                         cgcTGATTGTTTTTCAATAGCTTGTCATGCCCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTa                                                                                          >  1:33933/1‑144 (MQ=255)
                                                          gcTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATc                                                                                    <  2:103317/149‑1 (MQ=255)
                                                              aTTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATcgc                                                                                  <  8:147246/147‑1 (MQ=255)
                                                              aTTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGATGCAATGCGTGACCTGGACGCTAAAGATcgcg                                                                                 <  6:312208/148‑1 (MQ=255)
                                                                       cAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTg                                                                        <  1:97377/148‑1 (MQ=255)
                                                                        aaTAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATaa                                                                                                                                        >  7:162460/1‑83 (MQ=255)
                                                                        aaTAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATaa                                                                                                                                        <  8:162460/83‑1 (MQ=255)
                                                                        aaTAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAg                                                                      >  5:104559/1‑149 (MQ=255)
                                                                            gCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTAc                                                                  <  8:182025/149‑1 (MQ=255)
                                                                             cTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACt                                                                 >  1:344917/1‑149 (MQ=255)
                                                                                 tCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGAATTATGCCGGTTCCGACCGGGATAATGACGGCAAGCTTATCTAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTc                                                              <  8:226269/148‑1 (MQ=255)
                                                                                   aTGCGCCGCCGAAACGGACACCCTGCATACGTTCGGGCGCAACAGTCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGTTAAAGATCGCGAGATTCTTGAGGTACTTTCGaa                                                           <  1:249682/149‑1 (MQ=255)
                                                                                         cgccgAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGAcgc                                                     >  4:381191/1‑149 (MQ=255)
                                                                                         cgccgAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGAcgc                                                     <  7:79052/149‑1 (MQ=255)
                                                                                           ccgAAAAGGAGACCCTGCATACGGTCGGGCGCAACATGCATTATGCCGGTTCCGACCGGGATATTGACGGCAAGCTTATCTAGGGGCTGGAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGAcgcgc                                                   <  4:237616/149‑1 (MQ=255)
                                                                                                      acCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTc                                        <  2:290435/149‑1 (MQ=255)
                                                                                                           gCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTAtt                                              >  4:113588/1‑138 (MQ=255)
                                                                                                           gCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTAtt                                              <  3:113588/138‑1 (MQ=255)
                                                                                                              tACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGtt                                >  2:222948/1‑149 (MQ=255)
                                                                                                               aCGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGt                                 >  8:48455/1‑147 (MQ=255)
                                                                                                                cGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTa                               <  4:117405/148‑1 (MQ=255)
                                                                                                                cGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAg                              <  6:227725/149‑1 (MQ=255)
                                                                                                                 ggTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAgc                             >  5:370723/1‑149 (MQ=255)
                                                                                                                   tCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAgccg                           >  5:219692/1‑149 (MQ=255)
                                                                                                                           aaCAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCa                     >  5:153946/1‑147 (MQ=255)
                                                                                                                           aaCAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATa                   <  1:20291/149‑1 (MQ=255)
                                                                                                                             cAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCATGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGGGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAgg                 <  8:184405/149‑1 (MQ=255)
                                                                                                                             cAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAgg                 <  1:159147/149‑1 (MQ=255)
                                                                                                                             cAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTATCAAGGGGCTTCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGc                <  2:102391/149‑1 (MQ=255)
                                                                                                                             cAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATTGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGc                <  4:72889/149‑1 (MQ=255)
                                                                                                                                tcATTATGCCGGTCGGAACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCt              >  6:261362/2‑149 (MQ=255)
                                                                                                                                tcATTATGCCGGTCGGAACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGAGTGACCTGTACGCTAAATATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCt              >  8:284470/2‑149 (MQ=255)
                                                                                                                                   ttATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTc           <  6:219692/149‑1 (MQ=255)
                                                                                                                                    tATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTcc          >  2:274583/1‑149 (MQ=255)
                                                                                                                                         cGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTgcg                                           >  1:132647/1‑111 (MQ=255)
                                                                                                                                         cGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTgcg                                           <  2:132647/111‑1 (MQ=255)
                                                                                                                                          ggTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAGt    <  2:33933/149‑1 (MQ=255)
                                                                                                                                            ttCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTACAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAGTGc  <  8:132067/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAGTGC  >  NC_002947/3803879‑3804166

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