Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,535,189 C→T intergenic (‑1136/‑1408) topAI ← / → ythB toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,535,1890CT98.3% 177.4 / ‑4.2 59intergenic (‑1136/‑1408)topAI/ythBtoxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/32);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

GTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTAT  >  NC_000913/4534945‑4535430
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
gTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTa                                                                                                                                                                                                                                             <  1:587016/251‑1 (MQ=255)
   aCTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATg                                                                                                                                                                                                                                          <  2:209563/251‑1 (MQ=255)
              tATCCCGGTAAGTTTCAGCCTAAAATAGGTATAGAAAAAGGACGATTTTCCCCGCAAGGGGATAATTACATATTGATTGTAGAGGCATGAAGCCAAATGGAGTGAGCCGTCATAGCCCGCTTGAAGCTTAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAATTTATTGGTGTTGGGACGTCATTTTCTGCGTAGGTTAATGCTATCAAAAt                                                                                                                                                                                                                                <  2:442059/250‑1 (MQ=255)
              tATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGCGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATg                                                                                                                                                                                                                               >  1:405218/1‑251 (MQ=255)
                     gTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCGCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTTGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCAtt                                                                                                                                                                                                                          <  2:177069/249‑1 (MQ=255)
                              aggcTAAAAGAGGGAATGAAAAAGGACGATCGTTCCCGCAATGGTATTTTTACACATGGATGGTAGAGGCATGAAGACAAATGGAGGGAGTAGTCATAGCCCTGTTGAATCGGAATCAATATTAGCATGGCTGGCGGCACAGCAGGAGGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTGCGGCGTAGGTTAATGCTAGCAAAAGGGCATTGCGTAAGCTgg                                                                                                                                                                                                               <  2:476792/248‑1 (MQ=255)
                                      agagGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCACTCTTGAAGCTGAAACAATAGCAGCATGGCTGGGGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACAACTGTACAAAGTTGTTGGTGTTGTGAAGCAATTTTCGGCGTAGGTGAATGCTATCAAAATGTAATTGAGTAGGCTGGAAAGTAAt                                                                                                                                                                                                       >  2:328495/1‑251 (MQ=255)
                                        agGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATgg                                                                                                                                                                                                     >  1:332012/1‑251 (MQ=255)
                                              ggAAAAAGGACGATTTGCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTTTACAAAGTTGTTGGTGGTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAg                                                                                                                                                                                                <  2:220114/250‑1 (MQ=255)
                                                          aTTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGttt                                                                                                                                                                                   <  1:105824/251‑1 (MQ=255)
                                                                  cGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGAc                                                                                                                                                                             >  2:537858/1‑249 (MQ=255)
                                                                    cAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGACCCACTGTACACAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAGGCTGAAATTTAATGGCAAAGACGAGCAGGGTTTACTGACAGTg                                                                                                                                                                         >  2:400656/1‑251 (MQ=255)
                                                                            tatTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACtgt                                                                                                                                                                 >  1:178916/1‑251 (MQ=255)
                                                                                 aCATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACt                                                                                                                                                            <  1:537858/251‑1 (MQ=255)
                                                                                     atTGATGGTAGAGGCAAGAAGCCAAATGGAGGGTGCCATCATAGCCCACGTGAAGCTGAAACAATAGCATCATGGCTGGTTGCAAATAAGGATTATATAGAGCGCCATTACGGGGAACCCCTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTAGCAAAATGTCAGTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGTTTTACTGACAGTGGGGACTGTGCACTCCtt                                                                                                                                                        <  2:177826/251‑1 (MQ=255)
                                                                                        gATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGGTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTga                                                                                                                                                      <  2:391914/250‑2 (MQ=255)
                                                                                         agggTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCCTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAg                                                                                                                                                    <  2:47708/249‑1 (MQ=255)
                                                                                         aTGGTAGAGGCATGAAGCCAAATGGCGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGTAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGc                                                                                                                                                      >  2:156372/1‑249 (MQ=255)
                                                                                          tGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAgg                                                                                                                                                   >  1:644050/1‑251 (MQ=255)
                                                                                                       aaGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGc                                                                                                                                      <  1:452551/251‑1 (MQ=255)
                                                                                                       aaGCCAAATGGAGGGAGCAGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGc                                                                                                                                      >  1:113520/1‑251 (MQ=255)
                                                                                                            aaaTGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTg                                                                                                                                 >  1:228209/1‑251 (MQ=255)
                                                                                                              agggagtgaggcggcATAGCCCTCTTGAAGCTGAAACAATAGCAGTAGTGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTATAAAGTTGTTGGTGTTGTGACGCCATTTTCGGTGTAGGTTAATGCTATCAAAATTTCATTGCGTAATCTGTAAATTAATGGCAAAGACTAGAAGGGTTTACTGACAGTGTGGACTGTGCCCTCCGTGCAGGGGGCGGAAAGGGCGTTTGtt                                                                                                                               <  2:652751/237‑1 (MQ=255)
                                                                                                                  agggagCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAAGGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTctct                                                                                                                           <  1:400656/251‑1 (MQ=255)
                                                                                                                       gCCGACATAGCCATCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGTAATTTCATGGCAAAGAAGAGCAGGGTTTACTGCCAGTGGGGACTGTGCACTCCTAGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTtct                                                                                                                      >  2:126331/1‑251 (MQ=255)
                                                                                                                                 ccatctTTAAGCTGAAACAATAGCAGCATGGCTGGTGGTACATAAGGATCAGATAGAGCGCCATTACGGGGAACCACTGTACATAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGTACTGTGCACTCCTTGCAGGTGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCtcct                                                                                                            <  2:124785/248‑4 (MQ=255)
                                                                                                                                  ccTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTc                                                                                                           >  1:298047/1‑251 (MQ=255)
                                                                                                                                  ccTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGAGGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTGTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCa                                                                                                          <  2:669418/251‑1 (MQ=255)
                                                                                                                                        gAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGCCTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAAc                                                                                                      >  2:675318/1‑250 (MQ=255)
                                                                                                                                        gAAGCTGAAACAATAGCAGCAGGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCCCGGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAAc                                                                                                      <  2:194585/250‑1 (MQ=255)
                                                                                                                                              gaaccaagagcagcATGGCTGGTGGCATATACGGTTGATTTAGAGCGCCATTACGGGTAACCGCTGTACAAAGTTGTTGGTGTTGTGACGTCATGTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGTAAATTAATGGCAAAGGCGAGCAGGGTTGACTGACAGTGGGGTCTGTGCTCGCCTTGCAGGTGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGGTTACTCAAAACATGagg                                                                                                <  2:156016/242‑3 (MQ=255)
                                                                                                                                              gAAAAAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGaa                                                                                                 <  2:178916/249‑1 (MQ=255)
                                                                                                                                                 aCAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTTCGGGGAACCACTGTACAAAGTTGTTGGTATTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATgg                                                                                            <  1:610841/251‑1 (MQ=255)
                                                                                                                                                   aaTAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGc                                                                                           >  1:406527/1‑250 (MQ=255)
                                                                                                                                                    aTAGCAGCATGGATGGAGGCACCTAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGCTGTTGGTGTTGTGACGCCATTTTAGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGTGCCGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGAGATTGTTCTCTTTTCTCCTGTCTACTCAACACATGAAGATGGc                                                                                           >  2:70951/1‑249 (MQ=255)
                                                                                                                                                                     ggCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCGTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACa                                                                        <  1:312998/251‑1 (MQ=255)
                                                                                                                                                                      gCACCTAAGGATGATATAGAGCGACATTACGGGGAACCACTGTACAAAGTTGGTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCCGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTTTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAg                                                                       <  2:9426/251‑1 (MQ=255)
                                                                                                                                                                        acaTAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACAt                                                                                                    <  2:217555/220‑1 (MQ=255)
                                                                                                                                                                        acaTAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACAt                                                                                                    >  1:217555/1‑220 (MQ=255)
                                                                                                                                                                          aTAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACt                                                                   >  1:395836/1‑251 (MQ=255)
                                                                                                                                                                          aTAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACt                                                                   >  1:325612/1‑251 (MQ=255)
                                                                                                                                                                               gatgatATAGAGCGCCATTACGGGGAACCACTGTACACAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGCAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCAc                                                                    >  2:278483/1‑245 (MQ=255)
                                                                                                                                                                               gatgatATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCAc                                                                    <  1:278483/245‑1 (MQ=255)
                                                                                                                                                                                 ggatATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCAATTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCa                                                            <  2:332012/250‑1 (MQ=255)
                                                                                                                                                                                     ataGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGt                                                        >  1:174332/1‑251 (MQ=255)
                                                                                                                                                                                       agggcgcCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGtt                                                       <  1:575010/247‑1 (MQ=255)
                                                                                                                                                                                       agagCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAACCATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTg                                                      >  2:621474/1‑251 (MQ=255)
                                                                                                                                                                                       agagCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTg                                                      >  1:572511/1‑251 (MQ=255)
                                                                                                                                                                                             ccATTACGGGGAACCGCTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGtct                                                <  1:341241/251‑1 (MQ=255)
                                                                                                                                                                                                ttACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCAc                                             >  1:283542/1‑251 (MQ=255)
                                                                                                                                                                                                  aCGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCATTATCCTCAATGTTGCTGTCTCAcg                                            <  1:240990/250‑1 (MQ=255)
                                                                                                                                                                                                   cGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCAcgcg                                          <  2:228209/251‑1 (MQ=255)
                                                                                                                                                                                                    ggggAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTAGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACgcgc                                         >  1:162706/1‑251 (MQ=255)
                                                                                                                                                                                                               ggaCAAAGTTGTTGGTGTTGTGCCGCCTTGTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGAGCTAAAGATAGtt                              <  2:471365/249‑1 (MQ=255)
                                                                                                                                                                                                               gTACAAAGTTGTTGCTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTACCAACATGTCATTGCGTAAGCTGGACATTAATGGCAACGACGAGCAGGGTTTACTAACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTCCTCCAACCACGAAGATGGCAGGTTTCTTGTTAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATGGtt                              >  2:440184/1‑251 (MQ=255)
                                                                                                                                                                                                                  cAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTc                           >  1:511604/1‑251 (MQ=255)
                                                                                                                                                                                                                               gttgtGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCg               <  1:377333/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gtgACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGatat           <  1:14180/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tttCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGGAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTGCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCGCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGTCCTTAt  <  2:113520/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
GTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTAT  >  NC_000913/4534945‑4535430

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