Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,478,205 C→A 5.6% S284S (TCC→TCA PP_1291 → PhoH family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,478,2050CA 5.6% 256.7 / 2.5 80S284S (TCC→TCAPP_1291PhoH family protein
Reads supporting (aligned to +/- strand):  ref base C (42/33);  new base A (3/2);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTGCTGCTACTCGACCTGCACCAGGAGCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGAAGCAGAGAAGATGGAGCCG  >  NC_002947/1478056‑1478352
                                                                                                                                                     |                                                                                                                                                   
gaACTGCTGCTACTCGACCTGCACCAGGAGCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTcc                                                                                                                                                     >  1:1738345/1‑150 (MQ=255)
   ctgctgCTACTCGACCTGCACCAGGAGCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGt                                                                                                                                                  <  1:707322/150‑1 (MQ=255)
    tgctgcTACTCGACCTGCACCAGGAGCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCAGGATCTGCGCCAGGTGAACGTGACCGGGGCCGCCGGCTCCGCTa                                                                                                                                                 >  2:496335/1‑150 (MQ=255)
       tgctACTCGACCTGCACCAGGAGCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTaaaa                                                                                                                                              <  1:1352262/150‑1 (MQ=255)
         ctACTCGACCTGCACCAGGAGCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAAcc                                                                                                                                            >  2:1793348/1‑150 (MQ=255)
          tACTCGACCTGCACCAGGAGCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCa                                                                                                                                           >  2:1558173/1‑150 (MQ=255)
            cTCGACCTGCACCAGGAGCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATc                                                                                                                                         <  1:933798/150‑1 (MQ=255)
               gACCTGCACCAGGAGCCTTTGCTGCAACAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTg                                                                                                                                      >  2:869832/1‑150 (MQ=255)
                             gCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGc                                                                                                                        >  2:1741499/1‑150 (MQ=255)
                             gCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGc                                                                                                                        <  2:490847/150‑1 (MQ=255)
                              ccTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGcc                                                                                                                       >  2:1104009/1‑150 (MQ=255)
                              ccTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCAGGTAAAACCATCCTGGCACTGGCTGCGGcc                                                                                                                       >  2:56218/1‑150 (MQ=255)
                                tttGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCAt                                                                                                                     >  2:1734843/1‑150 (MQ=255)
                                   gctgcACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGa                                                                                                                  <  1:1171896/150‑1 (MQ=255)
                                     tgcACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTaaaa                                                                                                                                              <  2:56033/120‑1 (MQ=255)
                                     tgcACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTaaaa                                                                                                                                              >  1:56033/1‑120 (MQ=255)
                                     tgcACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTAAGTGAAAACCATCCTGGCACTGGCTGCGGCCATCGAGc                                                                                                                <  1:627607/150‑1 (MQ=255)
                                     tgcACCAGAAAGCCTGGGGCCTGAAACCACGCGACATCCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCCGGCGAACCTGACCGGCGCCGGCCGCTCCGGTAAAACCATCCCGGCACTGGCAGCCGCCATCGAGc                                                                                                                >  2:446787/1‑150 (MQ=255)
                                        tccAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGAGACTTCTTCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCaaa                                                                                                             <  1:1788991/149‑1 (MQ=255)
                                        aCCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCaaa                                                                                                             >  1:301071/1‑150 (MQ=255)
                                          cAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCAGCGAGCAAAca                                                                                                           >  2:1788991/1‑149 (MQ=255)
                                              aaGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCAGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATgg                                                                                                       <  1:56218/150‑1 (MQ=255)
                                                   tGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGc                                                                                                  <  1:1793348/150‑1 (MQ=255)
                                                              aCCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAAc                                                                                                                                             >  1:1343912/1‑96 (MQ=255)
                                                              aCCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAAc                                                                                                                                             <  2:1343912/96‑1 (MQ=255)
                                                                cACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATc                                                                                                                    >  2:379622/1‑119 (MQ=255)
                                                                cACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATc                                                                                                                    <  1:379622/119‑1 (MQ=255)
                                                                cACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTAc                                                                                         <  1:677749/146‑1 (MQ=255)
                                                                cACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTAc                                                                                         >  2:677749/1‑146 (MQ=255)
                                                                cACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGgc                                                                                     <  1:1734843/150‑1 (MQ=255)
                                                                 aCGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGgcg                                                                                    <  1:1438154/150‑1 (MQ=255)
                                                                  cgcgACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGgcgc                                                                                   <  2:301071/150‑1 (MQ=255)
                                                                     gACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCTCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCTTTGAGCAAACCATGGTCAGCCAGCGCTACCGGCGcatc                                                                                >  1:1645418/1‑150 (MQ=255)
                                                                        aTTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCatcatc                                                                             <  2:285712/150‑1 (MQ=255)
                                                                            atcaAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCa                                                                         >  1:1661311/1‑150 (MQ=255)
                                                                             tcaAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCGGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCAc                                                                        >  1:1036260/1‑150 (MQ=255)
                                                                             tcaAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCAc                                                                        >  2:673934/1‑150 (MQ=255)
                                                                               aaaGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAgc                                                                                              >  2:1134422/1‑126 (MQ=255)
                                                                               aaaGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAgc                                                                                              <  1:1134422/126‑1 (MQ=255)
                                                                                     tGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCg                                                                >  2:1130116/1‑150 (MQ=255)
                                                                                       gcgcTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTg                                                              <  1:673934/150‑1 (MQ=255)
                                                                                         gcTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCatcatc                                                                             <  2:1166209/133‑1 (MQ=255)
                                                                                         gcTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCatcatc                                                                             >  1:1166209/1‑133 (MQ=255)
                                                                                           tGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAgg                                                          >  1:1016128/1‑150 (MQ=255)
                                                                                               cGCCCTGCTGGCCCCGGAGATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGTCTGCGGCCATCGAGCACACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCt                                                      <  2:1004301/150‑1 (MQ=255)
                                                                                                gCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACACGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTg                                                     <  1:496335/150‑1 (MQ=255)
                                                                                                     gctgGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCa                                                <  1:1741499/150‑1 (MQ=255)
                                                                                                       tgGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCa                                                                                                    >  2:1390773/1‑96 (MQ=255)
                                                                                                       tgGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCa                                                                                                    <  1:1390773/96‑1 (MQ=255)
                                                                                                       tgGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAgg                                              >  1:256643/1‑150 (MQ=255)
                                                                                                         gACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGgag                                            >  1:1568114/1‑150 (MQ=255)
                                                                                                          aCCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGgaga                                           >  2:1568860/1‑150 (MQ=255)
                                                                                                          aCCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGgaga                                           >  1:1636846/1‑150 (MQ=255)
                                                                                                            ccGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTg                                                     >  2:192043/1‑138 (MQ=255)
                                                                                                            ccGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTg                                                     <  1:192043/138‑1 (MQ=255)
                                                                                                            ccGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGAtt                                         <  1:869832/150‑1 (MQ=255)
                                                                                                            ccGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGAtt                                         >  1:607751/1‑150 (MQ=255)
                                                                                                              ggATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGGTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTgg                                       >  2:1032865/1‑150 (MQ=255)
                                                                                                              ggATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTgg                                       <  2:607751/150‑1 (MQ=255)
                                                                                                              ggATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTgg                                       >  1:735524/1‑150 (MQ=255)
                                                                                                               gATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGc                                      <  1:1558173/150‑1 (MQ=255)
                                                                                                                  atCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTc                                   <  1:1194161/150‑1 (MQ=255)
                                                                                                                      aCCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGc                               >  2:1557975/1‑150 (MQ=255)
                                                                                                                      aCCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGc                               >  2:1921751/1‑150 (MQ=255)
                                                                                                                      aCCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGc                               >  1:1834608/1‑150 (MQ=255)
                                                                                                                       ccTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGcc                              >  2:1887753/1‑150 (MQ=255)
                                                                                                                         tGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCgg                            <  1:1032865/150‑1 (MQ=255)
                                                                                                                         tGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCGTCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCgg                            <  1:1860224/150‑1 (MQ=255)
                                                                                                                           gTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCACACCATGGTCAGCAAGCGCTACCGGCGCACCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGt                          >  2:214527/1‑150 (MQ=255)
                                                                                                                               aCCTTACCGGCGCGGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCg                      <  2:1645418/150‑1 (MQ=255)
                                                                                                                                ttttaCCGGAGCCGGCGGCGCCGGTAAAACCATCCTGGCACTGGCTGCGGCAATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGa                     >  1:1102932/5‑150 (MQ=255)
                                                                                                                                cggagccGGCGGCGCCGGTCAGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGaa                    >  2:1639101/6‑150 (MQ=255)
                                                                                                                                 cTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGaa                    >  2:150528/1‑150 (MQ=255)
                                                                                                                                 cTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGaa                    <  1:1568860/150‑1 (MQ=255)
                                                                                                                                  tGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGAAg                   >  1:195736/1‑150 (MQ=255)
                                                                                                                                     ccGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGAAGCag                <  2:1427115/150‑1 (MQ=255)
                                                                                                                                      cGGAGCCGGCGGCGCAGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGACATTTGCTTCCTGCCGGGTACCGAAGCaga               >  1:1735316/1‑150 (MQ=255)
                                                                                                                                       ggCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGAAGCagag              >  2:1742845/1‑150 (MQ=255)
                                                                                                                                          gccgccGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGAAGCAGagaag           <  1:1219052/150‑1 (MQ=255)
                                                                                                                                              ccagCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGAAGCAGAGAAGATgg       >  2:129995/4‑150 (MQ=255)
                                                                                                                                                   tCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGAGTTACATCAGGAGATTGGCTTCCTGCCGGGTACCGAAGCAGAGAAGATGGAGCCg  >  2:51279/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAACTGCTGCTACTCGACCTGCACCAGGAGCCTTTGCTGCACCAGGAAGCCTGGGGCCTGAAACCACGCGACATTCATCAAAGCCTGGCGCTGTTCGCCCTGCTGGACCCGGATATCCACCTGGTGAACCTGACCGGCGCCGCCGGCTCCGGTAAAACCATCCTGGCACTGGCTGCGGCCATCGAGCAAACCATGGTCAGCAAGCGCTACCGGCGCATCATCGCCACCCGCAGCGTGCAGGGGCTGGATCAGGAGATTGGCTTCCTGCCGGGTACCGAAGCAGAGAAGATGGAGCCG  >  NC_002947/1478056‑1478352

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: