Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,914,733 C→G 5.6% R33P (CGG→CCG)  PP_3454 ← transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,914,7330CG 5.6% 227.2 / 4.4 72R33P (CGG→CCG) PP_3454transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (39/29);  new base G (2/2);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GGGCCGTCAGCATGAGGATTGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGG  >  NC_002947/3914587‑3914880
                                                                                                                                                  |                                                                                                                                                   
gggCCGTCAGCATGAGGATTGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCggg                                                                                                                                                  <  1:875845/150‑1 (MQ=255)
     gTCAGCATGAGGATTGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTAAAGGATGACCAGGTCAGGACTTTGATGCGGGTATGCGGACACCACAAGAGCTCCCCGGGCGATc                                                                                                                                             >  2:674462/1‑150 (MQ=255)
      tCAGCATGAGGATTGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCa                                                                                                                                            >  2:168259/1‑150 (MQ=255)
       cAGCATGAGGATTGGCAGGGTCTGCGACTCCTGGCGCAGCAGACGCATTAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCGGGGCGATCat                                                                                                                                           <  1:966103/150‑2 (MQ=255)
          cATGAGGATTGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCGGGGCGATCatgtg                                                                                                                                        <  1:423753/150‑5 (MQ=255)
          cATGAGGATTGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTg                                                                                                                                        >  1:917726/1‑150 (MQ=255)
            tGAGGATTGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGa                                                                                                                                      <  1:168259/150‑1 (MQ=255)
             gAGGATTGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGAc                                                                                                                                     >  1:932455/1‑150 (MQ=255)
                  ttGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGa                                                                                                                                >  1:1021893/1‑150 (MQ=255)
                  ttGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGa                                                                                                                                >  2:485007/1‑150 (MQ=255)
                     gCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGt                                                                                                                             >  2:656248/1‑150 (MQ=255)
                      cAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTc                                                                                                                            <  1:875110/150‑1 (MQ=255)
                       aGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTcg                                                                                                                           >  2:908376/1‑150 (MQ=255)
                        gggTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTcgt                                                                                                                          >  1:548477/1‑150 (MQ=255)
                         ggTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGAGCACCTGGACATGGAAGTcgtt                                                                                                                         >  1:19571/1‑150 (MQ=255)
                                 aCTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGGACCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCa                                                                                                                 <  1:674462/150‑1 (MQ=255)
                                 aCTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCa                                                                                                                 <  1:797570/150‑1 (MQ=255)
                                  cTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAg                                                                                                                <  1:885417/150‑1 (MQ=255)
                                   tCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAgg                                                                                                               <  1:485007/150‑1 (MQ=255)
                                       ggCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGc                                                                                                           >  2:423852/1‑150 (MQ=255)
                                             gcagcCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGa                                                                                                     >  2:935863/1‑150 (MQ=255)
                                              cagcCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGAt                                                                                                    <  1:423852/150‑1 (MQ=255)
                                                gcCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCa                                                                                                  <  2:1021893/150‑1 (MQ=255)
                                                 ccGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAg                                                                                                 <  2:917726/150‑1 (MQ=255)
                                                      acaGCTGCAGGCCATCGATGCCCGGCAGCATCCGTTCGAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCgg                                                                                            >  1:932319/1‑150 (MQ=255)
                                                      acaGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGGCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCgg                                                                                            >  1:977157/1‑150 (MQ=255)
                                                           tGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCCGGTCCGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAgg                                                                                       >  1:22249/1‑150 (MQ=255)
                                                              aGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGCCAGGgcg                                                                                    >  1:391794/1‑150 (MQ=255)
                                                              aGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCGGGCGATCAGCTCGGACAGGgcg                                                                                    >  2:557870/1‑150 (MQ=255)
                                                              aGGCCATCGATGCCCGGCAGCATCAGGTCAAGGAACGCCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCACCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGgcg                                                                                    >  1:933444/1‑150 (MQ=255)
                                                                gCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGa                                                                                                                                      <  2:1075908/98‑1 (MQ=255)
                                                                gCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGa                                                                                                                                      >  1:1075908/1‑98 (MQ=255)
                                                                accATCGATGCACGACAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCt                                                                                  >  2:42465/2‑150 (MQ=255)
                                                                         gCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTc                                                                         <  1:42465/150‑1 (MQ=255)
                                                                           ccGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGa                                                                       >  1:526450/1‑150 (MQ=255)
                                                                             ggCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATgg                                                                                                                                 >  1:114631/1‑90 (MQ=255)
                                                                             ggCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATgg                                                                                                                                 <  2:114631/90‑1 (MQ=255)
                                                                              gcagcaTCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCa                                                                    <  2:932319/150‑1 (MQ=255)
                                                                               gatcttCAGGTCAAGGATGACCAGGTCAGGACTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAg                                                                   <  2:261898/145‑1 (MQ=255)
                                                                                   aTCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGa                                                                       <  2:992193/142‑1 (MQ=255)
                                                                                   aTCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGa                                                                       >  1:992193/1‑142 (MQ=255)
                                                                                    tCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACTTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGaaaa                                                              >  1:354488/1‑150 (MQ=255)
                                                                                     cAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAAt                                                             <  1:291507/150‑1 (MQ=255)
                                                                                     cAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGAAAGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGagatc                                                             >  1:261898/1‑146 (MQ=255)
                                                                                           aaGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCgg                                                       >  1:95288/1‑150 (MQ=255)
                                                                                           aaGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCgg                                                       >  1:373602/1‑150 (MQ=255)
                                                                                           aaGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAAGGTTCgg                                                       <  1:367432/150‑1 (MQ=255)
                                                                                                  aCCAGCTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATCGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCCACCACGAAATGTTCGCCCTACGTg                                               >  2:332169/1‑150 (MQ=255)
                                                                                                   ccAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTg                                               >  2:727940/1‑150 (MQ=255)
                                                                                                   ccAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTg                                               >  2:590660/1‑150 (MQ=255)
                                                                                                       gTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTg                                           <  1:727940/150‑1 (MQ=255)
                                                                                                              cTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCt                                    <  2:373602/150‑1 (MQ=255)
                                                                                                                     gACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCa                                                                                                                 <  2:1072567/66‑1 (MQ=255)
                                                                                                                     gACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCa                                                                                                                 >  1:1072567/1‑66 (MQ=255)
                                                                                                                      aCGGTACTCGTCCCGAAACTAATCGCACCAGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGTCAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCAAAGGTGCTTGGGCGTCTGGGTTGaa                            >  2:562174/1‑150 (MQ=255)
                                                                                                                      aCGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGaa                            <  2:1025987/150‑1 (MQ=255)
                                                                                                                      aCGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGaa                            <  2:251142/150‑1 (MQ=255)
                                                                                                                      aCGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCAACTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGaa                            <  1:1081719/150‑1 (MQ=255)
                                                                                                                        ggTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGCACCTGGAAGTCGTTGCGTTGCAGGTCGCTGGCGATCAGCTCGGACAGGGCGCTGTTGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTGTGGCGTCCGAGTTGAAGt                          >  2:607478/1‑150 (MQ=255)
                                                                                                                               gTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTgtc                                                                               <  2:44792/90‑1 (MQ=255)
                                                                                                                               gTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTgtc                                                                               >  1:44792/1‑90 (MQ=255)
                                                                                                                               gTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCtg                   >  2:984688/1‑150 (MQ=255)
                                                                                                                                  cAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTc                                                                                              >  1:582778/1‑72 (MQ=255)
                                                                                                                                  cAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTc                                                                                              <  2:582778/72‑1 (MQ=255)
                                                                                                                                  cAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTgg                <  2:22249/150‑1 (MQ=255)
                                                                                                                                     caggtGATCGCCCGGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGa             >  2:777727/5‑150 (MQ=255)
                                                                                                                                     cacaTGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCGTAGGTGGTGAGGCGGCTGGGGTGAAGGAAGCCTGTGGCGa             >  2:121565/1‑150 (MQ=255)
                                                                                                                                      acaTGATCGCCCCGGGCGATCACTTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAg            <  2:354488/150‑1 (MQ=255)
                                                                                                                                      acaTGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAg            >  2:804157/1‑150 (MQ=255)
                                                                                                                                          gATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCt        >  2:111702/1‑150 (MQ=255)
                                                                                                                                           aTCGCCCGGGGCGATCACCTGGACATGGAAGTCGTAGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCtt       >  2:593569/1‑150 (MQ=255)
                                                                                                                                            tCGCCCCGGGCGATCCAATACACATGGAAGTCGTTGCGTTGCAGGTAGCCGGCGATCAGCCCGGACAGGGCGCTGTCGTCTTCGCCCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTc      <  1:607478/150‑1 (MQ=255)
                                                                                                                                                cccGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCgg  >  2:160849/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGGCCGTCAGCATGAGGATTGGCAGGGTCTGCGACTCCTGGCGCAGCAGCCGGCACAGCTGCAGGCCATCGATGCCCGGCAGCATCAGGTCAAGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGG  >  NC_002947/3914587‑3914880

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