Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 634,122 A→T 5.7% Q142L (CAG→CTG)  PP_0546 → sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947634,1220AT 5.7% 216.0 / 3.4 70Q142L (CAG→CTG) PP_0546sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (26/40);  new base T (2/2);  total (28/42)
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.85e-01

GCTGACCGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATCGCGGCAAAGCCAGT  >  NC_002947/633973‑634271
                                                                                                                                                     |                                                                                                                                                      
gctgACCGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCa                                                                                                                                                        <  2:1620194/150‑1 (MQ=255)
 ctgACCGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAg                                                                                                                                                       <  1:345122/150‑1 (MQ=255)
      cGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCt                                                                                                                                                  <  1:248137/150‑1 (MQ=255)
      cGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCACCCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCt                                                                                                                                                  <  1:1261741/150‑1 (MQ=255)
      cGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCACCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCt                                                                                                                                                  <  1:789789/150‑1 (MQ=255)
       gACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGACTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTg                                                                                                                                                >  2:1896036/1‑150 (MQ=255)
       gACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGGGGCTGGGCGCTGACTGGAGCGGGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGGCGAGCGCCAGgcgccg                                                                                                                                                 >  1:1756056/1‑148 (MQ=255)
          gCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGAcc                                                                                                                                             <  2:672822/150‑1 (MQ=255)
          gCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGAcc                                                                                                                                             >  2:1089855/1‑150 (MQ=255)
           cccGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGAGTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCAGCGCCTGGCGCTCGAcca                                                                                                                                            <  2:85449/150‑1 (MQ=255)
              gcgGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGAccatcc                                                                                                                                         <  1:83202/150‑1 (MQ=255)
              gcgGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGAccatcc                                                                                                                                         >  1:1371749/1‑150 (MQ=255)
                ggggTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGAccatcctt                                                                                                                                       <  2:1918840/150‑3 (MQ=255)
                     gATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGaa                                                                                                                                  <  2:494574/150‑1 (MQ=255)
                      aTCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAAc                                                                                                                                 <  2:1843493/150‑1 (MQ=255)
                       tCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACg                                                                                                                                >  2:1371752/1‑150 (MQ=255)
                        ccTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGa                                                                                                                               <  1:358818/150‑1 (MQ=255)
                         cTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAg                                                                                                                              <  2:142093/150‑1 (MQ=255)
                         cTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAg                                                                                                                              >  2:1830945/1‑150 (MQ=255)
                          tCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTTTTGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGc                                                                                                                             >  2:268571/1‑150 (MQ=255)
                            aaTTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCAc                                                                                                                           >  2:910655/1‑150 (MQ=255)
                                gCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTcc                                                                                                                       <  1:164258/150‑1 (MQ=255)
                                gCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTcc                                                                                                                       <  1:161722/150‑1 (MQ=255)
                                   tCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTg                                                                                                                    >  1:1428526/1‑150 (MQ=255)
                                      gTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGcc                                                                                                                 <  1:1371752/150‑1 (MQ=255)
                                      gTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGcc                                                                                                                 <  2:694616/150‑1 (MQ=255)
                                      gTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGcc                                                                                                                 >  2:1554661/1‑150 (MQ=255)
                                          ccccGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGa                                                                                                                               <  2:1692273/132‑1 (MQ=255)
                                          ccccGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGa                                                                                                                               >  1:1692273/1‑132 (MQ=255)
                                                     ggcgCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGGGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGa                                                                                                  <  1:850397/150‑1 (MQ=255)
                                                            cTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGaa                                                                                                                                  <  1:1223469/111‑1 (MQ=255)
                                                            cTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGaa                                                                                                                                  >  2:1223469/1‑111 (MQ=255)
                                                            cTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTc                                                                                           <  1:808347/150‑1 (MQ=255)
                                                             ttCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGAc                                                                                                 >  2:1835728/1‑143 (MQ=255)
                                                             ttCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGAc                                                                                                 <  1:1835728/143‑1 (MQ=255)
                                                                  gCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCAGCGCCTGGCGCTCGAccagac                                                                                                                                         <  2:1050169/98‑4 (MQ=38)
                                                                  gCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCAGCGCCTGGCGCTCGAccagac                                                                                                                                         >  1:1050169/1‑95 (MQ=38)
                                                                      gCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCg                                                                                 <  1:428428/150‑1 (MQ=255)
                                                                       ccGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCgg                                                                                <  2:758857/150‑1 (MQ=255)
                                                                        cGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGgt                                                                               <  2:1728419/150‑1 (MQ=255)
                                                                         gggCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGgtg                                                                              <  2:1371749/150‑1 (MQ=255)
                                                                              gtgGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGa                                                                         >  1:2578/1‑150 (MQ=255)
                                                                               tgGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGAc                                                                        <  1:1089855/150‑1 (MQ=255)
                                                                                 gCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCTATCAGAACGAGCAATTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTCTTCGAccc                                                                      >  2:1039259/1‑150 (MQ=255)
                                                                                     ggCGCTGACTTGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCAt                                                                  <  1:454449/150‑1 (MQ=255)
                                                                                      gcgcTGACTGGAGCGAGGCGCGCGAGTGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATg                                                                 >  2:170660/1‑150 (MQ=255)
                                                                                        gcTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGc                                                               >  1:1052698/1‑150 (MQ=255)
                                                                                          tGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACGTCCGTGGCCGTCACACCGGGCTCACCTGCTCGGCCAGCCCGCTGTTCGACCCGCATGGCgg                                                             >  2:600474/1‑149 (MQ=255)
                                                                                            aCTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGAcc                                                           >  2:1415716/1‑150 (MQ=255)
                                                                                            aCTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGAcc                                                           <  1:1799990/150‑1 (MQ=255)
                                                                                             cTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACct                                                          <  2:1428526/150‑1 (MQ=255)
                                                                                                gAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTg                                                                                                                    <  1:467242/89‑1 (MQ=255)
                                                                                                gAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTg                                                                                                                    >  2:467242/1‑89 (MQ=255)
                                                                                                 aGCGGGGCGCGCGAGCGCACCAATGGCCTCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACctgctg                                                      <  1:600474/150‑1 (MQ=255)
                                                                                                    gAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCg                                                   <  1:874021/150‑1 (MQ=255)
                                                                                                         gcgcgAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCt                                              <  2:94517/150‑1 (MQ=255)
                                                                                                          cgcgAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTc                                             <  2:1861504/150‑1 (MQ=255)
                                                                                                                  cACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTc                                     <  1:910655/150‑1 (MQ=255)
                                                                                                                    ccAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGAc                                          <  2:399956/143‑1 (MQ=255)
                                                                                                                    ccAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGAc                                          >  1:399956/1‑143 (MQ=255)
                                                                                                                     cAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCctc                                  >  2:888162/1‑150 (MQ=255)
                                                                                                                          gCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTgcgc                             <  2:1862378/150‑1 (MQ=255)
                                                                                                                            aTCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTgcgcgc                           >  2:609148/1‑150 (MQ=255)
                                                                                                                             tCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCCTGGCGAACTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGcc                          >  1:679700/1‑150 (MQ=255)
                                                                                                                              cGGCACATGTCTGGTCGAGCGCCAGGCGCTGTCCATCCATCAGAACGAGCACTTCCGTGGCCGTCACCCCGGGCTGACCTGCTCTGCCAGCCTGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGccc                         <  2:450936/150‑1 (MQ=255)
                                                                                                                              cGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGccc                         <  1:673598/150‑1 (MQ=255)
                                                                                                                               ggCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCGCGACGTGTCCTCTGCGCGcccc                        >  2:259945/1‑150 (MQ=255)
                                                                                                                                       tCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATc                >  2:392900/1‑150 (MQ=255)
                                                                                                                                         tGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATcgc              <  2:2578/150‑1 (MQ=255)
                                                                                                                                          tggtcaGCGCCTGGCGCTCGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATcgc              >  2:570516/6‑150 (MQ=255)
                                                                                                                                                    cAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATCGCGGCAAAGCCAg   <  1:698823/150‑1 (MQ=255)
                                                                                                                                                     aGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATCGCGGCAAAGCCAGt  >  2:1817095/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCTGACCGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATCGCGGCAAAGCCAGT  >  NC_002947/633973‑634271

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: