Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,228 C→A 6.8% intergenic (+23/+29) 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,2280CA 6.8% 307.7 / 10.5 111intergenic (+23/+29)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (52/51);  new base A (5/3);  total (57/54)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

Only 100 of 113 total aligned reads displayed.

ATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCACCGGTT  >  NC_000913/2519081‑2519377
                                                                                                                                                   |                                                                                                                                                     
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)
        atcctcAACGCGATAGATAATATGACAGACATTGAGCTGGAAAAAGTCGCGCAATTTTTGTCACGTGAATGAATTCGCAATTTCCTCGACCTTCTCGATATTGACGAATCAGCATAAATCAAACACGGCAGGAGATAATATCCTGCCGtt                                                                                                                                             <  1:169405/150‑1 (MQ=255)
        atcatcAACGCGATAGATAATATGACATCCATTGACCAGTAAAATGTCGCGCATTTTTTGTAACTTGAAGTAATTCGCAATTTCCTGGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                                                                             <  2:1453005/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)
             caACGCGATAGATAATATGACAGACATTGAGCAAGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTAGATATTGAAGAATCATCATAAATAAAAAACGGCAGTATATTATCTCATGACGTAAATCt                                                                                                                                        >  1:202733/1‑150 (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)
                    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                                                                                                                    >  2:871459/1‑150 (MQ=255)
                                 cagacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                                                                    <  1:794860/150‑1 (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)
                                         gAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTc                                                                                                            >  1:499935/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)
                                                    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)
                                                                 ttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATACAgcgc                                                                                    >  1:1927835/1‑150 (MQ=255)
                                                                   gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCtt                                                                                  >  2:1746534/1‑150 (MQ=255)
                                                                     cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATTTTGACGAATCAGCATAAATAAAAAACTGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                                                                                <  1:21819/150‑1 (MQ=255)
                                                                       cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                                                                              <  1:1601363/150‑1 (MQ=255)
                                                                       cGTGAAGGAATTCGCAATTTCCTAGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                                                                              >  2:110345/1‑150 (MQ=255)
                                                                          gAAGTAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                                                                           <  2:2012047/150‑1 (MQ=255)
                                                                          gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGaa                                                                           >  2:1588979/1‑149 (MQ=255)
                                                                          gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                                                                           <  2:1492655/150‑1 (MQ=255)
                                                                                aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc                                                                     <  2:520933/150‑1 (MQ=255)
                                                                                aTTCGCAATTTCCTCGCCAGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc                                                                     >  2:1958214/1‑150 (MQ=255)
                                                                                   cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTTACTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                                                                  >  1:1337753/1‑150 (MQ=255)
                                                                                   cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                                                                  <  2:1997681/150‑1 (MQ=255)
                                                                                   cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                                                                  <  1:1039308/150‑1 (MQ=255)
                                                                                   cGCAATTTCCTCGCCCGTCTCGATATTGAATAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                                                                  <  1:757979/150‑1 (MQ=255)
                                                                                     cAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAAAGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc                                                                >  1:106661/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)
                                                                                                 ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCACGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAc                                                    <  2:259187/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)
                                                                                                         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)
                                                                                                               aCGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACg                                      <  2:554029/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)
                                                                                                                            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)
                                                                                                                                   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)
                                                                                                                                              attatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCAc       <  2:1798256/150‑1 (MQ=255)
                                                                                                                                               ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCAcc      <  2:1602565/150‑1 (MQ=255)
                                                                                                                                                tatCTCCTGCCGATTATCTTTTTACACGCAAATTACTGCAGATTTTCGCCTTCAGCAATAAAATCCAGCGCTTAGTTGATACGCTCCAAACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTCCTGGAGACTGCCCCGCAccc     >  2:1447459/1‑149 (MQ=255)
                                                                                                                                                 atCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGACTTACCAATTCCGTGAACGCTAACATCCAGTGCTGGAGACTGCCAAGCACCgg    >  2:835202/1‑150 (MQ=255)
                                                                                                                                                   ctcCTGCCGTTTTTCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCACCGGtt  >  1:743889/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCACCGGTT  >  NC_000913/2519081‑2519377

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: