Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,730,641 C→T 5.7% L142L (CTG→TTG)  PP_3297 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,730,6410CT 5.7% 252.6 / 3.2 87L142L (CTG→TTG) PP_3297hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/44);  new base T (2/3);  total (39/48)
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.67e-01

ACACTGTTTCACGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCAGTAATCGAGGCCCGTGGCCTGG  >  NC_002947/3730497‑3730781
                                                                                                                                                 |                                                                                                                                            
acacTGTTTCACGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATg                                                                                                                                         >  1:251762/1‑150 (MQ=255)
 cacTGTTTCACGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATgg                                                                                                                                        >  1:320715/1‑150 (MQ=255)
 cacTGTTTCACGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATgg                                                                                                                                        <  2:958412/150‑1 (MQ=255)
  acTGTTTCACGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGATCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATggg                                                                                                                                       >  2:2296554/1‑150 (MQ=255)
  acTGTTTCACGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATggg                                                                                                                                       >  2:872042/1‑150 (MQ=255)
  acTGTTTCACGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATggg                                                                                                                                       >  2:688833/1‑150 (MQ=255)
   cTGTTTCACGCGATGCCTGCAGCTATTCGTGGGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATgggg                                                                                                                                      >  1:1244922/1‑150 (MQ=255)
    tGTTTCACGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGt                                                                                                                                     <  1:1922186/150‑1 (MQ=255)
    tGTTTCACGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGt                                                                                                                                     >  1:1797969/1‑150 (MQ=255)
      tttCACGCGATGCCTGCAGCTGTTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTtg                                                                                                                                   >  2:812868/1‑150 (MQ=255)
          aCGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTtgctgc                                                                                                                               >  2:967986/1‑150 (MQ=255)
          aCGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTtgctgc                                                                                                                               <  2:251762/150‑1 (MQ=255)
            gcgATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCgg                                                                                                                             >  1:1798328/1‑150 (MQ=255)
            gcgATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCgg                                                                                                                             >  2:2263924/1‑150 (MQ=255)
             cgATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGt                                                                                                                            <  1:165868/150‑1 (MQ=255)
                tGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGa                                                                                                                         >  1:276870/1‑150 (MQ=255)
                     gCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTg                                                                                                                    <  1:2263924/150‑1 (MQ=255)
                     gCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTg                                                                                                                    <  1:1934732/150‑1 (MQ=255)
                              cgcgAGAGCCTGGCCCTGAACACCCAGGAACACTCATTCGAGCGGTATCGCGATGTTTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGtg                                                                                                           <  1:368554/150‑1 (MQ=255)
                                 gagaGCCTGGCCCTGAACAGCCAGGAACACTCATTCGATCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGtggtt                                                                                                        <  1:2296554/150‑1 (MQ=255)
                                 gagaGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGtggtt                                                                                                        <  1:2271281/150‑1 (MQ=255)
                                  agagCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGtggttg                                                                                                       >  1:1698671/1‑150 (MQ=255)
                                  agagCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGtggttg                                                                                                       <  1:95626/150‑1 (MQ=255)
                                   gagCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGtggttgg                                                                                                      >  1:1800624/1‑150 (MQ=255)
                                       cTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTtac                                                                                                  >  1:1431845/1‑150 (MQ=255)
                                       cTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTtac                                                                                                  <  1:2450363/150‑1 (MQ=255)
                                             cTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAg                                                                                            <  2:607726/150‑1 (MQ=255)
                                             cTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAg                                                                                            <  1:812868/150‑1 (MQ=255)
                                              tGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAgg                                                                                           >  2:2204886/1‑150 (MQ=255)
                                                aaCAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGaa                                                                                         <  1:967986/150‑1 (MQ=255)
                                                  cAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATggg                                                                                                                                       >  1:871133/1‑102 (MQ=255)
                                                  cAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATggg                                                                                                                                       <  2:871133/102‑1 (MQ=255)
                                                   aGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCACCGTGGTGCCCTACGACAATGAACGa                                                                                      <  2:1797969/150‑1 (MQ=255)
                                                    gCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAg                                                                                     >  2:517675/1‑150 (MQ=255)
                                                         gAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCg                                                                                >  2:409974/1‑150 (MQ=255)
                                                              cTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCa                                                                              <  1:2342228/147‑1 (MQ=255)
                                                              cTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCa                                                                              >  2:2342228/1‑147 (MQ=255)
                                                              cTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGa                                                                           >  1:729259/1‑150 (MQ=255)
                                                               tCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGTTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTta                                                                                                   >  1:1978499/1‑125 (MQ=255)
                                                               tCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGTTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTta                                                                                                   >  1:1978846/1‑125 (MQ=255)
                                                               tCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGTTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTta                                                                                                   <  2:1978846/125‑1 (MQ=255)
                                                               tCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGTTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTta                                                                                                   <  2:1978499/125‑1 (MQ=255)
                                                               tCATTCGAGCGGTATCGCGATGTGAACCGTGGCAGCCGGTATCAGTACGAGTCGGCCAGCTGGAACTTTGCCGATCCAGTGCTGATGGGGTTGCTGCAGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGCACGAGGCTCGCAGGAt                                                                          >  1:1881191/1‑150 (MQ=255)
                                                                 aTTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATcc                                                                        >  1:1658368/1‑150 (MQ=255)
                                                                     gAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACCTCCTCTGCAGACACATATCTGCCGATCCAGTGCTGATGGGGTTGCTGCAGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTc                                                                    <  2:1881191/150‑1 (MQ=255)
                                                                        cGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTg                                                                 <  2:1658368/150‑1 (MQ=255)
                                                                            aTCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGTTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCggt                                                             <  1:688833/150‑1 (MQ=255)
                                                                            aTCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCtgt                                                             <  1:872042/150‑1 (MQ=255)
                                                                             tCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCtgt                                                             >  1:1090833/1‑149 (MQ=255)
                                                                             tCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCtgt                                                             <  2:1090833/149‑1 (MQ=255)
                                                                             gcgcgATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTtacta                                                                                                <  1:2188255/113‑1 (MQ=255)
                                                                             gcgcgATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTtacta                                                                                                >  2:2188255/2‑114 (MQ=255)
                                                                                   tgtgTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACAt                                                      <  1:517675/150‑1 (MQ=255)
                                                                                       ttCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAGCGAGGCTCGCAGTATCCGTTCGTGCTGTGTTACATCGCt                                                  <  1:2204886/150‑1 (MQ=255)
                                                                                         cACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGGGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGa                                                >  1:1812390/1‑150 (MQ=255)
                                                                                           ccAACATCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAgg                                              <  1:1121031/150‑1 (MQ=255)
                                                                                           ccAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAgg                                              >  1:1053673/1‑150 (MQ=255)
                                                                                           ccAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAgg                                              <  1:709741/150‑1 (MQ=255)
                                                                                           ccAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAgg                                              >  2:1125367/1‑150 (MQ=255)
                                                                                            cAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGtggttgg                                                                                                      <  1:1116163/93‑1 (MQ=255)
                                                                                            cAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGtggttgg                                                                                                      >  2:1116163/1‑93 (MQ=255)
                                                                                                aGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCaa                                         <  1:2437582/150‑1 (MQ=255)
                                                                                                   cGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCAAACTCAGTGAGCACTGCATTGTTGGGTTGTTTGTTGGCTAATGGGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCt                                      <  1:1116895/150‑1 (MQ=255)
                                                                                                    ggTATCAGTACGAGTAGGACCGCTGGACAATTGGCGACCTAGTATAAATGTGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCtt                                     <  2:1204705/150‑1 (MQ=255)
                                                                                                     gTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGt                                                                                                                      <  1:2244803/68‑1 (MQ=255)
                                                                                                     gTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGt                                                                                                                      >  2:2244803/1‑68 (MQ=255)
                                                                                                       aTCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCtg                                                              >  2:229176/1‑122 (MQ=255)
                                                                                                       aTCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCtg                                                              <  1:229176/122‑1 (MQ=255)
                                                                                                       aTCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGAt                                  <  2:1431845/150‑1 (MQ=255)
                                                                                                          aGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGAt                                  >  2:1681235/1‑147 (MQ=255)
                                                                                                          aGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGAt                                  <  1:1681235/147‑1 (MQ=255)
                                                                                                           gTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTg                                                                                                                           <  1:84385/57‑1 (MQ=255)
                                                                                                           gTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTg                                                                                                                           >  2:84385/1‑57 (MQ=255)
                                                                                                                gAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGcc                          <  2:729259/150‑1 (MQ=255)
                                                                                                                  gTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTCGACGTTGGCCAAGGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTACTCTGCGGCCCa                        >  2:1328094/1‑150 (MQ=255)
                                                                                                                  gTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGATGCGGTGCACGTTGTCCAACGTGTTTGTTTACTCACAGGAACGAGGCTCAAAGCTTCTTGTCTTTCGTTTTGACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCa                        <  1:1863283/150‑1 (MQ=255)
                                                                                                                      gCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCAGTaa                    >  2:792722/1‑150 (MQ=255)
                                                                                                                       ccAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCAGTAAt                   <  1:1950529/150‑1 (MQ=255)
                                                                                                                       ccAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGAAGTAAGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCAGTAAt                   <  2:1244922/150‑1 (MQ=255)
                                                                                                                            tGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCAGTAATCGAgg              <  1:409974/150‑1 (MQ=255)
                                                                                                                            tGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCAGTAATCGAgg              <  1:2164963/150‑1 (MQ=255)
                                                                                                                                 aTTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGCCCTGCGGCCAAGTTATCGTGTcccag         >  2:510570/1‑148 (MQ=255)
                                                                                                                                    tGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCAGTAATCGAGGCCCGTGGc      <  2:276870/150‑1 (MQ=255)
                                                                                                                                      ccGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCAGTAATCGAGGCCCGTGGCCt    <  2:1053673/150‑1 (MQ=255)
                                                                                                                                       cGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCgg                            >  1:2051932/1‑125 (MQ=255)
                                                                                                                                       cGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCgg                            <  2:2051932/125‑1 (MQ=255)
                                                                                                                                        gATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCAGTAATCGAGGCCCGTGGCCTgg  <  1:1177519/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
ACACTGTTTCACGCGATGCCTGCAGCTATTCGCGAGAGCCTGGCCCTGAACAGCCAGGAACACTCATTCGAGCGGTATCGCGATGTGTTCACCAACAGCCGGTATCAGTACGAGTCGGCCAGCTGGAAATTTGCCGATCCAGTGCTGATGGGGTTGCTGCGGTGGACGTTGGCCAACGTGGTTGGTTACTACAAGGAACGAGGCTCGCAGGATCCGTTCGTGCTGTGTTACATCGCTGAGGTCCAAGCTTGATCTGCGGCCCAGTAATCGAGGCCCGTGGCCTGG  >  NC_002947/3730497‑3730781

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