Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,187,214 A→T D622V (GAC→GTC)  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,2140AT99.2% 409.9 / NA 124D622V (GAC→GTC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (58/65);  minor base G (1/0);  total (59/65)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

TGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTC  >  NC_000913/4187071‑4187353
                                                                                                                                               |                                                                                                                                           
 ggCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGAt                                                                                                                                       <  1:712774/149‑1 (MQ=255)
 ggCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGAt                                                                                                                                       <  1:712775/149‑1 (MQ=255)
 ggCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGAt                                                                                                                                       <  1:987767/149‑1 (MQ=255)
  gCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATc                                                                                                                                      <  2:858230/149‑1 (MQ=255)
  gCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATc                                                                                                                                      <  1:620584/149‑1 (MQ=255)
   ccGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTTCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCa                                                                                                                                     <  2:1190840/149‑1 (MQ=255)
   ccGTGCCATTCTGTGGATGATTGTACCGAAACGTCTGCCTTAATCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCa                                                                                                                                     <  1:545993/149‑1 (MQ=255)
    cGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCAt                                                                                                                                    <  2:864065/149‑1 (MQ=255)
       gCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTa                                                                                                                                 <  1:793370/149‑1 (MQ=255)
        ccATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCATTATTTTTGCGGTCCAGATCATGTac                                                                                                                                <  2:565362/149‑1 (MQ=255)
         cATTCTGTGGAAGATTGTACCGAAAGGTCTGCCTTACTCCATCGACAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGGAAGATTTGCGGTCCAGAACATGAaca                                                                                                                               >  1:39294/1‑149 (MQ=255)
             cTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCgg                                                                                                                           >  2:100834/1‑149 (MQ=255)
                tgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATTTACACCGGCtt                                                                                                                        >  2:220265/1‑149 (MQ=255)
                tgGATGATCGTTCCGAAAGGTCGGCCTTACTCCATCGTCAACCCGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGACACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCtt                                                                                                                        <  2:115237/149‑1 (MQ=255)
                   atgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGc                                                                                                                     >  1:1110512/1‑149 (MQ=255)
                   atgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGc                                                                                                                     >  1:83291/1‑149 (MQ=255)
                     gatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCt                                                                                                                   <  2:1010067/149‑1 (MQ=255)
                          tACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCa                                                                                                              >  2:40438/1‑149 (MQ=255)
                            ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGCCCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAgc                                                                                                            >  2:109556/1‑149 (MQ=255)
                             cGAAAGGTCTGCCTTACTCCATCGTCAACCATGCGCAGGGTAAAAAAGCAATCTACAAAATGCTGAACACCTGCTACAGCATTCTCGGTCTGAAACCGCCCGTTATTTTTGCGGTCAAGATCATGAACACCGGCTTCTCATATGCAgcg                                                                                                           >  1:627648/1‑149 (MQ=255)
                               aaaGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGcgcg                                                                                                         >  2:1101271/1‑149 (MQ=255)
                               aaaGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTAAACCGGCTTCGCCTATGCAGcgcg                                                                                                         >  2:633833/1‑149 (MQ=255)
                                 aGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGtt                                                                                                       <  1:1131287/149‑1 (MQ=255)
                                 aGGGCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGtt                                                                                                       <  1:697484/149‑1 (MQ=255)
                                  ggTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTc                                                                                                      >  1:1115111/1‑149 (MQ=255)
                                   gTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCt                                                                                                     >  2:346360/1‑149 (MQ=255)
                                      tGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCCGCGCTTTCTGGt                                                                                                  >  2:114495/1‑149 (MQ=255)
                                      tGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGAAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGt                                                                                                  >  1:485570/1‑149 (MQ=255)
                                               ccATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGACACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGtt                                                                                         >  1:114207/1‑149 (MQ=255)
                                               ccATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGtt                                                                                         <  1:100834/149‑1 (MQ=255)
                                                cATCGTCAACCAGGCGCTGGGTAAACAAGCAATCTCCAAAATGCTGAACACCTGCGACCGCTTTCTCGGTATGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                        <  1:618799/149‑1 (MQ=255)
                                                cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                        >  2:57249/1‑149 (MQ=255)
                                                cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                        >  2:620111/1‑149 (MQ=255)
                                                cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                        >  1:580176/1‑149 (MQ=255)
                                                cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACATGCTACCGCATTCTCGGTCTCAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTACTGGTGCATCTGTTg                                                                                        >  2:759183/1‑149 (MQ=255)
                                                    gTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGGCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTAt                                                                                    <  1:1023410/149‑1 (MQ=255)
                                                        aCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                                <  2:622457/149‑1 (MQ=255)
                                                          cAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGa                                                                              <  2:757848/149‑1 (MQ=255)
                                                          cAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGAATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGATTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGa                                                                              <  2:789020/149‑1 (MQ=255)
                                                           aGGCGCTGGGTAAAAAAGCAATCACCAAAATGCTGCACACCTGCTACCTCATTCTCGGTCTGAATCCGACCGTTATTTTTGCGGTCCAGATCAAGTACACCGGCTTCGCCTATGCAGAGCGTTCTGGTGCAGCTGTTCGTATCGATGAc                                                                             >  1:1053827/1‑149 (MQ=255)
                                                               gcTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATgg                                                                         <  2:1110512/149‑1 (MQ=255)
                                                                   ggTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCaa                                                                     >  1:444264/1‑148 (MQ=255)
                                                                    gTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATc                                                                    <  1:620111/149‑1 (MQ=255)
                                                                     tAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAACCCGACCGGTATATTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATcc                                                                   <  2:100255/149‑1 (MQ=255)
                                                                       aaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTTCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGAGTTCTGGTGCTTGTGTTGGTTTCGATGACATGGTCATCCCg                                                                 >  2:510476/1‑149 (MQ=255)
                                                                       aaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCg                                                                 >  2:715814/1‑149 (MQ=255)
                                                                        aaacGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTAGACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCgg                                                                <  2:1016132/149‑1 (MQ=255)
                                                                        aaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCgg                                                                >  2:398028/1‑149 (MQ=255)
                                                                         aaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGga                                                               >  1:558505/1‑149 (MQ=255)
                                                                         aaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGga                                                               <  2:113164/149‑1 (MQ=255)
                                                                         aaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGga                                                               >  1:1251411/1‑149 (MQ=255)
                                                                         aaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGga                                                               >  2:445775/1‑149 (MQ=255)
                                                                            gCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaa                                                            <  1:755420/149‑1 (MQ=255)
                                                                              aaTCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaaga                                                          >  1:153799/1‑149 (MQ=255)
                                                                                 ctcCAAACTGCTGACCACCTGCGACCGCATTCTCGGTCTGAAACCTACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTATGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAac                                                       <  1:54838/149‑1 (MQ=255)
                                                                                  tcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                      >  2:433774/1‑149 (MQ=255)
                                                                                  tcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                      <  1:1201180/149‑1 (MQ=255)
                                                                                  tcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                      >  2:1249203/1‑149 (MQ=255)
                                                                                    cAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACg                                                    <  1:348709/149‑1 (MQ=255)
                                                                                    cAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACg                                                    >  2:1169194/1‑149 (MQ=255)
                                                                                     aaaaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTCTGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGa                                                   <  1:219456/149‑1 (MQ=255)
                                                                                     aaaaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCAGGAGAAGAAACACGa                                                   >  1:1072715/1‑149 (MQ=255)
                                                                                     aaaaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCACGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGa                                                   >  1:144705/1‑149 (MQ=255)
                                                                                       aaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGaaa                                                 <  2:910110/149‑1 (MQ=255)
                                                                                           cTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatcat                                             >  2:578397/1‑149 (MQ=255)
                                                                                            tGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatcatc                                            >  1:849895/1‑149 (MQ=255)
                                                                                            tGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatcatc                                            <  1:1249203/149‑1 (MQ=255)
                                                                                             gAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGATCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtct                                           >  2:32417/1‑149 (MQ=255)
                                                                                                cacCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCg                                        >  1:464678/1‑149 (MQ=255)
                                                                                                 acCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGa                                       <  1:300020/149‑1 (MQ=255)
                                                                                                   cTGCTACCGCATTCTCGGGCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTTCACCTGATTCGCTTATGCATCGCGTTCTGGTTCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                     <  1:188109/149‑1 (MQ=255)
                                                                                                     gCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgca                                   >  1:965292/1‑149 (MQ=255)
                                                                                                       tACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcaga                                 >  1:441610/1‑149 (MQ=255)
                                                                                                        aCCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGAGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaa                                <  1:531004/149‑1 (MQ=255)
                                                                                                         ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaac                               >  1:526824/1‑148 (MQ=255)
                                                                                                             aTTCTCGGTCTGAAACCGACCGTTTTTTTTGCGGTCCAGATCATGTACACCGGCGTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAACTCATCTCCGAGgcagaagcaga                           <  2:663127/149‑1 (MQ=255)
                                                                                                              ttCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagcagaa                          <  1:1101271/149‑1 (MQ=255)
                                                                                                                ctcGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCAAGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGt                        <  2:367106/149‑1 (MQ=255)
                                                                                                                 tcGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGCAACACGAAATCATCTCCGAGGCAGAAGCAGAAGtt                       >  1:521343/1‑149 (MQ=255)
                                                                                                                  cGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTg                      <  1:1066981/149‑1 (MQ=255)
                                                                                                                  cGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTg                      >  1:1185475/1‑149 (MQ=255)
                                                                                                                   ggTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGc                     <  2:797605/149‑1 (MQ=255)
                                                                                                                    gTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCt                    <  1:471428/149‑1 (MQ=255)
                                                                                                                     tCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGAATATGCAGCGCGTTCTCGTGCATCTGTTGGTATCGAGGACATTGTCATCCCGTAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTg                   <  1:66739/149‑1 (MQ=255)
                                                                                                                       tGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaa                 >  1:1035667/1‑149 (MQ=255)
                                                                                                                       tGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCACGGAGAAGAAACACGAAATCATCTCCGTGGCAGAAGCAGAAGTTGCTGaa                 >  1:907578/1‑149 (MQ=255)
                                                                                                                        gAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaaa                <  1:346360/149‑1 (MQ=255)
                                                                                                                        gAAAACTACCGTTATTTTAGCTGTCCAGATCAAGTACACCGGCTTGGCCTATGCAGCGCGTTCTGGTGCATGTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaaa                <  1:268363/149‑1 (MQ=255)
                                                                                                                             cgaccgTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg           <  2:83291/149‑1 (MQ=255)
                                                                                                                              gaccgTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCACGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAgg          <  2:144705/149‑1 (MQ=255)
                                                                                                                               accgTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGa         >  1:783618/1‑149 (MQ=255)
                                                                                                                                ccgTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAg        <  1:774252/149‑1 (MQ=255)
                                                                                                                                  gTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCa      >  2:339300/1‑149 (MQ=255)
                                                                                                                                  gTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCa      <  1:349134/149‑1 (MQ=255)
                                                                                                                                   ttATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGCCCTGGTCATCCCGGAGAAGAAACCCGCAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAg     >  2:801934/1‑149 (MQ=255)
                                                                                                                                   ttATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAg     >  2:1270294/1‑149 (MQ=255)
                                                                                                                                      tttttGCGGTCCAGATCCTGTACACCGGCTTCGCCTAAGCAGCGCGTGCTGGTGCATCTGGGGGTATCGATGACATGTTCATCCCGGAGAAGATACACGACATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCCGTTc  <  2:816114/149‑1 (MQ=255)
                                                                                                                                      tttttGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTc  <  2:485570/149‑1 (MQ=255)
                                                                                                                                      tttttGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTc  <  2:857727/149‑1 (MQ=255)
                                                                                                                                      tttttGCGGTAGAGATCATGTACACCCTCTTCTCCTATGCAGCTAGTTATTGTGCTTCTGTTGGTATCGATGACATGGTCATCCCGTAGAAGAAACACGCAATCAGCTCCGAGGCAGAATCAGAAGTTGCTGAAATTCAGGAGCAGTTc  <  1:57249/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                           
TGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTC  >  NC_000913/4187071‑4187353

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: ‑