Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,222 T→G 6.3% intergenic (+17/+35) yfeD → / ← gltX DUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,519,2220TG 6.3% 310.6 / 8.7 106intergenic (+17/+35)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (51/47);  new base G (4/3);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 111 total aligned reads displayed.

CCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGC  >  NC_000913/2519074‑2519370
                                                                                                                                                    |                                                                                                                                                    
ccGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCATGata                                                                                                                                                     >  2:426686/1‑150 (MQ=255)
ccGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGata                                                                                                                                                     >  1:64095/1‑150 (MQ=255)
 cGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGatat                                                                                                                                                    <  1:1065228/150‑1 (MQ=255)
 cGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGatat                                                                                                                                                    <  1:12092/150‑1 (MQ=255)
       aTCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtctc                                                                                                                                              >  1:931687/1‑150 (MQ=255)
        tCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACTGCAGGATATTATCTcc                                                                                                                                             >  1:554029/1‑150 (MQ=255)
          gACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                                                                           >  2:381550/1‑150 (MQ=255)
            cAAATCATCAACGCGAGAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                                                                         <  2:104622/150‑1 (MQ=255)
            aaaaTCATCAACGCGATAGATAATATGACAGAAATTGAGCAGTAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCACGTCTCGTTATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                                                                         <  1:996834/149‑1 (MQ=255)
             aaaTCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                                                                        <  1:852280/150‑1 (MQ=255)
               atcatcAACGCGATAGATAATATGACATCCATTGACCAGTAAAATGTCGCGCATTTTTTGTAACTTGAAGTAATTCGCAATTTCCTGGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                                                                      <  2:1453005/150‑1 (MQ=255)
               atcatcAACGCGATAGATAATATGACAGACATTGAGCAGTAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                                                                      <  2:726271/150‑1 (MQ=255)
                tcatcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGttt                                                                                                                                     >  2:1484860/1‑150 (MQ=255)
                 aatcaACGCTATATATAATATGACAGACATTTAGCAGGAAAAAGTCGCGAAATTTTTGTCACGTGAAGTAATTCGCAATTTCCTCGCCCGTCTCTATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTa                                                                                                                                    <  2:1906612/149‑1 (MQ=255)
                  atcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttt                                                                                                                                   <  2:1928392/150‑3 (MQ=255)
                    caACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCt                                                                                                                                 >  2:1553583/1‑150 (MQ=255)
                    caACGCGATAGATAATATGACAGACATTGAGCAAGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTAGATATTGAAGAATCATCATAAATAAAAAACGGCAGTATATTATCTCATGACGTAAATCt                                                                                                                                 >  1:202733/1‑150 (MQ=255)
                    caACGCGATAGATAAAATGAAAGACATTGAGCAGGACAAAGTCGCGCAATTTTTGTCACTTGAAGTAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttttta                                                                                                                                 <  1:1674136/150‑5 (MQ=255)
                    caACGAGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAAAGGCAGGATATTATCTCCTGGCGTTTATCt                                                                                                                                 >  2:628196/1‑150 (MQ=255)
                      aCGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttt                                                                                                                               >  1:443832/1‑150 (MQ=255)
                         cgATAGATAATATGACAGACATTGATCAGTAAAAAGTCTCTCAATTTTTGTCACTTGAATGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTa                                                                                                                            <  1:1423418/150‑1 (MQ=255)
                         cgATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGAGATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTa                                                                                                                            >  2:1818146/1‑150 (MQ=255)
                          gatagataTTATGACAGACATTGAGCCGGAAAAAGTCGCGCAATTTTTGACACTTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGGTATTATCTCCTGCCGTTTATCTTTTTac                                                                                                                           <  1:358639/150‑1 (MQ=255)
                           atagataATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                                                                          >  2:1810243/1‑150 (MQ=255)
                              gataATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTTTCTCCTGCCGTTTATCTTTTTACACGc                                                                                                                       >  1:682814/1‑150 (MQ=255)
                              gataATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                                                                       <  1:1846343/150‑1 (MQ=255)
                               ataataTGACAGCCATTGAGCAGCAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCt                                                                                                                      <  2:280071/150‑1 (MQ=255)
                               ataataTGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCt                                                                                                                      <  2:1223293/150‑1 (MQ=255)
                                taataTGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTa                                                                                                                     >  2:845957/1‑150 (MQ=255)
                                  atatGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAt                                                                                                                   >  1:1125238/1‑150 (MQ=255)
                                      tacagacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAcc                                                                                                               <  1:1065011/149‑2 (MQ=255)
                                        cagacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                                                             <  1:794860/150‑1 (MQ=255)
                                        cagacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                                                             >  2:871459/1‑150 (MQ=255)
                                         agacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                                                                            <  1:1676826/150‑1 (MQ=255)
                                              ttGAGCAGGAAAAATTCTCGCAATTTTTGTCACGTGAAGTAATTCGCAATTTCCTCGCCCGTCTCGATCTGGCCGACTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAttt                                                                                                       <  1:1541403/150‑1 (MQ=255)
                                              ttGAGCAGGAAAAAGTAGCGCAATTTTTGTCACGTGAAGGAATTAGCAATTTCCTCGCCCTTCTCGAGATTGACGAATCAGCATAAATAAAAAAACTCAGGATATTATCTCCTGCCGTTAATATTTTTACACCCTAATTACTGCTGAttt                                                                                                       >  2:887522/1‑150 (MQ=255)
                                                taccAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAtttta                                                                                                     >  2:898195/4‑149 (MQ=255)
                                                gAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTc                                                                                                     >  1:1999228/1‑150 (MQ=255)
                                                  gCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATCTCCTCGCTCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCAGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgc                                                                                                   >  1:1466843/1‑150 (MQ=255)
                                                     ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACTGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                                                                                >  2:813242/1‑150 (MQ=255)
                                                       aaaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                                                                              <  1:648624/150‑1 (MQ=255)
                                                         aaaGTTGCGCAATTTTTGTCACGTGAATGAATTCTCAATTTCATCTCCCGTCTCGATATTGACGAATCAGCATAATTAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAga                                                                                            <  1:183106/150‑2 (MQ=255)
                                                         aaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                                                                            <  2:1119414/150‑1 (MQ=255)
                                                           aGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCaa                                                                                          <  2:1470368/150‑1 (MQ=255)
                                                             tCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCGCTTCTCGATATTGAAGAATCATCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATa                                                                                        <  1:1694147/150‑1 (MQ=255)
                                                              cgcgCATTTTTTGTCAAGTGAAGTAATTCGCAATTTCCTCGCCCGTCTCGATATTTACGAATCATCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTAGCGTTCAGCAATaa                                                                                       <  2:485555/150‑1 (MQ=255)
                                                              cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACTAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                                                                       <  1:1703882/150‑1 (MQ=255)
                                                                cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaaa                                                                                     >  1:950816/1‑150 (MQ=255)
                                                                  cAATTTTTGTCACGGGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTATTACAAGCTAATTAATGCTGATTTTCGCGTTCAGCAATAAAATc                                                                                   >  1:1115432/1‑150 (MQ=255)
                                                                   aaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                                                                  <  2:1167436/150‑1 (MQ=255)
                                                                    aTTTTTGTCACGTGATGGAATTCGCAATATACACGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGTCAGTAAATTAACTACTGCAGATTATCTTTTTCCACGATAATTAGTGCTGATTTTCGCGAGCATCAATAAAATcaa                                                                                 >  2:377100/1‑148 (MQ=255)
                                                                          gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCtt                                                                           >  2:1746534/1‑150 (MQ=255)
                                                                            cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATTTTGACGAATCAGCATAAATAAAAAACTGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                                                                         <  1:21819/150‑1 (MQ=255)
                                                                            cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGAGTTAAGCAATAAAATCCAGCGCtttt                                                                         >  1:566458/1‑149 (MQ=255)
                                                                              cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                                                                       <  1:1601363/150‑1 (MQ=255)
                                                                                 gAAGTAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                                                                    <  2:2012047/150‑1 (MQ=255)
                                                                                 gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                                                                    <  2:1492655/150‑1 (MQ=255)
                                                                                       aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc                                                              <  2:520933/150‑1 (MQ=255)
                                                                                       aTTCGCAATTTCCTCGCCAGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc                                                              >  2:1958214/1‑150 (MQ=255)
                                                                                          cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                                                           <  1:1039308/150‑1 (MQ=255)
                                                                                          cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                                                           <  2:1997681/150‑1 (MQ=255)
                                                                                            cAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAAAGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc                                                         >  1:106661/1‑150 (MQ=255)
                                                                                              aTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCTCTTTGTTGATACGCTCGATACTg                                                       >  2:1846357/1‑150 (MQ=255)
                                                                                                 tCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTAAACGCTAATTACTGCTGATTTTCGCGTTAAGCAATAAAATCCAGCGCTTTGTTGATACGCACGATACTGCgg                                                    >  1:637253/1‑150 (MQ=255)
                                                                                                   cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATACTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGt                                                  >  2:619169/1‑150 (MQ=255)
                                                                                                   cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATACTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGt                                                  >  2:1571055/1‑150 (MQ=255)
                                                                                                        ccGTCTCGATATTGACGAATCATCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAc                                             <  1:108951/150‑1 (MQ=255)
                                                                                                        ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAc                                             <  1:460397/150‑1 (MQ=255)
                                                                                                         cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAcc                                            <  2:1464962/150‑1 (MQ=255)
                                                                                                         cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAcc                                            <  1:1378926/150‑1 (MQ=255)
                                                                                                         cGTCTCGATATTGACGAAACAGCATAAATAAAAAACGGCATGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAcc                                            <  1:1240326/150‑1 (MQ=255)
                                                                                                              cGATATTGACGAATCATCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTg                                       <  1:845957/150‑1 (MQ=255)
                                                                                                                atatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCg                                     <  1:381550/150‑1 (MQ=255)
                                                                                                                   ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGa                                  >  1:1923014/1‑150 (MQ=255)
                                                                                                                    tGACGAATCAGCATAAATAAAAAACGGCAGGATATTTTCTCCTGCCGTTTATCTTGTTACACGCTAATTACTGCTGATTTTCGCGTTGAGCAATAACATCCAGAGCTTTGTTGATACGCTCGATCCTGCGGGTCTTACGAATGGCGTGaa                                 >  2:1228944/1‑150 (MQ=255)
                                                                                                                     gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAAc                                <  2:1126804/150‑1 (MQ=255)
                                                                                                                        gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGt                             >  1:428537/1‑150 (MQ=255)
                                                                                                                         aaTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTa                            <  1:1938510/150‑1 (MQ=255)
                                                                                                                             ggCATAAATAAAAAACGGCAGGATATTATATCCTGCCGTTTATATTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGTGTCTTACCAATTGCGTGAACGGTAACAt                        >  1:1347359/2‑150 (MQ=255)
                                                                                                                             aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCATGAACGGTAACAt                        <  1:515577/150‑1 (MQ=255)
                                                                                                                               cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCATTGTTGATACGCTCTATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATcc                      >  2:1264700/1‑150 (MQ=255)
                                                                                                                                ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAGTTGCGTGAACGGTAACATCCa                     >  2:1468259/1‑150 (MQ=255)
                                                                                                                                  aaataaaAAACGGCAGGTTATTATCTCCTGCCGTTTATCTTTTTACACTCTAATAACTGCTGATTATCGCGTTCAGCAATAAAATCCAGCGCTTTGATGATACGCTCGAAACTGCGGGTCTTACCAATTGAGTGAACGGTAACATCCatt                   >  1:1430992/1‑148 (MQ=255)
                                                                                                                                   aataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTTTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTg                  >  2:370234/1‑150 (MQ=255)
                                                                                                                                   aaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTg                  >  2:1158084/7‑150 (MQ=255)
                                                                                                                                    ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGc                 >  2:552211/1‑150 (MQ=255)
                                                                                                                                    ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGc                 >  2:1144398/1‑150 (MQ=255)
                                                                                                                                     taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTAGCGTTCAGCAATAAAATCCAGCGCTTTCTTGATACGATCGAAACTGCGTGTCTTACCAATTTCGTGAACGGTAACATCCAGTGCt                >  1:793739/1‑150 (MQ=255)
                                                                                                                                      aaaaatCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTg               <  1:110345/150‑1 (MQ=255)
                                                                                                                                      aaaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAGTGCTGATTTTGGCGTTCAGCAATAAAATCCATCGGATTGTTAATACCCTCGATACTGCTGTTCTTAAAAATTGCGTGAACTGTAAGATCCAGAGCTg               >  2:1246108/1‑150 (MQ=255)
                                                                                                                                      aaaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTg               >  2:1127517/1‑150 (MQ=255)
                                                                                                                                       aaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTgg              <  1:1746534/150‑1 (MQ=255)
                                                                                                                                        aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGga             >  1:1905524/1‑150 (MQ=255)
                                                                                                                                         aaaCTGCAGGATATTATCTACTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAATTCCAGCGCTTTTTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGgag            <  1:628196/150‑1 (MQ=255)
                                                                                                                                          aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAACTCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGgaga           >  1:1881517/1‑150 (MQ=255)
                                                                                                                                          aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGgaga           >  1:930145/1‑150 (MQ=255)
                                                                                                                                          aaCGGCAGGAGATAATATCCTGCCGTTTATCATTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAAAAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAACAATTGCGTGAACGGTAACATCCAGAGCTGgaga           >  1:1738787/1‑150 (MQ=255)
                                                                                                                                             ggCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAAACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTg        >  2:334025/1‑150 (MQ=255)
                                                                                                                                              gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGc       >  2:89963/1‑150 (MQ=255)
                                                                                                                                                 ggATATTATCTCCTGCCGTTTATCTTTTTACCCGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGcccc    >  1:513183/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGC  >  NC_000913/2519074‑2519370

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