Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,603,124 G→T 5.6% W13L (TGG→TTG)  PP_5513 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,603,1240GT 5.6% 338.0 / 3.5 105W13L (TGG→TTG) PP_5513hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (50/49);  new base T (3/3);  total (53/52)
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.45e-01

Only 100 of 107 total aligned reads displayed.

GAGCGATAAGCGCTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGTCCGTTCTGGTCACTGATT  >  NC_002947/2602981‑2603272
                                                                                                                                                |                                                                                                                                                    
gAGCGATAAGCGCTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGtggatgga                                                                                                                                                <  1:206647/150‑1 (MQ=255)
 aGCGATAAGCGCTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATggag                                                                                                                                               <  2:901484/150‑1 (MQ=255)
 aGCGATAAGCGCTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATggag                                                                                                                                               <  2:1559966/150‑1 (MQ=255)
   cGATTAGCGCTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCGGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCCACTTGATggagaa                                                                                                                                             <  1:248132/150‑3 (MQ=255)
    gATAAGCGCTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGATGAACACTTCTCCATCAAGTGGATGGAGGaa                                                                                                                                            <  2:1405249/150‑1 (MQ=255)
      tAAGCGCTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGt                                                                                                                                          <  2:198611/150‑1 (MQ=255)
         gcgcTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGGCGCCTGGCCTCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCCACTTGATGGAGAAGTgttc                                                                                                                                      <  2:682787/150‑4 (MQ=255)
          cgcTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGAt                                                                                                                                      <  2:942637/150‑1 (MQ=255)
          cgcTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGAt                                                                                                                                      <  2:1653919/150‑1 (MQ=255)
                    gTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGact                                                                                                                            >  2:1353027/1‑150 (MQ=255)
                         ccGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTgg                                                                                                                       >  1:1707408/1‑150 (MQ=255)
                         ccGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTgg                                                                                                                       >  1:1063112/1‑150 (MQ=255)
                          cGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTggg                                                                                                                      <  1:1676969/150‑1 (MQ=255)
                                      cccGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGt                                                                                                                                 >  2:1930358/1‑127 (MQ=255)
                                      cccGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGt                                                                                                                                 <  1:1930358/127‑1 (MQ=255)
                                        cGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTc                                                                                                        <  1:97157/150‑1 (MQ=255)
                                         gAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGa                                                                                                                                       <  2:403117/118‑1 (MQ=255)
                                         gAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGa                                                                                                                                       <  2:410903/118‑1 (MQ=255)
                                         gAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGa                                                                                                                                       >  1:410903/1‑118 (MQ=255)
                                         gAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGa                                                                                                                                       >  1:403117/1‑118 (MQ=255)
                                           gTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGc                                                                                                     <  1:1053134/150‑1 (MQ=255)
                                              cAGGCTAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCCCCTATCCTCTAc                                                                                                  >  2:1996119/1‑150 (MQ=255)
                                              cAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTAc                                                                                                  >  2:1473004/1‑150 (MQ=255)
                                               aGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGc                                                                                                                                   <  2:1063894/116‑1 (MQ=255)
                                               aGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGc                                                                                                                                   >  1:1063894/1‑116 (MQ=255)
                                               aGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACCCATGGCAGCAATACGGAACAGAACAATAAGGGGCTGGTTTTACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACt                                                                                                 >  1:45912/1‑150 (MQ=255)
                                                 gCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGa                                                                                                                              >  2:381402/1‑119 (MQ=255)
                                                 gCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGa                                                                                                                              <  1:381402/119‑1 (MQ=255)
                                                  ccAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGa                                                                                              <  2:1756731/150‑1 (MQ=255)
                                                   cAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGAt                                                                                             <  1:1353027/150‑1 (MQ=255)
                                                   cAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGAt                                                                                             >  1:2352427/1‑150 (MQ=255)
                                                    aaGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATc                                                                                            >  1:1189903/1‑150 (MQ=255)
                                                          ctccGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTccc                                                                                      >  2:2046079/1‑150 (MQ=255)
                                                          ctccGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTccc                                                                                      >  2:2271216/1‑150 (MQ=255)
                                                           tccGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTTACCGACGAACACTTCTCCATCAAGCGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGCTCGGTCCCg                                                                                     >  1:214460/1‑150 (MQ=255)
                                                           tccGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCg                                                                                     >  1:2449863/1‑150 (MQ=255)
                                                              gTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGc                                                                                                                     <  2:2405596/115‑1 (MQ=255)
                                                              gTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGc                                                                                                                     >  1:2405596/1‑115 (MQ=255)
                                                               ttCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAgg                                                                                 <  1:388719/150‑1 (MQ=255)
                                                                tCCAACGTCTTCTGACACCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAggg                                                                                >  2:1223780/1‑150 (MQ=255)
                                                                 ccAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGa                                                                               >  2:65221/1‑150 (MQ=255)
                                                                 ccAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGa                                                                               <  2:1707408/150‑1 (MQ=255)
                                                                  cAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAg                                                                                  >  1:1583137/1‑146 (MQ=255)
                                                                  cAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAg                                                                                  <  2:1583137/146‑1 (MQ=255)
                                                                   aaCGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACt                                                                             >  2:1013724/1‑150 (MQ=255)
                                                                      gTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGc                                                                          <  2:1568442/150‑1 (MQ=255)
                                                                         ttctGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCgg                                                                       >  2:778932/1‑150 (MQ=255)
                                                                         ttctGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCgg                                                                       >  2:1816166/1‑150 (MQ=255)
                                                                          tctGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCCTCCATCCACTTGATGGAGAAGTgttc                                                                                                                                      >  1:106095/1‑83 (MQ=21)
                                                                          tctGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCCTCCATCCACTTGATGGAGAAGTgttc                                                                                                                                      <  2:106095/86‑4 (MQ=21)
                                                                              aCGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGt                                                                  <  1:65221/150‑1 (MQ=255)
                                                                               cGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTg                                                                 <  1:475980/150‑1 (MQ=255)
                                                                                gcgtGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGAGGATAGTGAGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTCCGTgg                                                                >  1:283260/4‑150 (MQ=255)
                                                                                gCCTGGCCGCAATAAGGAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTAc                                                                                                  >  1:254642/1‑116 (MQ=255)
                                                                                  cTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGCCCACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCCTATTCGCGTATCCGCTACTGGATCGCTCCGGAAGGTACTGGCCGGTACGTGGAt                                                              <  1:1996119/150‑1 (MQ=255)
                                                                                   tGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATg                                                             >  1:1805071/1‑150 (MQ=255)
                                                                                      ccGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCt                                                                                                                 <  2:329265/95‑1 (MQ=255)
                                                                                      ccGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCt                                                                                                                 >  1:329265/1‑95 (MQ=255)
                                                                                      ccGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTcc                                                          <  1:914131/150‑1 (MQ=255)
                                                                                       cGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGAt                                                                                                                                      >  2:806908/1‑73 (MQ=255)
                                                                                       cGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGAt                                                                                                                                      <  1:806908/73‑1 (MQ=255)
                                                                                       cGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTccc                                                         >  2:7638/1‑150 (MQ=255)
                                                                                        gCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTcccc                                                        >  1:1852045/1‑150 (MQ=255)
                                                                                          aaTAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGt                                                                  >  2:1140301/1‑138 (MQ=255)
                                                                                          aaTAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGt                                                                  <  1:1140301/138‑1 (MQ=255)
                                                                                            tAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTg                                                    >  2:2302121/1‑150 (MQ=255)
                                                                                              aGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGAGAATGTAGGAAGTGGAAGGCGTGCACTACGTGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCGGAAGAGACTGACCGCTACGTGTATGTCCCCGCTgcc                                                  >  2:742289/1‑148 (MQ=255)
                                                                                               gAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCt                                                 >  2:2202922/1‑150 (MQ=255)
                                                                                               gAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCt                                                 >  1:1274919/1‑150 (MQ=255)
                                                                                                         aaTGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGa                                       >  1:2367070/1‑150 (MQ=255)
                                                                                                         aaTGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGa                                       <  1:2271216/150‑1 (MQ=255)
                                                                                                         aaTGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGa                                       <  1:1814523/150‑1 (MQ=255)
                                                                                                          aTGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGAt                                      <  1:2302121/150‑1 (MQ=255)
                                                                                                          aTGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGAt                                      <  1:2046079/150‑1 (MQ=255)
                                                                                                              gCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAg                                                                                  <  1:945544/102‑1 (MQ=255)
                                                                                                                 gctgcACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGAt                                                              <  2:1341065/115‑1 (MQ=255)
                                                                                                                 gCTGCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGAt                                                              >  1:1341065/1‑119 (MQ=255)
                                                                                                                  cTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTc                              <  1:1816166/150‑1 (MQ=255)
                                                                                                                  cTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTc                              <  1:1013724/150‑1 (MQ=255)
                                                                                                                     cACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGTCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGt                           <  1:1535206/150‑1 (MQ=255)
                                                                                                                        cgacgaACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTAcc                        <  1:1421667/150‑1 (MQ=255)
                                                                                                                          acgaACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACt                                         >  2:713228/1‑131 (MQ=255)
                                                                                                                          acgaACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACt                                         <  1:713228/131‑1 (MQ=255)
                                                                                                                             aaCACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGAATACCGGGCTCCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCccgtc                   <  2:2276054/150‑1 (MQ=255)
                                                                                                                              acacTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTAc                                                                                                  <  2:2183766/70‑1 (MQ=255)
                                                                                                                                acTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTAc                                                                                                  <  2:1257037/68‑1 (MQ=255)
                                                                                                                                acTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTAc                                                                                                  >  1:1257037/1‑68 (MQ=255)
                                                                                                                                acTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTg                                                                 >  2:65624/1‑101 (MQ=255)
                                                                                                                                acTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTg                                                                 <  1:65624/101‑1 (MQ=255)
                                                                                                                                acTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCccgt                    <  2:1915876/146‑1 (MQ=255)
                                                                                                                                acTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCccgt                    >  1:1915876/1‑146 (MQ=255)
                                                                                                                                  ttcctcCATCCACTTGATGGAGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGTCCGTTc              >  1:432154/4‑150 (MQ=255)
                                                                                                                                  ttCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGTCCGTTc              >  2:2375982/1‑150 (MQ=255)
                                                                                                                                     ctcCATCCACTTGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGTCCGTTCTg            >  1:765674/1‑150 (MQ=255)
                                                                                                                                     ctcCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGTCCGTTCTg            <  1:1473004/150‑1 (MQ=255)
                                                                                                                                     ctcCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGCCCGTTCTg            <  2:1852045/150‑1 (MQ=255)
                                                                                                                                         aTCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGTCCGTTCTGGTCa        >  2:1586633/1‑150 (MQ=255)
                                                                                                                                            aaGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGTCCGTTCTGGTCACTg     <  1:2375982/150‑1 (MQ=255)
                                                                                                                                              gTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGTCCGTTCTGGTCACTGAt   >  1:2119299/1‑150 (MQ=255)
                                                                                                                                               tggatggaGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGTCCGTTCTGGTCACTGAtt  <  1:1171366/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GAGCGATAAGCGCTACCGCAGTGACCCGGCATACCGTGCCCGAGTCCAGGCCAAGCTCCTCCGTTCCAACGTCTTCTGACGCCTGGCCGCAATAAGAAAAACAACAATGAGCGGCTTCACCGACGAACACTTCTCCATCAAGTGGATGGAGGAAGTGGATGGCGTGGACTACTGGGCATCTTCGCCTTTCCGCTACTGGATCGGTCCCGAAGGGACTGGCCGGTACGTGGATGTCCCCGCTGGCTTCATCACTGATGGAGCTTCTGTACCCCGTCCGTTCTGGTCACTGATT  >  NC_002947/2602981‑2603272

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: