Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 62,789 A→G 5.2% L159P (CTC→CCC)  rapA ← RNA polymerase‑binding ATPase and RNAP recycling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091362,7890AG 5.2% 236.1 / 2.2 76L159P (CTC→CCC) rapARNA polymerase‑binding ATPase and RNAP recycling factor
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base G (2/2);  total (37/39)
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.87e-01

CGGGACGATAATTAGCACACGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTT  >  NC_000913/62641‑62936
                                                                                                                                                    |                                                                                                                                                   
cGGGACGATAATTAGCACACGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCACGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCATCAGGACGCGCGGCGCGTGGCGGCGAACAACATCATGAGCGATGTTGAg                                                                                                                                                    >  2:577682/1‑150 (MQ=255)
  ggACGATAATTAGCACACGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCt                                                                                                                                                  <  1:1540642/150‑1 (MQ=255)
   gACGATAATTAGCACACGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTg                                                                                                                                                 <  2:1252092/150‑1 (MQ=255)
     cGATAATTAGCACACGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGAt                                                                                                                                               >  1:783103/1‑150 (MQ=255)
      gATAATTAGCACACGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATg                                                                                                                                              <  1:707844/150‑1 (MQ=255)
          aTTAGCACACGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCggg                                                                                                                                          <  1:270892/150‑1 (MQ=255)
             aGCACACGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATc                                                                                                                                       <  1:1365800/150‑1 (MQ=255)
              gCACACGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCa                                                                                                                                      >  2:105077/1‑150 (MQ=255)
               cacaCGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCTGGATCAg                                                                                                                                     <  2:1407907/150‑1 (MQ=255)
                acacGTTAAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATAAgg                                                                                                                                    >  2:789366/1‑150 (MQ=255)
                  acGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCt                                                                                                                                  <  1:1561312/150‑1 (MQ=255)
                   cGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTg                                                                                                                                 >  1:1563431/1‑150 (MQ=255)
                   cGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATCCCGCATCAgctca                                                                                                                                 <  2:780882/150‑5 (MQ=255)
                     ttCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGGCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGt                                                                                                                               <  1:1032512/150‑1 (MQ=255)
                         gcagcGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGc                                                                                                                           >  1:1333164/1‑150 (MQ=255)
                         gcagcGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGc                                                                                                                           >  2:1246537/1‑150 (MQ=255)
                          cagcGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCt                                                                                                                          <  1:883143/150‑1 (MQ=255)
                                cAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACcgc                                                                                                                    >  1:1620900/1‑150 (MQ=255)
                                 agagagCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACcgcg                                                                                                                   <  2:1426448/150‑1 (MQ=255)
                                  tagagCAGTTGCTGATGCAGGATCATCCCTTATTAAATGGTTTTCCCTAAAACAACTTCGTCAGACAGCAGTACGCGCGGCGCGTGGCGGCGACCAACATCATGAGAGATGTTGAGCTGATGCTTTATCATGCTGTTACTCTGACcgctc                                                                                                                  <  2:347205/149‑3 (MQ=255)
                                    gagCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCTGCAGGACGCGCTGCGCGTGGCGGCGACGAACATCATGAGCGATGTTGAGGTGATGGGGGATGAGTCTTGTACTCTTACCGGGCAg                                                                                                                <  2:1253487/150‑1 (MQ=255)
                                       cAGGTGCTGATGCAGGAGCATCCCGGCTGCAAGGTTTTTCCCAAAACCCCCTTCGGCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGcc                                                                                                             <  2:112713/150‑1 (MQ=255)
                                        aGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCggg                                                                                                                                          <  2:138740/120‑1 (MQ=255)
                                        aGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCggg                                                                                                                                          >  1:138740/1‑120 (MQ=255)
                                         gTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGc                                                                                                           <  1:140435/150‑1 (MQ=255)
                                         gTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGc                                                                                                           >  1:259808/1‑150 (MQ=255)
                                          ttGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCt                                                                                                          <  2:1133567/150‑1 (MQ=255)
                                             cTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTa                                                                                                       >  2:928216/1‑150 (MQ=255)
                                                aTGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACgg                                                                                                    <  2:1563431/150‑1 (MQ=255)
                                                 tGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGc                                                                                                   >  2:1627656/1‑150 (MQ=255)
                                                 tGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGc                                                                                                   <  1:1485566/150‑1 (MQ=255)
                                                 tGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGc                                                                                                   >  1:443246/1‑150 (MQ=255)
                                                     ggATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATAc                                                                                               >  2:842123/1‑150 (MQ=255)
                                                     ggATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATAc                                                                                               >  2:951767/1‑150 (MQ=255)
                                                       atcatcCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCTATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACgg                                                                                             <  2:1369676/150‑1 (MQ=255)
                                                         catcCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGaa                                                                                           >  2:760492/1‑150 (MQ=255)
                                                          atcCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAAc                                                                                          <  2:1343562/150‑1 (MQ=255)
                                                           tcCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACt                                                                                         <  1:52648/150‑1 (MQ=255)
                                                                    cAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATGAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGGTgg                                                                                <  1:626247/150‑1 (MQ=255)
                                                                    cAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGACCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTgg                                                                                <  1:238819/150‑1 (MQ=255)
                                                                     aaTGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGa                                                                               >  2:1341766/1‑150 (MQ=255)
                                                                      aTGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGaa                                                                              <  1:5421/150‑1 (MQ=255)
                                                                           tttCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATAttt                                                                         <  2:1620900/150‑1 (MQ=255)
                                                                              cccTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTAcg                                                                      >  2:1607188/1‑150 (MQ=255)
                                                                              cccTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTAcg                                                                      >  2:1389785/1‑150 (MQ=255)
                                                                                    aCCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTAcgcgcgcg                                                                <  1:1213085/150‑1 (MQ=255)
                                                                                      ccACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGAt                                                              <  1:1607188/150‑1 (MQ=255)
                                                                                           tCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGc                                                         >  1:1228071/1‑150 (MQ=255)
                                                                                                     gcagGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGATCTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGc                                               >  2:77656/1‑150 (MQ=255)
                                                                                                     gcagGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGc                                               <  1:1244159/150‑1 (MQ=255)
                                                                                                        ggACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGt                                            >  1:1084309/1‑150 (MQ=255)
                                                                                                        ggACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGt                                            >  2:1207054/1‑150 (MQ=255)
                                                                                                        ggACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGt                                            >  1:180065/1‑150 (MQ=255)
                                                                                                         gACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTc                                           <  1:1530026/150‑1 (MQ=255)
                                                                                                           cgcgcgGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCa                                         >  1:1149125/1‑150 (MQ=255)
                                                                                                             cgcggcgcgTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATa                                       >  1:779760/1‑150 (MQ=255)
                                                                                                                  cgcgTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTc                                  <  1:842123/150‑1 (MQ=255)
                                                                                                                      tGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATc                              <  1:366902/150‑1 (MQ=255)
                                                                                                                      tGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATc                              <  2:779760/150‑1 (MQ=255)
                                                                                                                      tGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATc                              <  1:711923/150‑1 (MQ=255)
                                                                                                                            cgACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCg                        >  2:243491/1‑150 (MQ=255)
                                                                                                                            agcccAACATCATGAGCGATGTTGGGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCg                        >  1:1392246/4‑150 (MQ=255)
                                                                                                                            agcccAACATCATGAGCGATGTTGGGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCg                        >  1:1386530/4‑150 (MQ=255)
                                                                                                                              aCCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGc                      <  1:951767/150‑1 (MQ=255)
                                                                                                                               ccAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCa                     >  2:115245/1‑150 (MQ=255)
                                                                                                                                 aaCATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCaaa                   >  1:1574938/1‑150 (MQ=255)
                                                                                                                                 aaCATCATGAGCGATCTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATATGGAACTGTTAGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCACGGCaaa                   >  2:472293/1‑150 (MQ=255)
                                                                                                                                  aCATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAAc                  <  2:259808/150‑1 (MQ=255)
                                                                                                                                  aCATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAAc                  >  2:295134/1‑150 (MQ=255)
                                                                                                                                    atcatGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAg                >  1:754261/1‑150 (MQ=255)
                                                                                                                                        tGAGCGATGTTGGGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACgg            <  2:1386530/150‑1 (MQ=255)
                                                                                                                                        tGAGCGATGTTGGGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACgg            <  2:1392246/150‑1 (MQ=255)
                                                                                                                                           gCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTcc         >  2:1170432/1‑150 (MQ=255)
                                                                                                                                              aTGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGc      >  2:99218/1‑150 (MQ=255)
                                                                                                                                                 ttGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGt   <  1:1627656/150‑1 (MQ=255)
                                                                                                                                                  tGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGtt  >  1:997914/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGGACGATAATTAGCACACGTTCAGCAGCGCCAGAGAGCAGTTGCTGATGCAGGATCATCCCGGCTTCAATGGTTTTCCCTAAACCCACTTCGTCAGCCAGCAGGACGCGCGGCGCGTGGCGGCGACCAACATCATGAGCGATGTTGAGCTGATGCGGGATCAGGCTGGTACGCTGACCGCGCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTT  >  NC_000913/62641‑62936

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