Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,031,221 T→A 9.8% intergenic (‑290/‑146) ygfT ← / → uacT putative oxidoreductase, Fe‑S subunit/nucleotide‑binding subunit/uric acid permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,031,2210TA 9.8% 266.9 / 17.3 82intergenic (‑290/‑146)ygfT/uacTputative oxidoreductase, Fe‑S subunit/nucleotide‑binding subunit/uric acid permease
Reads supporting (aligned to +/- strand):  ref base T (41/33);  new base A (4/4);  total (45/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

GCATATATTAGGTAGCATGACCGTTTTAATTATAATTTTGCGACTGAGTTCAAATTATTGCCCCTCAATCTATGCCGCTTTCATTATTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGGTTGACCGCATATTATCACCAGGAAAGCTGATCATACTCGGTCTGCAACACGTCCTTGTCATGTACG  >  NC_000913/3030990‑3031493
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                
gCATATATTAGGTAGCATGACCGTTTTAATTATAATTTTGCGACTGAGTTCAAATTATTGCCCCTCAATCTATGCCGCTTTCATTATTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGACTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATCCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCCATTCCTGTCTCTAACt                                                                                                                                                                                                                                       >  2:17147/1‑275 (MQ=255)
      atTAGGTAGCATGACCGTTTTAATTATAATTTTGCGACTGAGTTCAAATTATTGCCCCTCAATCTATGCCGCTTTCATTATTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGAGACATGCCAGAATGATTAATGAAATATCAATATGAGAAACTTACAGATAATTCCATTGTCAGGTTCTCTTCTAATTTTCGTGGTTTAACCACGATTAGGTGTTTGCCTCTGCAATTCCTGTCTCTAACTCCCCt                                                                                                                                                                                                                                  >  4:395297/1‑274 (MQ=255)
               cATGACCGTTTTAATTATAATTTTGCGACTGAGTTCAAATTATTGCCCCTCAATCTATGCCGCTTTCATTATTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCaa                                                                                                                                                                                                                         >  2:127335/1‑274 (MQ=255)
                                  aTTTTGCGACTGAGTTCAAATTATTGCCCCTCAATCTATTCCGCTTTCATTAGTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCa                                                                                                                                                                                                     <  1:96320/275‑1 (MQ=255)
                                                               cTCAATCTATGCCGCTTTCATTATTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAAc                                                                                                                                                                        <  2:479040/275‑1 (MQ=255)
                                                                                                    aTCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGATCTCTTCTCATTTTCGTGGTTTACCCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTCCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAacac                                                                                                                                   >  2:193526/1‑275 (MQ=255)
                                                                                                     tCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGAt                                                                                                                                                                                                                                                                 <  6:1072749/148‑1 (MQ=255)
                                                                                                          agagCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCt                                                                                                                                                                                                                                                            <  7:1033751/148‑1 (MQ=255)
                                                                                                                        aTAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtc                                                                                                                                                                                                                                              >  5:451899/1‑148 (MQ=255)
                                                                                                                         tAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtct                                                                                                                                                                                                                                             <  6:451899/148‑1 (MQ=255)
                                                                                                                          aaaaaCGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTg                                                                                                                                                                                                                                                <  7:385559/144‑1 (MQ=255)
                                                                                                                          aaaaaCGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTg                                                                                                                                                                                                                                                >  8:385559/1‑144 (MQ=255)
                                                                                                                                                   aCATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                                                                                                                                                                   <  8:517275/148‑1 (MQ=255)
                                                                                                                                                   aCATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTAAACCACGaa                                                                                                                                                                                                                                                                         >  7:559093/1‑93 (MQ=255)
                                                                                                                                                   aCATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTAAACCACGaa                                                                                                                                                                                                                                                                         <  8:559093/94‑2 (MQ=255)
                                                                                                                                                             aTGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACg                                                                                                                                                                                                         >  7:963920/1‑148 (MQ=255)
                                                                                                                                                                   aaTGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAgg                                                                    >  4:57296/1‑275 (MQ=255)
                                                                                                                                                                     tGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtctc                                                                                                                                                                                                                                            <  8:837325/105‑1 (MQ=255)
                                                                                                                                                                     tGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtctc                                                                                                                                                                                                                                            >  7:837325/1‑105 (MQ=255)
                                                                                                                                                                        aaTATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGttt                                                                                                                                                                                              <  8:963920/148‑1 (MQ=255)
                                                                                                                                                                                 atGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGAc                                                                                                                                                                                          <  7:9815/143‑1 (MQ=255)
                                                                                                                                                                                 atGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGAc                                                                                                                                                                                          >  8:9815/1‑143 (MQ=255)
                                                                                                                                                                                              aGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtct                                                                                                                                                                                                                                             <  8:765079/79‑1 (MQ=255)
                                                                                                                                                                                              aGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtct                                                                                                                                                                                                                                             >  7:765079/1‑79 (MQ=255)
                                                                                                                                                                                              aGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAg                                                                                                                                                                                            >  5:213542/1‑128 (MQ=255)
                                                                                                                                                                                              aGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAg                                                                                                                                                                                            <  6:213542/128‑1 (MQ=255)
                                                                                                                                                                                               tataagggacaggTCAGGTTCTCTTCTCATTTTCGTTGGTTAACCACGATTATGGGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACGCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGGTTGACCGCATATTATCACCAGGAAAGc                                         <  4:401304/262‑1 (MQ=255)
                                                                                                                                                                                               gATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACg                                                                                                                                                                       >  7:398934/1‑148 (MQ=255)
                                                                                                                                                                                                aTAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGGTTGACCGCATCTTATCACCAGGAAAg                                          >  2:45890/1‑272 (MQ=255)
                                                                                                                                                                                                 tAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGAt                                                                                                                                                                                                                                                                         <  5:790400/48‑1 (MQ=255)
                                                                                                                                                                                                 tAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGAt                                                                                                                                                                                                                                                                         >  6:790400/1‑48 (MQ=255)
                                                                                                                                                                                                 tAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGttt                                                                                                                                                                                              >  7:839783/1‑123 (MQ=255)
                                                                                                                                                                                                 tAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGttt                                                                                                                                                                                              <  8:839783/123‑1 (MQ=255)
                                                                                                                                                                                                    ttCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGa                                                                                                                                                                                           >  2:174586/1‑123 (MQ=255)
                                                                                                                                                                                                    ttCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGc                                                                                                                                                                  <  8:398934/148‑1 (MQ=255)
                                                                                                                                                                                                    ttCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGTTTGACCGCATATTATCACCAGGAAAGCTGATCa                                   >  2:536365/1‑275 (MQ=255)
                                                                                                                                                                                                     tCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACg                                                                                                                                                                       >  7:256571/1‑142 (MQ=255)
                                                                                                                                                                                                     tCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACg                                                                                                                                                                       <  8:256571/142‑1 (MQ=255)
                                                                                                                                                                                                        aTTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGt                                                                                                                                                                                       <  4:174585/123‑1 (MQ=255)
                                                                                                                                                                                                        aTTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtttt                                                                                                                                                              >  6:825223/1‑148 (MQ=255)
                                                                                                                                                                                                         ttGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtct                                                                                                                                                                                                                                             >  6:766468/1‑68 (MQ=255)
                                                                                                                                                                                                         ttGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtct                                                                                                                                                                                                                                             <  5:766468/68‑1 (MQ=255)
                                                                                                                                                                                                         ttGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                                                                                                                                                                   <  5:788475/94‑1 (MQ=255)
                                                                                                                                                                                                         ttGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                                                                                                                                                                   >  6:788475/1‑94 (MQ=255)
                                                                                                                                                                                                             cAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCa                                                                                                                                                                                                                >  5:987309/1‑93 (MQ=255)
                                                                                                                                                                                                             cAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCa                                                                                                                                                                                                                <  6:987309/93‑1 (MQ=255)
                                                                                                                                                                                                             cAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACg                                                                                                                                                                                                         <  7:1034654/100‑1 (MQ=255)
                                                                                                                                                                                                             cAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACg                                                                                                                                                                                                         >  8:1034654/1‑100 (MQ=255)
                                                                                                                                                                                                                gTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                                                                                                                                                                   >  8:498495/1‑87 (MQ=255)
                                                                                                                                                                                                                gTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                                                                                                                                                                   <  7:498495/87‑1 (MQ=255)
                                                                                                                                                                                                                 ttctcttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGtt                                                                                                                                                                                      <  7:1003483/115‑1 (MQ=255)
                                                                                                                                                                                                                 ttctcttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGtt                                                                                                                                                                                      >  8:1003483/1‑115 (MQ=255)
                                                                                                                                                                                                                  tctcttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTgt                                                                                                                                                      >  8:211447/1‑145 (MQ=255)
                                                                                                                                                                                                                  gctcttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTgg                                                                                                                                                      <  7:211447/145‑1 (MQ=255)
                                                                                                                                                                                                                     cttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTggcgggc                                                                                                                                                 >  6:112317/1‑148 (MQ=255)
                                                                                                                                                                                                                      ttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCa                                                                                                   >  2:495654/1‑193 (MQ=255)
                                                                                                                                                                                                                      ttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGGTTGACCGCATATTATCACCAGGAAAGCTGATCATACTCGGTCTGCAacac                  >  2:158219/1‑274 (MQ=255)
                                                                                                                                                                                                                      ttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGGTTGACCGCACATTATCACCAGGa                                             >  2:401305/1‑247 (MQ=255)
                                                                                                                                                                                                                        ctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGaa                                                                                                                                              <  5:825223/148‑1 (MQ=255)
                                                                                                                                                                                                                          cATTTTCGTTTTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTccc                                                                                                                                                                                                                                    <  8:3733/60‑1 (MQ=255)
                                                                                                                                                                                                                          cATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAg                                                                                                                                            >  5:870886/1‑148 (MQ=255)
                                                                                                                                                                                                                          cATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAacac                                                                                                                                   >  3:91095/1‑157 (MQ=255)
                                                                                                                                                                                                                          cATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCAtctt                                                                                               <  4:495653/193‑1 (MQ=255)
                                                                                                                                                                                                                          cATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTATAACTccc                                                                                                                                                                                                                                    >  7:3733/1‑60 (MQ=255)
                                                                                                                                                                                                                            ttttCGTGGTTTAACCGCGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGt                                                                                                                                                                     <  6:840783/121‑1 (MQ=255)
                                                                                                                                                                                                                            ttttCGTGGTTTAACCGCGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGt                                                                                                                                                                     >  5:840783/1‑121 (MQ=255)
                                                                                                                                                                                                                            ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTc                                                                                                                                                                                    <  5:835556/106‑1 (MQ=255)
                                                                                                                                                                                                                            ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTc                                                                                                                                                                                    >  6:835556/1‑106 (MQ=255)
                                                                                                                                                                                                                            ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtt                                                                                                                                                                >  6:489326/1‑126 (MQ=255)
                                                                                                                                                                                                                            ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtt                                                                                                                                                                <  5:489326/126‑1 (MQ=255)
                                                                                                                                                                                                                            ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGga                                                                                                                                          >  7:1016555/1‑148 (MQ=255)
                                                                                                                                                                                                                            ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGga                                                                                                                                          <  5:112317/148‑1 (MQ=255)
                                                                                                                                                                                                                             aatCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtt                                                                                                                                                                <  6:733581/123‑1 (MQ=255)
                                                                                                                                                                                                                             aatCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtt                                                                                                                                                                >  5:733581/3‑125 (MQ=255)
                                                                                                                                                                                                                              ttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACtgat                                                                                                                               <  1:91095/157‑1 (MQ=255)
                                                                                                                                                                                                                              ttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGGTTGACCGCATATTATCACCAGGAAAGCTGATCATACTCGGTCTGCACCACGTCCTTGTc         >  2:336731/1‑275 (MQ=255)
                                                                                                                                                                                                                              ttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGGTTGACCGCATATTATCACCAGGAAAGCTGATCATACTCGGTCTGCAACACGTCCTTGt          >  1:4362/1‑274 (MQ=255)
                                                                                                                                                                                                                              atCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                                                                                                                                                                 <  6:694589/90‑1 (MQ=255)
                                                                                                                                                                                                                              atCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                                                                                                                                                                 <  6:453284/90‑1 (MQ=255)
                                                                                                                                                                                                                              atCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                                                                                                                                                                 >  5:694589/2‑91 (MQ=255)
                                                                                                                                                                                                                              atCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                                                                                                                                                                 >  5:453284/2‑91 (MQ=255)
                                                                                                                                                                                                                                cGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGGTTGACCGCATATTATCACCAGGAAAGCTGATCATACTCGGTCTGCAACACGTCCTTGTCAt       >  4:24915/1‑275 (MQ=255)
                                                                                                                                                                                                                                     tttAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGGTTGACCGCATATTATCACCAGGAAAGCTGATCATACTCGGTCTGCAACACGTCCTTGTCATGTACg  <  2:57297/275‑1 (MQ=255)
                                                                                                                                                                                                                                      ttAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAg                                                                                                                                                                                            <  7:1020508/88‑1 (MQ=255)
                                                                                                                                                                                                                                      ttAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAg                                                                                                                                                                                            >  8:1020508/1‑88 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                
GCATATATTAGGTAGCATGACCGTTTTAATTATAATTTTGCGACTGAGTTCAAATTATTGCCCCTCAATCTATGCCGCTTTCATTATTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGATGAGCGCCATAGATTCCCAACTTCCCTCATCTTCTGGGCAAGACCGCCCAACTGATGAGGTTGACCGCATATTATCACCAGGAAAGCTGATCATACTCGGTCTGCAACACGTCCTTGTCATGTACG  >  NC_000913/3030990‑3031493

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: