Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,031,221 T→A 6.6% intergenic (‑290/‑146) ygfT ← / → uacT fused putative oxidoreductase, Fe‑S subunit and nucleotide‑binding subunit/urate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,031,2210TA 6.6% 481.3 / 15.3 158intergenic (‑290/‑146)ygfT/uacTfused putative oxidoreductase, Fe‑S subunit and nucleotide‑binding subunit/urate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (72/75);  new base A (4/7);  total (76/82)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 161 total aligned reads displayed.

ATTATTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTG  >  NC_000913/3031072‑3031366
                                                                                                                                                     |                                                                                                                                                 
attattTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGttt                                                                                                                                                   <  2:1428440/150‑1 (MQ=255)
attattTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGttt                                                                                                                                                   >  1:853942/1‑150 (MQ=255)
 ttattTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTa                                                                                                                                                  <  1:979833/150‑1 (MQ=255)
   attTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAAc                                                                                                                                                >  1:1717722/1‑150 (MQ=255)
    tttCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAAcc                                                                                                                                               >  1:107061/1‑150 (MQ=255)
      tCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCAc                                                                                                                                             >  1:462383/1‑150 (MQ=255)
      tCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCAc                                                                                                                                             >  1:155078/1‑150 (MQ=255)
      tCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCAc                                                                                                                                             >  1:1623281/1‑150 (MQ=255)
          cATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGAtt                                                                                                                                         >  2:757068/1‑150 (MQ=255)
          cATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGAtt                                                                                                                                         >  2:219868/1‑150 (MQ=255)
          cATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGAtt                                                                                                                                         <  1:519533/150‑1 (MQ=255)
          cATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTAAACCACGaaa                                                                                                                                         <  1:2272457/150‑3 (MQ=255)
           aTTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTa                                                                                                                                        <  1:956740/150‑1 (MQ=255)
           aTTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTa                                                                                                                                        <  2:920669/150‑1 (MQ=255)
            ttGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTAAACCACGaaaat                                                                                                                                       <  1:1160946/150‑5 (MQ=255)
             tGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTAAACCACGaaaatg                                                                                                                                      <  2:519995/150‑6 (MQ=255)
              gCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTAtgt                                                                                                                                     >  2:1809632/1‑150 (MQ=255)
              gCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTAtgt                                                                                                                                     <  2:1509970/150‑1 (MQ=255)
              gCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTAAACCACGaaaatga                                                                                                                                     <  1:1789378/150‑7 (MQ=255)
               ccGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTAtgtg                                                                                                                                    >  2:1567562/1‑150 (MQ=255)
                cGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACACATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGa                                                                                                                                   >  1:735970/1‑150 (MQ=255)
                   tCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTg                                                                                                                                <  1:749274/150‑1 (MQ=255)
                    cAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGc                                                                                                                               >  1:265910/1‑150 (MQ=255)
                           gcgcTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCa                                                                                                                        >  2:584892/1‑150 (MQ=255)
                              cTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAAtt                                                                                                                     >  1:1890271/1‑150 (MQ=255)
                                aTCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTcc                                                                                                                   <  2:1143275/150‑1 (MQ=255)
                                aTCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGAAAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTcc                                                                                                                   <  1:2389461/150‑1 (MQ=255)
                                    aGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtc                                                                                                               <  1:2228034/150‑1 (MQ=255)
                                    aGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtc                                                                                                               <  1:584892/150‑1 (MQ=255)
                                    aGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtc                                                                                                               <  2:1260158/150‑1 (MQ=255)
                                    aGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtc                                                                                                               <  2:1678451/150‑1 (MQ=255)
                                      aTAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtctc                                                                                                             <  2:112362/150‑1 (MQ=255)
                                      aTAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtctc                                                                                                             <  1:132983/150‑1 (MQ=255)
                                      aTAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtctc                                                                                                             >  1:1908985/1‑150 (MQ=255)
                                        aaaaaCGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTa                                                                                                           <  2:1063600/150‑1 (MQ=255)
                                         aaaaCGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTaa                                                                                                          >  2:1332147/1‑150 (MQ=255)
                                         aaaaCGATTAAATGCGAGAGTGCGACATGCCAGAATGATAAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTaa                                                                                                          >  2:1144636/1‑150 (MQ=255)
                                           aaCGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACt                                                                                                        >  1:431160/1‑150 (MQ=255)
                                             cGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTcc                                                                                                      >  2:1865975/1‑150 (MQ=255)
                                             cGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTcc                                                                                                      >  1:215334/1‑150 (MQ=255)
                                                 tAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCtt                                                                                                  <  1:1050351/150‑1 (MQ=255)
                                                 tAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCtt                                                                                                  <  1:2015635/150‑1 (MQ=255)
                                                      gCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCg                                                                                             <  2:941825/150‑1 (MQ=255)
                                                      gCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCg                                                                                             <  2:155078/150‑1 (MQ=255)
                                                      gCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACACATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCg                                                                                             <  2:735970/150‑1 (MQ=255)
                                                       cGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGc                                                                                            >  2:859061/1‑150 (MQ=255)
                                                       cGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGc                                                                                            >  1:2023775/1‑150 (MQ=255)
                                                       cGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGc                                                                                            >  2:1402542/1‑150 (MQ=255)
                                                          gagTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCaaa                                                                                         <  1:1243526/150‑1 (MQ=255)
                                                               cGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                                    <  1:219868/150‑1 (MQ=255)
                                                                gACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTgg                                                                                   <  2:1902955/150‑1 (MQ=255)
                                                                 aCATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGc                                                                                  <  2:2064337/150‑1 (MQ=255)
                                                                  cATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCa                                                                                 <  2:1446543/150‑1 (MQ=255)
                                                                     gCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTc                                                                              >  1:1980342/1‑150 (MQ=255)
                                                                      ccAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCACATTCCTCGCAAAAACTGGCACTcc                                                                             >  1:975518/1‑150 (MQ=255)
                                                                           aTGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAg                                                                        <  1:2344490/150‑1 (MQ=255)
                                                                              aTTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTAGCAAAAACTGGCACTCCACGAGCAt                                                                     >  2:767360/1‑150 (MQ=255)
                                                                               ttAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtg                                                                    >  1:685202/1‑150 (MQ=255)
                                                                                tAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgt                                                                   >  1:840827/1‑150 (MQ=255)
                                                                                 aaTGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                                  >  1:1540800/1‑150 (MQ=255)
                                                                                 aaTGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                                  >  2:2253205/1‑150 (MQ=255)
                                                                                 aaTGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                                  <  1:1822317/150‑1 (MQ=255)
                                                                                 aaTGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAAACAAGATAATGTGATTGCCTCCGCAATTCCAGTCTCAAACTCCCCTTCCTCGCAAAAACTCGCACTCCACGAGCAagtg                                                                  >  1:55437/1‑150 (MQ=255)
                                                                                  aTGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATgtgt                                                                 >  1:267029/1‑150 (MQ=255)
                                                                                     aaaTATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTa                                                              <  2:2485523/150‑1 (MQ=255)
                                                                                      aaTATCAATATGTGAAAATTACAGATATTTCCATTGTCAGGTGCTGTTATCATTTTCGTTGTTTAACCACGATTATGTGATTTAATCCGAAATTACTGTCTGTAACTCCGCTTCCTCGCAAAAACTGTCACTCCACGAGCATGTGTTTAg                                                             <  2:1717722/150‑1 (MQ=255)
                                                                                         atCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAAc                                                                                      >  1:1447103/1‑122 (MQ=255)
                                                                                            aaTATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCATTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGtt                                                       >  2:1722851/1‑150 (MQ=255)
                                                                                               atGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCa                                                    <  1:614575/150‑1 (MQ=255)
                                                                                                tGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCAt                                                   <  2:1625685/150‑1 (MQ=255)
                                                                                                tGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCAt                                                   <  2:518957/150‑1 (MQ=255)
                                                                                                 gagaAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCAtt                                                  >  1:2355094/1‑150 (MQ=255)
                                                                                                 gagaAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCAtt                                                  <  2:462383/150‑1 (MQ=255)
                                                                                                 gagaAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCAtt                                                  <  2:2371302/150‑1 (MQ=255)
                                                                                                    aaaaTTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGGTTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAAc                                               <  2:2308076/150‑1 (MQ=255)
                                                                                                     aaaTTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACg                                              >  2:294943/1‑150 (MQ=255)
                                                                                                     aaaTTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACg                                              <  1:634317/150‑1 (MQ=255)
                                                                                                      aaTTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGt                                             >  1:412043/1‑150 (MQ=255)
                                                                                                      aaTTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACAACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGt                                             <  2:2080838/150‑1 (MQ=255)
                                                                                                             gATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGt                                      <  1:1319751/150‑1 (MQ=255)
                                                                                                             gATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGt                                      >  1:1861221/1‑150 (MQ=255)
                                                                                                               tAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTTGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTACTTTTACTCGATAATACTTTCACTCCAAGATAATGTGTTTATAAAGTTTCATTTACTTTAAGGGttt                                    <  2:1950725/150‑2 (MQ=255)
                                                                                                                   tCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCttt                                >  1:1409256/1‑150 (MQ=255)
                                                                                                                     cATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCttttt                              <  1:57579/150‑1 (MQ=255)
                                                                                                                     cATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCttttt                              >  2:737467/1‑150 (MQ=255)
                                                                                                                       ttGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTAc                            <  2:1623281/150‑1 (MQ=255)
                                                                                                                         gTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTActc                          >  2:213058/1‑150 (MQ=255)
                                                                                                                           cAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTg                        >  1:1444938/1‑150 (MQ=255)
                                                                                                                            aGGTTCTCTTCTCATTTTTGTTTTTTAACCACGATTATGTGATTGCCTACGCAATTACTGTCTCTAACTCCCCTTCCTAGCAAAAACTGGAACTCCTCGAGCATGTGTTTAGACAGTTTCATTAACTTAAACGGTTTCTTTTTTCTGTgg                       <  2:712573/150‑1 (MQ=255)
                                                                                                                             ggTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTggc                      <  1:1763655/150‑1 (MQ=255)
                                                                                                                              gTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTggcg                     <  1:2150424/150‑1 (MQ=255)
                                                                                                                              gTTATCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCAATAACGTAAACGGTTGCTTTTTACTCTggcg                     >  2:600017/1‑150 (MQ=255)
                                                                                                                                      ctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAg             <  1:1972582/150‑1 (MQ=255)
                                                                                                                                        cataatCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGga           >  1:1594406/6‑150 (MQ=255)
                                                                                                                                         aTTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGgag          <  2:100336/150‑1 (MQ=255)
                                                                                                                                          ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGgaga         >  1:318316/1‑150 (MQ=255)
                                                                                                                                          taatCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGgaga         >  1:1127727/4‑150 (MQ=255)
                                                                                                                                           tttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGaa        >  2:286894/1‑150 (MQ=255)
                                                                                                                                             tCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAac      <  2:1762624/150‑1 (MQ=255)
                                                                                                                                             tCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGATTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAac      >  2:979771/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ATTATTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTG  >  NC_000913/3031072‑3031366

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