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→ATThrpBRNA‑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

CATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG  >  NC_000913/164150‑164642
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
catgtATCAGGCACTACGCGGATTACTTGCTTGGGGAATGCAGCAATGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAAGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGGGCTGGAGGTGCTTTCACCTTCCCAAAGGCCCTTACAAATTacac                                                                                                                                                                                                                                                    >  2:574011/5‑251 (MQ=255)
cATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTacac                                                                                                                                                                                                                                                    <  2:292126/251‑1 (MQ=255)
     aTCAGGCACAACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCTCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGGTTGGCGAGGCCACCAAGCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCGCCTGCCCAAATGCCATTACAAATTACACgaga                                                                                                                                                                                                                                                <  2:515809/250‑1 (MQ=255)
          gCACTACGCGGATTACTTGATTGGGGAATGCAGAAACGTCTGGCTAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACATGCCCAAAGGCCATTAAAAATTACACGTGATTTg                                                                                                                                                                                                                                            >  2:645837/1‑249 (MQ=255)
               aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCTCCCAGGGGCGCTTGCCGCTGGTGCTTGAGTTTCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGAc                                                                                                                                                                                                                                      <  2:447163/250‑1 (MQ=255)
                cgcgGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAGTTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGCCGATCGCCCAGGGGCGCGGGCCGCTGGTGCTGGCGTTGCTTTCACCTGCCCACAGGCCACTACCAATTACACGAGGTTTGAGCGACtt                                                                                                                                                                                                                                    >  1:277318/1‑251 (MQ=255)
                   ggATTACTTGATTGGGGAATGCCGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCTGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTg                                                                                                                                                                                                                                 >  1:137998/1‑251 (MQ=255)
                   ggATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTg                                                                                                                                                                                                                                 <  1:13400/251‑1 (MQ=255)
                     aTTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCTCATTACACTGGGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGCTAACCCGACCGCGCTGGCGGTGAGCATGCAAGAGAAGGTTGACGAGGCCACCAATCCGACGATCACCCAGGGGCGCGTGCCGCTGGGGCTGGAGTTGCTTTCACCTGCCCCAAGGCCATTACACATTACCCCAGATTTGAGCGACGTCTgg                                                                                                                                                                                                                                >  2:547239/1‑250 (MQ=255)
                        aCTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGaaa                                                                                                                                                                                                                             <  1:183235/250‑1 (MQ=255)
                        aCTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGaa                                                                                                                                                                                                                              >  1:253945/1‑249 (MQ=255)
                               tGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGt                                                                                                                                                                                                                      >  1:29106/1‑250 (MQ=255)
                                 gggAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTAcc                                                                                                                                                                                                                   <  1:63051/251‑1 (MQ=255)
                                               gTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGGGCaaaa                                                                                                                                                                                                      >  1:401658/1‑250 (MQ=255)
                                                       aGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGa                                                                                                                                                                                              <  2:600976/250‑1 (MQ=255)
                                                            aTTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATATCCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGc                                                                                                                                                                                        <  2:185629/251‑1 (MQ=255)
                                                                     gcATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATccc                                                                                                                                                                                >  1:219728/1‑250 (MQ=255)
                                                                                    gACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGCCTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACCTGTCTGGCCgg                                                                                                                                                                >  2:126031/1‑251 (MQ=255)
                                                                                      cGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGCCGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGgac                                                                                                                                                              >  1:589282/1‑251 (MQ=255)
                                                                                      cGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGgac                                                                                                                                                              >  1:72242/1‑251 (MQ=255)
                                                                                           aGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCACGGGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTGCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTACCCGAAACATGTCTGGGCGGACGAccc                                                                                                                                                         <  2:277318/251‑1 (MQ=255)
                                                                                            gCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCg                                                                                                                                                        <  2:29106/251‑1 (MQ=255)
                                                                                                     ccATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGAGTTCTGGAAAGGAGCGTACCGGGAGGTGTAAAAAGAGATGAAAGGGCGTTAGCCCAAACATGTCTGGCCGGACGACCCTGCAAAttct                                                                                                                                               <  2:61284/251‑4 (MQ=255)
                                                                                                         ccacgacgctcttccgagcccgcccgCGCTGGCGGTGAGAAGGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGTCAAATACTGCAc                                                                                                                                           <  2:489471/233‑1 (MQ=255)
                                                                                                            ttATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGAGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGCGATTTGAGCGACTTCTGGCAAGGAGCGTACCGAGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACcg                                                                                                                                         >  2:274769/1‑250 (MQ=255)
                                                                                                              agcATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACcgacg                                                                                                                                      <  1:694652/249‑1 (MQ=255)
                                                                                                              aTCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACcgac                                                                                                                                       <  1:126031/250‑1 (MQ=255)
                                                                                                                cATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATTCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGAGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGGGATTTGAGCGACTTCTGGAAAGGAGCGCACCGTGAGGTGAAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGACGGACGACCCGGCAAATACTGCACcgacgcg                                                                                                                                    >  1:679645/1‑251 (MQ=255)
                                                                                                                      gagaaCCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGAGTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTGATCCCAAACAt                                                                                                                                                                          <  2:120525/204‑1 (MQ=255)
                                                                                                                      gATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCAACTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACAt                                                                                                                                                                          >  1:120525/1‑207 (MQ=255)
                                                                                                                            ccgcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGaaaaa                                                                                                                        <  1:645837/251‑1 (MQ=255)
                                                                                                                             cgcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCAc                                                                                                                                           >  1:489471/1‑231 (MQ=255)
                                                                                                                                         cGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAggt                                                                                                           <  1:526460/251‑1 (MQ=255)
                                                                                                                                             gagaATGCAAGAGATGTTTGGCGAGGCCCCCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTAAggttgg                                                                                                        <  2:471878/250‑1 (MQ=255)
                                                                                                                                              agaATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAggttggtt                                                                                                      <  2:210804/251‑1 (MQ=255)
                                                                                                                                                 aTGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAGACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTtc                                                                                                   <  1:687420/251‑1 (MQ=255)
                                                                                                                                                     aagtgaggtgtGGCGAGGCCGCGCATTCGACGGTCGCCCAGGTGCGCGTGCGGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCGAGTATAAATTACACGAGAGTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGATATGAAAGGGCGGGATCCCAAACATGTCTGGCCGGACGACCCGGGATATGCTGTAGGGTCTCGACGGACGAAGAAGTATGCTTAAGGTTGGTTTtctccc                                                                                               <  2:622961/241‑1 (MQ=255)
                                                                                                                                                     aaGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTTGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTtctccc                                                                                               <  1:675064/251‑1 (MQ=255)
                                                                                                                                                     aaGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTtctccc                                                                                               <  1:536167/251‑1 (MQ=255)
                                                                                                                                                              ttGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTtctccctctccatg                                                                                       <  1:574011/250‑4 (MQ=255)
                                                                                                                                                                ggCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTgg                                                                                    >  2:345326/1‑251 (MQ=255)
                                                                                                                                                                   gAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGTACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGgag                                                                                 <  2:104253/251‑1 (MQ=255)
                                                                                                                                                                   gAGGCCACCAAGCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgag                                                                                 <  1:23787/251‑1 (MQ=255)
                                                                                                                                                                    aGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTggg                                                                                   >  2:227365/1‑248 (MQ=255)
                                                                                                                                                                    aGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                <  2:253945/251‑1 (MQ=255)
                                                                                                                                                                           caATCCGCCGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                         <  2:343562/251‑1 (MQ=255)
                                                                                                                                                                           caATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                         >  1:155808/1‑251 (MQ=255)
                                                                                                                                                                               ccGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAACGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg                                                                     >  2:139219/1‑251 (MQ=255)
                                                                                                                                                                                      ccGCCCAGGGGCGTGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCAAAAGGCCATTAAAAATTACACGAGATTTGAGCGACTTCGGGAAAGGAGCGTACCGTGAGGTCCAAAAAGAGATAAAAGGGCGTTATCCCAAACATGTCTGGCTGGTCGACCCGGAAAATCCTGCACCGACGCGACGGAGGAAAAAGTATGCGTATGGTGGGTTTTCCCCCGCTCCCTGGGGGAGGGGGCCGGGGTGAGGGCAt                                                              <  2:60005/250‑1 (MQ=255)
                                                                                                                                                                                               ggCGCGTTCCCCGGGTGCTGGGTTGGCTTTCACCTGCCCAAATGCCATTAAAAATTACACGAGATTTGGGCGATTTCGGGAAAGGAGCGTACCGTGGGGTGAAAAAAGGGAGGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCCCCGCCGCGGCGGACGAAAAAGTATGCGTAGGGTTGGTTTGCTCCCTCTCCCTGTGGGAGAGGGCCGGGTTGAGGGCATCAGCGCGCa                                                     <  2:476936/251‑1 (MQ=255)
                                                                                                                                                                                                 ggcgTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCAATTACAAATTCCACGAGATTTGAGCGCCTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                   <  1:500473/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggAGTTCCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                  <  1:26237/251‑1 (MQ=255)
                                                                                                                                                                                                                   gAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATTACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                 >  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
CATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG  >  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: