Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 164,396 C→T I764I (ATC→ATT)  hrpB → RNA‑dependent NTPase HrpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,3960CT96.9% 208.2 / ‑0.9 65I764I (ATC→ATT) hrpBRNA‑dependent NTPase HrpB
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (26/37);  total (27/38)
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.51e-01

CATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG  >  NC_000913/164150‑164642
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
catgtATCAGGCACTACGCGGATTACTTGCTTGGGGAATGCAGCAATGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAAGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGGGCTGGAGGTG‑CTTTCACCTTCCCAAAGGCCCTTACAAATTacac                                                                                                                                                                                                                                                    >  2:574011/5‑251 (MQ=255)
cATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTacac                                                                                                                                                                                                                                                    <  2:292126/251‑1 (MQ=255)
     aTCAGGCACAACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCTCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGGTTGGCGAGGCCACCAAGCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCGCCTGCCCAAATGCCATTACAAATTACACgaga                                                                                                                                                                                                                                                <  2:515809/250‑1 (MQ=255)
          gCACTACGCGGATTACTTGATTGGGGAATGCAGAAACGTCTGGCTAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACATGCCCAAAGGCCATTAAAAATTACACGTGATTTg                                                                                                                                                                                                                                            >  2:645837/1‑249 (MQ=255)
               aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCTCCCAGGGGCGCTTGCCGCTGGTGCTTGAGTTT‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGAc                                                                                                                                                                                                                                      <  2:447163/250‑1 (MQ=255)
                cgcgGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAGTTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGCCGATCGCCCAGGGGCGCGGGCCGCTGGTGCTGGCGTTG‑CTTTCACCTGCCCACAGGCCACTACCAATTACACGAGGTTTGAGCGACtt                                                                                                                                                                                                                                    >  1:277318/1‑251 (MQ=255)
                   ggATTACTTGATTGGGGAATGCCGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCTGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTg                                                                                                                                                                                                                                 >  1:137998/1‑251 (MQ=255)
                   ggATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTg                                                                                                                                                                                                                                 <  1:13400/251‑1 (MQ=255)
                     aTTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCTCATTACACTGGGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGCTAACCCGACCGCGCTGGCGGTGAGCATGCAAGAGAAGGTTGACGAGGCCACCAATCCGACGATCACCCAGGGGCGCGTGCCGCTGGGGCTGGAGTTG‑CTTTCACCTGCCCCAAGGCCATTACACATTACCCCAGATTTGAGCGACGTCTgg                                                                                                                                                                                                                                >  2:547239/1‑250 (MQ=255)
                        aCTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGaaa                                                                                                                                                                                                                             <  1:183235/250‑1 (MQ=255)
                        aCTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGaa                                                                                                                                                                                                                              >  1:253945/1‑249 (MQ=255)
                               tGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGt                                                                                                                                                                                                                      >  1:29106/1‑250 (MQ=255)
                                 gggAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTAcc                                                                                                                                                                                                                   <  1:63051/251‑1 (MQ=255)
                                               gTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGGGCaaaa                                                                                                                                                                                                      >  1:401658/1‑250 (MQ=255)
                                                       aGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGa                                                                                                                                                                                              <  2:600976/250‑1 (MQ=255)
                                                            aTTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATATCCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGc                                                                                                                                                                                        <  2:185629/251‑1 (MQ=255)
                                                                     gcATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATccc                                                                                                                                                                                >  1:219728/1‑250 (MQ=255)
                                                                                    gACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGCCTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACCTGTCTGGCCgg                                                                                                                                                                >  2:126031/1‑251 (MQ=255)
                                                                                      cGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGCCGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGgac                                                                                                                                                              >  1:589282/1‑251 (MQ=255)
                                                                                      cGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGgac                                                                                                                                                              >  1:72242/1‑251 (MQ=255)
                                                                                           aGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCACGGGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTGCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTACCCGAAACATGTCTGGGCGGACGAccc                                                                                                                                                         <  2:277318/251‑1 (MQ=255)
                                                                                            gCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCg                                                                                                                                                        <  2:29106/251‑1 (MQ=255)
                                                                                                     ccATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGAGTTCTGGAAAGGAGCGTACCGGGAGGTGTAAAAAGAGATGAAAGGGCGTTAGCCCAAACATGTCTGGCCGGACGACCCTGCAAAttct                                                                                                                                               <  2:61284/251‑4 (MQ=255)
                                                                                                         ccacgacgctcttccgagcccgcccgCGCTGGCGGTGAGAAGGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGTCAAATACTGCAc                                                                                                                                           <  2:489471/233‑1 (MQ=255)
                                                                                                            ttATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGAGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGCGATTTGAGCGACTTCTGGCAAGGAGCGTACCGAGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACcg                                                                                                                                         >  2:274769/1‑250 (MQ=255)
                                                                                                              agcATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACcgacg                                                                                                                                      <  1:694652/249‑1 (MQ=255)
                                                                                                              aTCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACcgac                                                                                                                                       <  1:126031/250‑1 (MQ=255)
                                                                                                                cATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATTCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGAGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGGGATTTGAGCGACTTCTGGAAAGGAGCGCACCGTGAGGTGAAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGACGGACGACCCGGCAAATACTGCACcgacgcg                                                                                                                                    >  1:679645/1‑251 (MQ=255)
                                                                                                                      gagaaCCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGAGTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTGATCCCAAACAt                                                                                                                                                                          <  2:120525/204‑1 (MQ=255)
                                                                                                                      gATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCAACTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACAt                                                                                                                                                                          >  1:120525/1‑207 (MQ=255)
                                                                                                                            ccgcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGaaaaa                                                                                                                        <  1:645837/251‑1 (MQ=255)
                                                                                                                             cgcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCAc                                                                                                                                           >  1:489471/1‑231 (MQ=255)
                                                                                                                                         cGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAggt                                                                                                           <  1:526460/251‑1 (MQ=255)
                                                                                                                                             gagaATGCAAGAGATGTTTGGCGAGGCCCCCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTAAggttgg                                                                                                        <  2:471878/250‑1 (MQ=255)
                                                                                                                                              agaATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAggttggtt                                                                                                      <  2:210804/251‑1 (MQ=255)
                                                                                                                                                 aTGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAGACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTtc                                                                                                   <  1:687420/251‑1 (MQ=255)
                                                                                                                                                     aagtgaggtgtGGCGAGGCCGCGCATTCGACGGTCGCCCAGGTGCGCGTGCGGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCGAGTATAAATTACACGAGAGTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGATATGAAAGGGCGGGATCCCAAACATGTCTGGCCGGACGACCCGGGATATGCTGTAGGGTCTCGACGGACGAAGAAGTATGCTTAAGGTTGGTTTtctccc                                                                                               <  2:622961/241‑1 (MQ=255)
                                                                                                                                                     aaGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTTGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTtctccc                                                                                               <  1:675064/251‑1 (MQ=255)
                                                                                                                                                     aaGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTtctccc                                                                                               <  1:536167/251‑1 (MQ=255)
                                                                                                                                                              ttGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTtctccctctccatg                                                                                       <  1:574011/250‑4 (MQ=255)
                                                                                                                                                                ggCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTgg                                                                                    >  2:345326/1‑251 (MQ=255)
                                                                                                                                                                   gAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGTACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGgag                                                                                 <  2:104253/251‑1 (MQ=255)
                                                                                                                                                                   gAGGCCACCAAGCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgag                                                                                 <  1:23787/251‑1 (MQ=255)
                                                                                                                                                                    aGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTggg                                                                                   >  2:227365/1‑248 (MQ=255)
                                                                                                                                                                    aGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                <  2:253945/251‑1 (MQ=255)
                                                                                                                                                                           caATCCGCCGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                         <  2:343562/251‑1 (MQ=255)
                                                                                                                                                                           caATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                         >  1:155808/1‑251 (MQ=255)
                                                                                                                                                                               ccGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAACGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg                                                                     >  2:139219/1‑251 (MQ=255)
                                                                                                                                                                                      ccGCCCAGGGGCGTGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCAAAAGGCCATTAAAAATTACACGAGATTTGAGCGACTTCGGGAAAGGAGCGTACCGTGAGGTCCAAAAAGAGATAAAAGGGCGTTATCCCAAACATGTCTGGCTGGTCGACCCGGAAAATCCTGCACCGACGCGACGGAGGAAAAAGTATGCGTATGGTGGGTTTTCCCCCGCTCCCTGGGGGAGGGGGCCGGGGTGAGGGCAt                                                              <  2:60005/250‑1 (MQ=255)
                                                                                                                                                                                               ggCGCGTTCCCCGGGTGCTGG‑GTTGGCTTTCACCTGCCCAAATGCCATTAAAAATTACACGAGATTTGGGCGATTTCGGGAAAGGAGCGTACCGTGGGGTGAAAAAAGGGAGGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCCCCGCCGCGGCGGACGAAAAAGTATGCGTAGGGTTGGTTTGCTCCCTCTCCCTGTGGGAGAGGGCCGGGTTGAGGGCATCAGCGCGCa                                                     <  2:476936/251‑1 (MQ=255)
                                                                                                                                                                                                 ggcgTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCAATTACAAATTCCACGAGATTTGAGCGCCTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                   <  1:500473/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggAGTTC‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                  <  1:26237/251‑1 (MQ=255)
                                                                                                                                                                                                                   gAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                 >  1:669700/1‑251 (MQ=255)
                                                                                                                                                                                                                           tttCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGCCTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGCGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCCCCCTCACCCTa                                    >  2:404026/1‑241 (MQ=255)
                                                                                                                                                                                                                           tttCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGGGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCCGGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTa                                    <  1:404026/241‑1 (MQ=255)
                                                                                                                                                                                                                             tCACCTGCCCAAAGTCCATAACAAATTAAACGAGATTTGAGCGACTTCTGGACAGGAGCGTACCGAGAGGTGCAAAAAGAGATGAAAGGGCGTTATCGCAAACATGTCTGGCCGGACGAACCGGCAAATACTGCACCGACGCGACGGACGAAAAAGAATTCGTAAGGTTGGTTTTCTACCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGAACGTTTACCCTCACCCTAACCCTCTCCCTc                        >  2:704794/1‑251 (MQ=255)
                                                                                                                                                                                                                             tCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTtctccctctcc                                                                                          >  1:849/1‑185 (MQ=255)
                                                                                                                                                                                                                             tCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTtctccctctcc                                                                                          <  2:849/185‑1 (MQ=255)
                                                                                                                                                                                                                                        aaGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCaa                                                                                                                                                    <  2:517557/116‑1 (MQ=255)
                                                                                                                                                                                                                                        aaGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCaa                                                                                                                                                    >  1:517557/1‑116 (MQ=255)
                                                                                                                                                                                                                                         aGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                      >  2:539155/1‑241 (MQ=255)
                                                                                                                                                                                                                                         aGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                      <  1:539155/241‑1 (MQ=255)
                                                                                                                                                                                                                                               ttACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCACTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGGGAGGGGACCGAt      >  2:88801/1‑251 (MQ=255)
                                                                                                                                                                                                                                                tACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt      <  1:462020/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaaTTACAGGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATAAAAGGTCGTGATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg  <  2:589282/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
CATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTG‑CTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG  >  NC_000913/164150‑164642

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