Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,187,214 A→C D622A (GAC→GCC)  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,2140AC96.3% 374.2 / ‑3.0 107D622A (GAC→GCC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base C (50/53);  total (52/55)
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

Only 100 of 110 total aligned reads displayed.

TTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTC  >  NC_000913/4187070‑4187358
                                                                                                                                                |                                                                                                                                                
ttGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGccc                                                                                                                                                <  6:375802/147‑1 (MQ=255)
  ggCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGa                                                                                                                                             <  1:49011/148‑1 (MQ=255)
        gCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTa                                                                                                                                      >  4:195116/1‑149 (MQ=255)
          cATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTaca                                                                                                                                    >  3:33400/1‑149 (MQ=255)
           aTTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTacac                                                                                                                                   >  6:357356/1‑149 (MQ=255)
                gtgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCt                                                                                                                              >  5:109400/1‑149 (MQ=255)
                 tgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGc                                                                                                                               <  1:194188/147‑1 (MQ=255)
                   gatgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTc                                                                                                                            <  4:387524/148‑1 (MQ=255)
                   gatgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGGCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTc                                                                                                                            <  2:63858/148‑1 (MQ=255)
                    atgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGc                                                                                                                          <  5:380326/149‑1 (MQ=255)
                    atgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGc                                                                                                                          >  8:332574/1‑149 (MQ=255)
                    atgagTGTCCCGAAAGGTCTGCCTTACTACAGCGTCAACCAGGCGCTGCGTAAAAAAGCAACCTCCAAAATGCGGAACACCTGCTACCGCATTCTCTGGCTGAACCCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGc                                                                                                                          <  5:271539/149‑1 (MQ=255)
                     tgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGcc                                                                                                                         >  3:69552/1‑149 (MQ=255)
                      gatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGcc                                                                                                                         >  6:406520/1‑148 (MQ=255)
                       atTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTa                                                                                                                       <  5:406520/149‑1 (MQ=255)
                        ttGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCAt                                                                                                                                         >  1:216219/1‑130 (MQ=255)
                        ttGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCAt                                                                                                                                         <  2:216219/130‑1 (MQ=255)
                           tACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGCAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATg                                                                                                                     >  1:198128/1‑147 (MQ=255)
                             ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATg                                                                                                                                        <  8:327792/126‑1 (MQ=255)
                             ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATg                                                                                                                                        >  7:327792/1‑126 (MQ=255)
                             ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAgc                                                                                                                 <  6:59343/149‑1 (MQ=255)
                              cGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAgcg                                                                                                                >  1:18115/1‑149 (MQ=255)
                              cGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAgcg                                                                                                                >  3:303714/1‑149 (MQ=255)
                              cGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAgcg                                                                                                                >  8:394544/1‑149 (MQ=255)
                                aaaGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGcgcg                                                                                                              >  4:31344/1‑149 (MQ=255)
                                     tCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTg                                                                                                         <  3:31344/149‑1 (MQ=255)
                                     tCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTg                                                                                                         >  2:175273/1‑149 (MQ=255)
                                     tCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTg                                                                                                         >  6:36525/1‑149 (MQ=255)
                                       tGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTg                                                                                                         >  4:91328/1‑147 (MQ=255)
                                       tGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTg                                                                                                         <  3:91328/147‑1 (MQ=255)
                                          cTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCa                                                                                                    <  7:100757/149‑1 (MQ=255)
                                            tACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTg                                                                                                         <  1:119129/142‑1 (MQ=255)
                                               tCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTg                                                                                                         >  7:252823/1‑139 (MQ=255)
                                               tCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGt                                                                                               >  3:181486/1‑149 (MQ=255)
                                               tCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGATGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGt                                                                                               >  3:18/1‑149 (MQ=255)
                                                 cATCGTCAAGCAGGCGCTGGGTTAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATTCAGCGCGTTCTGGTGCATCTGtt                                                                                              <  3:174246/148‑1 (MQ=255)
                                                 cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGccc                                                                                                                                                >  8:22858/1‑98 (MQ=255)
                                                 cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                             <  6:163550/149‑1 (MQ=255)
                                                      tCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATc                                                                                        >  4:339370/1‑149 (MQ=255)
                                                        aaCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCg                                                                                       <  2:26426/148‑1 (MQ=255)
                                                        aaCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCg                                                                                       <  6:196332/148‑1 (MQ=255)
                                                         aCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCAGAACCCCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                                     <  8:321847/149‑1 (MQ=255)
                                                          ccAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                                     >  5:176098/1‑148 (MQ=255)
                                                          ccAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                                     >  1:293662/1‑148 (MQ=255)
                                                          ccAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATg                                                                                    <  7:63498/149‑1 (MQ=255)
                                                           cAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATGCTCGGTCTGAACCCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCGTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGa                                                                                   <  4:189230/149‑1 (MQ=255)
                                                              gcgcTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACAt                                                                                >  8:313570/1‑149 (MQ=255)
                                                               cgcTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATg                                                                               <  5:357356/149‑1 (MQ=255)
                                                                 cTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGt                                                                             <  5:36525/149‑1 (MQ=255)
                                                                  tGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                                     >  7:416535/1‑140 (MQ=255)
                                                                    ggTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCa                                                                           <  3:209773/148‑1 (MQ=255)
                                                                          aaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGga                                                                    <  5:49403/149‑1 (MQ=255)
                                                                            aGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGgag                                                                   >  8:416236/1‑148 (MQ=255)
                                                                            aGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGgag                                                                   >  5:242564/1‑148 (MQ=255)
                                                                             gCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaa                                                                 <  4:69552/149‑1 (MQ=255)
                                                                              cAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGgag                                                                   >  7:406132/1‑146 (MQ=255)
                                                                              cAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaag                                                                <  6:176098/149‑1 (MQ=255)
                                                                                   tcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAACCCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaa                                                                 <  6:245267/143‑1 (MQ=255)
                                                                                   tcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaa                                                                 >  5:245267/1‑143 (MQ=255)
                                                                                   tcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                           >  4:62673/1‑149 (MQ=255)
                                                                                    ccAAAATGCTGAACACCTGCTACCTCAGTCTCGGTCTGACACCGACCGTTATTTTTGCGGACCAGCTAATTTCCACCGGCTGAGCATCTGCCGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCAGCCCGGAGAAGCAacac                                                          >  2:139901/1‑149 (MQ=255)
                                                                                    ccAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGAc                                                                                  >  6:145899/1‑125 (MQ=255)
                                                                                    ccAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAacac                                                          >  2:139898/1‑149 (MQ=255)
                                                                                     cAAAATGCTGACCACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAac                                                            <  6:109400/146‑1 (MQ=255)
                                                                                     cAAAAAGCCGTACACCTTCTACCCGATTCTCGGGCCGAAACCGAACGTTATTTTTGCCGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACg                                                         <  2:209104/149‑1 (MQ=255)
                                                                                       aaaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGaa                                                       <  5:213233/149‑1 (MQ=255)
                                                                                        aaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGaaa                                                      <  4:33400/149‑1 (MQ=255)
                                                                                          tGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatc                                                    <  3:392710/149‑1 (MQ=255)
                                                                                          tGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatc                                                    <  7:332574/149‑1 (MQ=255)
                                                                                                  acCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGa                                            <  7:416236/149‑1 (MQ=255)
                                                                                                    catatcccGCATTATCGTTATCAAACCTCCATTTATTTTTTCGACCAGATAATGTACACCTGCTTAGCCTATTCATCGCGTTCTGGTGCATATGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACTAAATCATCTCCGAgg                                          <  7:322238/142‑1 (MQ=255)
                                                                                                    cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                          >  4:81639/1‑149 (MQ=255)
                                                                                                    cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                          <  4:303714/149‑1 (MQ=255)
                                                                                                    cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                          >  3:42387/1‑149 (MQ=255)
                                                                                                       cTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                          >  5:156528/1‑146 (MQ=255)
                                                                                                        tACCGCATTCTCGGTCTGAAACCGCCCGTTATTTTTGATCCCAGATAATTTAAACCAGAGTCGCCTATTCATCGCGGTAGGGGTCATCGGTTTGTATAGGATAAAAGGGCATACCGGAGAAGAAACAAGAAATACTATACGAGgcc                                        >  8:322238/1‑145 (MQ=255)
                                                                                                          ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaag                                    >  8:363313/1‑149 (MQ=255)
                                                                                                          ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaa                                     <  3:195116/148‑1 (MQ=255)
                                                                                                           cGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGGAGCGATGACATGGGCATCCCGGAGAAGAACCACGAAATCATCTCCGAGgcagaag                                    >  5:312339/1‑148 (MQ=255)
                                                                                                            gCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagca                                  >  3:101625/1‑149 (MQ=255)
                                                                                                             cATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaag                                    >  2:25255/1‑146 (MQ=255)
                                                                                                             cATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaag                                    <  1:25255/146‑1 (MQ=255)
                                                                                                                  tcGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGtt                            >  2:280216/1‑149 (MQ=255)
                                                                                                                      tCTGAACCCGACCGGTAGTTTTGCGGCCCCGATCATGTCCACCGGCTGCGCCTATGCAGCGCGTTCTGGTGCATCTGTGGGTCTCGAGGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTg                        <  6:312339/149‑1 (MQ=255)
                                                                                                                       cTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGa                       <  2:18115/149‑1 (MQ=255)
                                                                                                                        tGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGa                       >  1:185062/1‑148 (MQ=255)
                                                                                                                         gAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaaa                     >  5:94047/1‑149 (MQ=255)
                                                                                                                             ccgaccgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCa                 <  6:373345/149‑1 (MQ=255)
                                                                                                                             ccgaccgTTATTTTTGCGGCCAAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGGTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCa                 <  2:293662/149‑1 (MQ=255)
                                                                                                                              cgaccgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                >  4:88887/1‑149 (MQ=255)
                                                                                                                              cgaccgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                <  4:42387/149‑1 (MQ=255)
                                                                                                                                accgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGa              <  6:94047/149‑1 (MQ=255)
                                                                                                                                 ccgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACg                                                         >  5:132524/1‑105 (MQ=255)
                                                                                                                                 ccgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACg                                                         <  6:132524/105‑1 (MQ=255)
                                                                                                                                  cgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGc            <  3:88887/149‑1 (MQ=255)
                                                                                                                                   gTTATTTTTTCGGCCCAGTTCATGTACACCGGCTTCGCCTATGCAGCGGGTTCTGGTGCAGCTGTTGGTTTCTATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCGGAGGATGAAGCAGAAGTTGGTTAAATTCAGGAGCa           <  3:92820/149‑1 (MQ=255)
                                                                                                                                     tATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGt         >  8:318772/1‑149 (MQ=255)
                                                                                                                                       tttttGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTc       >  7:261670/1‑149 (MQ=255)
                                                                                                                                            gCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgt   <  8:261670/148‑1 (MQ=255)
                                                                                                                                             cGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtc  <  1:280216/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTC  >  NC_000913/4187070‑4187358

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