Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 989,959 A→C 18.8% V133G (GTA→GGA)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,9590AC18.8% 122.5 / 4.6 71V133G (GTA→GGA) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (33/21);  new base C (2/12);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 2.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAA  >  NC_000913/989723‑990196
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
gggAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTAc                                                                                                                                                                                                                                                 <  1:992390/238‑1 (MQ=255)
     aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACCAATTTGCTTTCCAGCGTGTCGGAGGGCTGGGCCAACTCGGCCTGCGGTGAGCGCTAGCGCTGTACCActt                                                                                                                                                                                                                                            >  2:897698/1‑238 (MQ=255)
     aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTAGACTAATTTGCTTTCCAGCGTGTCGAGGGGTTGCGCCAAGTCGGCCTGCGGGGAGCGAGAGCGATGTAcccattt                                                                                                                                                                                                                                           >  1:148373/1‑234 (MQ=255)
                    tCGTATTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGACGTTTAACTATGGTTTCTTTTACTACGAGAAAAAAGCGGCGACGGGTGCCAAAACCCATCAGATGGAAGCGCGACATATCAAGCCCGGGGGTTAACGCCGAGAGGGCGACTAATTTGCTTCCCAGCGTGTCGAGGGTTGGCGCAACGCCGGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGcc                                                                                                                                                                                                                            <  2:270343/237‑1 (MQ=255)
                            gcggtggCCCACAAACCAGGATTCCGGTTGCAGACGTTTAACGATGTTTTCTTTTACTTCGCGAAAAAAAGGGCGACGGGTCCCAAAATCCACCAGATGGAAGCCCGACATACAAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGaa                                                                                                                                                                                                                  <  2:148373/233‑1 (MQ=255)
                             cTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCCGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGGGATGGACCATTTTACCGATAACCGCCAGCAAAAGa                                                                                                                                                                                                                   >  1:1209956/1‑237 (MQ=255)
                              tGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGGGGGCTTGCGCCACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTTACTGATAACCCCCCGCACAAGAACt                                                                                                                                                                                                                >  2:428772/1‑239 (MQ=255)
                               ggTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAATTGCTTTCCAGCGTGTCGGGGGGGTGCGCAAACTCGGCCTGCGGGGGGCGATAGCGCTGTCGCAATTTACTGATAAACGCCCACCAAGGGAACt                                                                                                                                                                                                                >  2:44251/1‑239 (MQ=255)
                                  gCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGGACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTACGCCGACAAATCGACTAATTTTCTTTCCCGCCTGGTGAGGGCTTAAGCAAAGTCCGCCTGCGGTGAGCGAAAAGGGTTGACCCTTTTCCTGAGACACGCCAGCAAAGGGAATTcgc                                                                                                                                                                                                            >  2:649585/1‑237 (MQ=255)
                                  gCCCACAAACCAGGATTCCTGTGGCAGACGCTTAACGATGTTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAAGCCATCAGATGGAAGCGCGACATATCAAGCCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTcc                                                                                                                                                                                                             <  1:1212552/238‑1 (MQ=255)
                                       caAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTTCAACGGGGCGAGGGCCTGCGCAACGCCGGCCCGCGGTGAGCGATAACGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCCGAgg                                                                                                                                                                                                       >  1:136552/1‑239 (MQ=255)
                                         aaCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCAAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTACCGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGAt                                                                                                                                                                                                     <  1:428772/239‑1 (MQ=255)
                                          aCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAAGTGACTAATTTGCTTTTCAACCGGTCGAGGGGCTGCGCAAACGTGGGCCGCGGGGGGTGGAAGAGGTGTACCATTTAACTGATAACCCCCAGGAAAGAAACATCCCAGAGGaag                                                                                                                                                                                                    >  2:934169/1‑237 (MQ=255)
                                               gATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTTCAACGTGTCGAGGGCTTGCGCAACGTCGGGCCGCGGGGGGCGATAGCGGTTTACCATTTCACCGATAAACGCCCGCAAAGGAAATTTCCCGAGGGTGACtt                                                                                                                                                                                                >  2:519363/1‑238 (MQ=255)
                                                aTTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGGGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                                               <  1:484931/238‑1 (MQ=255)
                                                 tttccgTTGCAGGCCGCTTAAGGATGGTTTCTTCTACGTCGCGAAAAAACCGGCGACGGGTGCAAAAAGCCATAAGATGGAAGGCCGAAATATCAAGACCGGGGTTTAACGCCGAAAAGTCGATTAATTTGTTTCCCAGCGTGTCGAGGGCTTGCGCACCGTCGGCCTGCGGTGGGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                                               <  2:784158/233‑1 (MQ=255)
                                                    cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTTCTTTTCAGCGGGTCGAGGGCTTGCGCAACGGCGGCCTGCGGTGAGCGGTAACGATGTAACATTTTACTTATAAACGCCCGCAAAGGAACTTCCCCGAGGGTGACTTTACGc                                                                                                                                                                                           >  1:532900/1‑238 (MQ=255)
                                                    cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCCAGAAGTCGACTCATTTTCTTTCCCGCGGGTCGGGGGGTTGCGCAACGCCCGCCCGCGGGGGGGGGTAGCGGTGTACCATTTTCATGATAAACCCCCGCAAAAGGACATCCCAGGGGGAGACtttccgc                                                                                                                                                                                           >  2:182688/1‑233 (MQ=255)
                                                                           ttttCTTGTACTGCGGAAGAAAAAGGGCGAGGGGTGCCAAAATCCACAGAAGGGAAGCGCGACATATCAAGACGGGCGTTTAACGCGGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGAAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAatat                                                                                                                                                                     <  1:967773/237‑1 (MQ=255)
                                                                                 tGTCCTTGGGGAAAAAAGCGGGGACGGGGCCAAAAACCCATCAGATGGAAGCCGGACATATCAAGCCCGGCGGTTAACGCGGAGAAGTCGACTAATTGGCTTTCCAGCGTGTCGAGGGCTTGCGAACCGTCGGCCTGCGGGGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCtt                                                                                                                                                              <  2:961553/239‑1 (MQ=255)
                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTCGGGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCATTTCAATGATAACCGCCCGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGc                                                                                                                                                                                    >  1:1235500/1‑210 (MQ=255)
                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCCGCGGTGGGCGATAGCGATGTAACATTTCACTGCTAAACGCCAGCAAAAGAACTTCCCCCAGGATGACTTCACGCCACCGGCCGCTTAAACGCATATCCAACTTTCCCttt                                                                                                                                                       >  1:1001106/1‑238 (MQ=255)
                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGGACAACAAATTGCTTTACAGAGGGTCGGAGGGCTGCGCGAAATCGGGCTTCAGGGGGCGATAGACATTTCACATTTCCTCGCAAAACGCCCGCCAAGGGACTTACCAGAAGAAGGGTTCCAGCCACCGGGCGCCTAAACCGATATTCCACTTGCCACTa                                                                                                                                                       >  2:210005/1‑239 (MQ=255)
                                                                                                aaccGGCGAGGGGGGCAAAAACCCATCAGATGGAAGCGCGACATACCAAGACCGGCGGTAAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGGGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGa                                                                                                                                               <  1:296234/236‑1 (MQ=255)
                                                                                                              gcaaaaacacATCAGATGGAAGCGCGACATATCAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGCTTCCCGGCGTGTCGAGGGCTTGCGCAACGTCGGCCGGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGtt                                                                                                                                <  2:532900/230‑1 (MQ=255)
                                                                                                              gacAAAACCCATAAAGTGGAAGGCCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTCTTTTCCAGCGTGTCGGGGGTTGGCGCACCGCCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGtt                                                                                                                                <  2:849192/237‑1 (MQ=255)
                                                                                                               ccAAAATCCATCAGAGAAAAGCGGAAAATATAAGGACCGGCTGTTAACGCGAATAAGTCGACTATTTGGCTTCCCGGCGTGTGTGGGGCTTGCGCACCGTGCGCCGGCGGTGGTCGAGAGCGAGTTCCCATTTCACTGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGtt                                                                                                                                <  2:943213/238‑1 (MQ=255)
                                                                                                               ccAAAACCCATCAGATGGAAGCGCGACATATCAAGGCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGtt                                                                                                                                <  1:519363/238‑1 (MQ=255)
                                                                                                                         tCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCCTTTCACTGATAACCGCCCGCCAAGGAACTTCCCAAGAGATGACTTCACGCCAACGGCCGCTTAAACGAAAATCCAGCTTGCCATTATCATTGGACACAGGGGCTTGTTCCGGGGTaaa                                                                                                                      >  2:572585/1‑238 (MQ=255)
                                                                                                                                ggAAGCGCGAAAGATCAAGACCGGGGGTTAACGCCGAGAAGCCGACTAATTTGCTTCCCAGCGTGTCGAGGTCTTGCCCAACGTCGGCCGGGGGTGAGCAATAGCGATGTACCTTTCCCCTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGaa                                                                                                                <  2:1209956/238‑1 (MQ=255)
                                                                                                                                        aaCATATCAAGCCCGGGGGTTAACGCCGGGAGGTCGACTATTTTGTTTCCCAGCGTGTCGGGGGTTTGCGCAACGGCGGCCTGCGGGGGGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                                       <  2:289057/238‑1 (MQ=255)
                                                                                                                                         aCATATCAAGCCCGGCGGTTAACGCGGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                                       <  1:44251/238‑1 (MQ=255)
                                                                                                                                            tatCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAATTTCTTTCCACCGTGTCGAGGGCTTGCGCCACGTCGGCCTGCGGTGAGCGATAGCGATGTACCCTTTCACTGATAACCGCCCGCAAACGCACTTCCCCGAGGATGACTTCCCCCCACCGGCCGCTTAAACGAATATCCCGCCTGCCCTTATTCTTGTACACGGTGGCTTGTTTCCGGGTTAAACGGACAACGCGGAACCac                                                                                                   >  2:689178/1‑238 (MQ=255)
                                                                                                                                               cAAGCCCGGGGTTTAACGTCGAAACGTCGGCTAGTTGGCTTCCCAGGGTGTGGAGGGTTTGCCCAAGGTCGGCGCGGGGGGACGGATAGCGTGGTCCCATTTCTCGAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTa                                                                                                <  2:1226803/239‑1 (MQ=255)
                                                                                                                                                      gggggTAAACGCCGGGAAGGCGAATAATTTGCTTTCCAGCTGGTCGAGGGCATGCCAACCGTCGGCTGGCGGTGACGAATAGCGATGTACAATTGCCCTGAGAACCGCCAGCAAAGGAACTGCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCGTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                                         <  2:157560/236‑1 (MQ=255)
                                                                                                                                                       gcggTTTAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTCGGGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCCTTTTACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCCCGGGCCGCTTAAACGAAAATCCAGCTTGCCATTATCGTTGGACACGGGGGATTGTTCCGGGGTAAAGCGGAACTCCGGGAAACAGTTaaaaaaaac                                                                                        >  2:686493/1‑232 (MQ=255)
                                                                                                                                                       gcggTTAACGCCGAGAAGTCGACTAATTTTCTTTCCAGCGTGTCGGGGGCTTGCGCAACGGCGGCCTGCGGGGAGCGAGAGCGGTGTACCATTTTACTGATAAACGCCCGCCAAGGAACTTCCCCGAGGGTGACTTTCCCCCACCGGCCGCTTAAACGAAAATACACCTTGCCCTTATCGTTTGACACGGTGACCTGTTTCGGGTTCAAGCGGAACTCCCGTAAACCGTTACGATAGt                                                                                         >  2:631283/1‑238 (MQ=255)
                                                                                                                                                        cggTTACCGCCGAGAAGTCGACTAATTTGCTTCCCAGGGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGAGATTTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                                        <  1:740722/238‑1 (MQ=255)
                                                                                                                                                        cggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCGGCGTTGCGTGTGCTTGCGAACCGGCTGCGTGCGGGGATCGATAGCGATGTCCCATTTCCCTAAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                                        <  2:136552/238‑1 (MQ=255)
                                                                                                                                                                    agaAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                                           >  2:379034/1‑239 (MQ=255)
                                                                                                                                                                    agaAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                                           >  2:116645/1‑239 (MQ=255)
                                                                                                                                                                      aaGTCGACTAATTGGCTTTCCAGCGTGTCGGGGGCTTGCCAACCGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                                         <  1:901034/239‑1 (MQ=255)
                                                                                                                                                                       agtcatagaaTTTGTTTCCCACTGTGCCGAGGGTTTGCGCCCCGTGGGCCTGCGGGGAGAGTTAGTGATGTACCATTGCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                                         <  1:1158402/230‑1 (MQ=255)
                                                                                                                                                                        gTCGGCTAATTTGCTTTCGTTGGTGGGGAGGTCTTGCGCATGGCCGCCCTGCGTGGAGGGAGAGCGATGTCCTTTCTCATTGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggc                                                                        <  2:379445/238‑1 (MQ=255)
                                                                                                                                                                         tCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcag                                                                      >  2:1218794/1‑239 (MQ=255)
                                                                                                                                                                            aCTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcagggc                                                                   >  1:204556/1‑239 (MQ=255)
                                                                                                                                                                                   tgcttccgggCGTGTGGGGGGCTTGCGAACGGCGGGCGAGCGGTGGGGGATGGCGGTGTTCCATTTCGCTAATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                                                            <  2:20286/230‑1 (MQ=255)
                                                                                                                                                                                       ttCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGcc                                                           >  1:707889/1‑236 (MQ=255)
                                                                                                                                                                                          cAGTGTGTCGGGGGTTTCGGCAATCGCGGCCTGCGTGGGCCAATAGCTATGTACCATTTCACTTATATCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGat                                                     <  1:865714/239‑1 (MQ=255)
                                                                                                                                                                                              gtgtCGAGGGTTTGCGCAACGCCGGCCGGCGGGGAGCGATAGCGATGTACAATTTCACAGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                                                  <  1:416752/238‑1 (MQ=255)
                                                                                                                                                                                              gtggCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                                                  <  1:1218794/238‑1 (MQ=255)
                                                                                                                                                                                               tgtCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtca                                                >  1:819315/1‑239 (MQ=255)
                                                                                                                                                                                                           tGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                                     >  1:731040/1‑239 (MQ=255)
                                                                                                                                                                                                            gcgcAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                                     >  1:1091819/1‑238 (MQ=255)
                                                                                                                                                                                                            gcgcAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                                     >  2:249316/1‑238 (MQ=255)
                                                                                                                                                                                                                         gCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgct                          >  2:334793/1‑238 (MQ=255)
                                                                                                                                                                                                                           gggcgtgaggggatagcgatGTCCCTTTCACCGAAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:548900/228‑1 (MQ=255)
                                                                                                                                                                                                                           cggCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:1152536/237‑1 (MQ=255)
                                                                                                                                                                                                                           cgcgggggAGCGGAAGCGATTTAACTTTCCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:600998/232‑1 (MQ=255)
                                                                                                                                                                                                                           cTGGGGTGAGCGATACCGTTGTGCCATTTCACTAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:680435/239‑1 (MQ=255)
                                                                                                                                                                                                                           cTGCGGTGGGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:182688/239‑1 (MQ=255)
                                                                                                                                                                                                                           cTGCGGTGAGGGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgctgc                        <  1:1117046/238‑1 (MQ=255)
                                                                                                                                                                                                                              cGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc                     >  2:678394/1‑238 (MQ=255)
                                                                                                                                                                                                                              cGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc                     >  1:467246/1‑238 (MQ=255)
                                                                                                                                                                                                                                   agcgatagcgatGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaa                                                                               >  1:490444/1‑175 (MQ=255)
                                                                                                                                                                                                                                        tagcgatGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGaa              >  1:988364/1‑235 (MQ=255)
                                                                                                                                                                                                                                           cgatGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTc       >  2:1020312/1‑239 (MQ=255)
                                                                                                                                                                                                                                            gttGTACCATTTCCCTGACAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTc       <  2:386943/236‑1 (MQ=255)
                                                                                                                                                                                                                                             atGTACCTTTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                                                                                   <  1:747319/161‑1 (MQ=255)
                                                                                                                                                                                                                                             atGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                                                                                   >  2:747319/1‑161 (MQ=255)
                                                                                                                                                                                                                                                tACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGaa  >  1:290137/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 aCCATTTCACTGAAAACCGCCATCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGaa  <  2:669835/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
GGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAA  >  NC_000913/989723‑990196

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