Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,187,507 A→C N720H (AAC→CAC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,187,5070AC98.0% 150.8 / NA 50N720H (AAC→CAC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (26/22);  total (27/23)
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.92e-01

CCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACC  >  NC_000913/4187215‑4187797
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                   
ccAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAAATGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCGGGTTTCCTTCCACAGCATc                                                                                                                                                                                                                                                                                             >  2:436298/1‑301 (MQ=255)
  aGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCt                                                                                                                                                                                                                                                                                            >  2:459195/1‑300 (MQ=255)
      actgtagcccGGGTTCGCTCATGAAGGGGGTTCGGGGGCATCTGTTGGTATCGATGACATGGGCATCCCGGAGAAGAAAAACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATAGCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACa                                                                                                                                                                                                                                                                                         <  1:350749/291‑1 (MQ=255)
                             cAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGCCAGGAAGAGAAGAAGGTTTCCTTCCACAACATCAACATGATGGCCGCCGTGGGTGCGCGTgg                                                                                                                                                                                                                                                                >  1:153985/1‑301 (MQ=255)
                                        ggTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAACTCATCTCCGAGGCAGAAGCAGAAGTTGGTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGcac                                                                                                                                                                                                                                                     <  1:178779/301‑1 (MQ=255)
                                                tgatGGTACCGTTGACCATGGCCACCCGGAAAAAAAACACGAAATCAGCCCCGAGGAAGAAGCAGAATTTTCGGAAATTCGGGGGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGGTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTc                                                                                                                                                                                                                                               <  1:336043/296‑1 (MQ=255)
                                                tGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGGGAAGCAGGTTTCCTTCCACAGCACCTACAAGATGGCCGACTCAGGTGCGCGGGGTTCTGGGGCACGGAtt                                                                                                                                                                                                                                                >  1:97645/1‑298 (MQ=255)
                                                tGTTGGTATCGATGACATGGTCAGCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAtt                                                                                                                                                                                                                                                <  2:97669/298‑1 (MQ=255)
                                                 gTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTAAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACGCCGGTGCGCGTGGGTCTGCGGCACAGATTcagg                                                                                                                                                                                                                                            >  2:336136/1‑298 (MQ=255)
                                                     ggtaatgatgcactGGTCATACTGAAGAAGAAACACGAAATCGTCTCCGAGGCGGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCAAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGc                                                                                                                                                                                                                                        <  1:69734/288‑1 (MQ=255)
                                                     gaaTCGATGAAATGGTCAACCCGGAGAAGAAACACGAAATCATTTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGGAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt                                                                                                                                                                                                                                       <  1:459071/299‑1 (MQ=255)
                                                             gACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTAACGGTAAGGAAGAAAAGCAGGTTTCCTTACACAGCATCTAAATGATGGCCGACTCCGGTGCGAGTGGTTCTGCGGCACAGATTATGTAGCTTGCTGGtt                                                                                                                                                                                                                               >  2:367769/1‑300 (MQ=255)
                                                                              agaagaAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACCCCGGTGCGATGGGTTCGGCGGCCCAGAGTCGTCAGCTTGCTGGTAGGCGTGGTCTGAgggg                                                                                                                                                                                                               >  1:367906/1‑300 (MQ=255)
                                                                                 agaaACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGTTCAGGAAGAGAAGCAGGTTTACTGCCACAGCATATACATGATGGCCGACTCTGGTTCGCGAGGTTCTGCGGCCCGGCTCGTGCAGCTTGCTGGTAGGCGGGGTTTGATGGCGa                                                                                                                                                                                                             >  1:144574/1‑300 (MQ=255)
                                                                                     acacGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACCCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGGGGGGTCTGTTGGCGAAGCCg                                                                                                                                                                                                        >  1:429185/1‑301 (MQ=255)
                                                                                                tctcCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTTTCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCAAAGCATCTACATGATGGCAGACTCCGGTGACCGTGTTTTGGAGGAACAGGTTCGTCAGCTTGATGGTATCGGTGGTCTAATGGCGAATACGGCTGGCTCcat                                                                                                                                                                                             >  2:133688/1‑301 (MQ=255)
                                                                                                   ccGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCCACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCcctcat                                                                                                                                                                                          >  2:29426/1‑301 (MQ=255)
                                                                                                                   agTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGTTGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCTTCAGCTTGCTGGTATGCGTGGTCTGTTGGCGATGCCGGATGGCTCCATCATCGAACCGCCACTcccc                                                                                                                                                                          >  2:47449/1‑301 (MQ=255)
                                                                                                                      tGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAATTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTTATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAAAAGCCTCGAAATGATGGCCGATATCGGTACGCGTGCTTCCGAGGCACAGATTAGTCAGCTAGCTGGTATGCGTGGTGTGGTGGGGAACGCGGTAGGCCCCTTCACGGAACCGGAAATCACtgcg                                                                                                                                                                       >  2:108934/1‑301 (MQ=255)
                                                                                                                          aaaaaTCAGGAGCAGTTCCAGTCTGGTCTGGAAATTGCGGGCGAACGCTACAAAAAAGTTATCGATATCTGGGAGGCGGCGAACGATCGTGGATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACt                                                                                                                                                                   <  1:367672/300‑1 (MQ=255)
                                                                                                                              gtCAGGAGCAGATCGAGTCTGTGCGGTTACAGGAGGACGAACTATACAAAAAAGGTAAGGATTACTCGGGTGTGGCGAACGATCGGTGATCCAAAGCGCTGATGGATAACCTGCAAACTGCAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCg                                                                                                                                                               <  1:436182/300‑1 (MQ=255)
                                                                                                                                                      gggtAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTa                                                                                                                                       <  2:154027/300‑1 (MQ=255)
                                                                                                                                                              gCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGGTGCGGCGAACGATCGTTTATCCAAAGCGATGATGGATAACCAGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGAGGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtct                                                                                                                               <  2:144610/301‑1 (MQ=255)
                                                                                                                                                                         cTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAAGCGTGGTCTTATGGCGAAGCCGGATGGCCCCATCAACGAAAAGCCAACCACCGCGAAATCCCGTGCAGGTCTGAACGTACCCCAGTCCTTCATCTCCCCCCCCggc                                                                                                                    >  1:425074/1‑300 (MQ=255)
                                                                                                                                                                               caaAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAgg                                                                                                                                                         <  2:395011/258‑1 (MQ=255)
                                                                                                                                                                                  aGTAATCGATAGCGGTGTGGCGGCGAACGAGCGTGTGTCAGAAGCGATGATGGATAACCTGCAAACAGATGCCGTGATTAACCGTGACGGTCAGGAAGAGAAACAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTaaa                                                                                                           <  2:199383/301‑1 (MQ=255)
                                                                                                                                                                                     tATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAgg                                                                                                         <  2:429296/300‑1 (MQ=255)
                                                                                                                                                                                          atatCTGGGGTGGGGCGCGCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGc                                                                                                   <  1:166047/301‑1 (MQ=255)
                                                                                                                                                                                          aaacggcggCTGCGGCGAAAGAGCGTGTAGTCAAGACGATGATGGAGAACCTGCAAACTGAAACCGGGATTAACCGTGACGGTCAGTAAGATAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAGGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGc                                                                                                   <  1:47437/294‑1 (MQ=255)
                                                                                                                                                                                           aaacagggacgggacgAACGTTTGTGTAGCAAAAGCGATTATGGAGTACCTGAAAACTGAAAGCGTTATTTACCGTGACGGGCCGGAAGACAAGCAGGTTTCCTCCGACAGCGTCTTCATGCTGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGCCTCCATCATCGAAACGCCAATCACCGCGATCTTCCGTGAAGGTCTGAACGTACTCCAGTACTTTATTTCCACCCACGGTGCTCGTAAAGGTCTGGCg                                                                                                  <  2:81758/287‑1 (MQ=255)
                                                                                                                                                                                            atCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCgg                                                                                                 >  2:454207/1‑301 (MQ=255)
                                                                                                                                                                                                             cgagcgttgacccAAAGCGATGATGGATAACCTGCAAACTGAGACCGTGATGAACCGTTACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACGGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACt                                                                                <  1:454086/290‑1 (MQ=255)
                                                                                                                                                                                                             cGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATACCCACCCACGGTGCTCGTACAGGTCTGGCGGATACCGCACTGAAACCt                                                                                >  1:63629/1‑301 (MQ=255)
                                                                                                                                                                                                                    gtttCCAAAGCGATTATGGACAACCTGAGAAAGTAAACCTGTATGACCCGTGACGGTCAGGAATAACAGCAGGTTTCCTTCCAAAGCATCTACATGTTGGCCGACATCGGTAGGCGTGGTTCTTAGCCAACGATTCTACCGCTTGCTGGGATGCGTCGTCTGAGGGAGAAGCCGGATGGCTCCATCATCGAAACTCCCATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAAc                                                                          <  1:133655/297‑1 (MQ=255)
                                                                                                                                                                                                                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTTATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCCCCATCATCGAAACGCCAATCACCGCGCAATTCCGGGAGGGTCTGAACGTACCCCAGCACTTCATCTCCACCCCCGCTGCTCGTCACGGCCTGGCGGATCCCGCACTCACAACAGCGAAccc                                                                        >  1:25988/1‑299 (MQ=255)
                                                                                                                                                                                                                        caaagaccattgtgGATAGCCTGCGATCTGATCCCGTGATTATCCTGGACGGGCAGGAGGACACGCAGCTTTTCTTGCACCGCAGCTAGATGAGTGCCGACGCCTGTGCGCGTGGTTCTGCGGCACAGATCCGTCAGCTTGCTGGGATGCGTGGTCTTATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTCCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCgg                                                                     <  1:49586/289‑1 (MQ=255)
                                                                                                                                                                                                                              cGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCt                                                               <  2:368003/301‑1 (MQ=255)
                                                                                                                                                                                                                                              gCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCCCAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGCACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCAACCCACGGTGCTCGTAACgg                                                                                                         >  2:353223/1‑243 (MQ=255)
                                                                                                                                                                                                                                               cAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGCCTGGTTg                                              >  2:413634/1‑301 (MQ=255)
                                                                                                                                                                                                                                               aaaaCGGAAACCGTGATGAACCGGGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCAACGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTg                                              <  2:425185/300‑1 (MQ=255)
                                                                                                                                                                                                                                                   cTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGGAGCAGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCCCTGAAAACTGCACACTCCGGTTCCCTgctgctccggccggttgccgg                                          >  1:308466/1‑281 (MQ=255)
                                                                                                                                                                                                                                                    tGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTATGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACGTCATCTCCACCCACGGTTCTCGTAAAGGTCTGGCGGATACGGCACTGCAAACTGCGAACGCCGGTTACCTGCCTCGTAGTATGGtagtagtg                                         >  2:136215/1‑294 (MQ=255)
                                                                                                                                                                                                                                                     gacaCCGTGAGTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACAGGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTgg                                        <  2:345122/298‑1 (MQ=255)
                                                                                                                                                                                                                                                      aaaCCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTATAGGTCTGGCGGCTACCGCCTTGAAAACTGAGCACCCCggagagcggtctcgcgtcctgtcaggagagcg                                       >  1:257532/1‑271 (MQ=255)
                                                                                                                                                                                                                                                      aaaCCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCgg                                                                     <  2:257608/271‑1 (MQ=255)
                                                                                                                                                                                                                                                             gATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTTACTCGTCGCCTGGTTGACGTGGCACAGGAc                                >  1:178767/1‑301 (MQ=255)
                                                                                                                                                                                                                                                             gATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGGCTGGCGGATACCGCCCTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGgcccgggcc                                >  1:43129/1‑294 (MQ=255)
                                                                                                                                                                                                                                                                   ccGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGAGACCGCACTGAAAACTGCGCACTCCGGTTACCTGACTCGTCGTCTGGCTGACGTGGCGCAGGACCtggtg                          >  1:371342/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                gAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGGGCAGGTCTGAACGTACTCCAGTACTTCATCCCCACCCACGGTGCCCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGCACTCCGGTTACCTGAACCGTCCTCTGGTTGACGTTGCGACGGCCGTGGGGGTTTCCGAAGACg             >  1:368/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           ggTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTAcc  >  2:312398/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                   
CCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACC  >  NC_000913/4187215‑4187797

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: