Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,603,440 G→A 88.7% intergenic (‑161/+46) yneE ← / ← uxaB bestrophin family putative inner membrane protein/altronate oxidoreductase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,603,4400GA88.7% 969.8 / 96.0 310intergenic (‑161/+46)yneE/uxaBbestrophin family putative inner membrane protein/altronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (14/21);  new base A (141/130);  total (159/151)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 311 total aligned reads displayed.

CGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTG  >  CP009273/1603293‑1603587
                                                                                                                                                   |                                                                                                                                                   
cGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACATACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTaaa                                                                                                                                                   >  2:2087095/1‑150 (MQ=255)
  tGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAAt                                                                                                                                                 >  2:2217250/1‑150 (MQ=255)
  tGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAAt                                                                                                                                                 >  2:1551452/1‑150 (MQ=255)
  tGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAAt                                                                                                                                                 >  2:4231737/1‑150 (MQ=255)
  tGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAAt                                                                                                                                                 >  2:183172/1‑150 (MQ=255)
  tGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAAt                                                                                                                                                 >  2:2667633/1‑150 (MQ=255)
      aaaGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGa                                                                                                                                             <  1:2685796/150‑1 (MQ=255)
       aaGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTCAAATAAGaa                                                                                                                                            >  1:471594/1‑150 (MQ=255)
       aaGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGaa                                                                                                                                            >  1:4818101/1‑150 (MQ=255)
        aGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGaaa                                                                                                                                           >  1:439660/1‑150 (MQ=255)
        aGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGaaa                                                                                                                                           <  1:2667633/150‑1 (MQ=255)
            gCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCg                                                                                                                                       >  2:4645728/1‑150 (MQ=255)
            gCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCg                                                                                                                                       >  1:1449784/1‑150 (MQ=255)
             cccTATTAAACGTGATTACGATCACATTCTCCAGACGTTACTTACAAATTACCTACAAGCTTACAACATGAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCgg                                                                                                                                      <  2:3780361/150‑1 (MQ=255)
             cccTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCgg                                                                                                                                      <  1:1029825/150‑1 (MQ=255)
             cccTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCgg                                                                                                                                      >  1:2968041/1‑150 (MQ=255)
              ccTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGt                                                                                                                                     >  2:388400/1‑150 (MQ=255)
              ccTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGt                                                                                                                                     >  2:2969780/1‑150 (MQ=255)
                  ttAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGaa                                                                                                                                 >  1:3278959/1‑150 (MQ=255)
                     aaCGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAAcc                                                                                                                              <  1:1326570/150‑1 (MQ=255)
                     aaCGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAAcc                                                                                                                              <  1:4322103/150‑1 (MQ=255)
                        gTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCggg                                                                                                                           >  2:3215578/1‑150 (MQ=255)
                          gATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTCATATGGTAAAAGAAGAAACCCGGTCGAAACCGGGtt                                                                                                                         >  2:2959936/1‑150 (MQ=255)
                            ttACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCa                                                                                                                       >  2:2778042/1‑150 (MQ=255)
                              aCGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGa                                                                                                                     <  1:1361882/150‑1 (MQ=255)
                              aCGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGa                                                                                                                     <  1:1710265/150‑1 (MQ=255)
                               cGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGaa                                                                                                                    <  1:3823979/150‑1 (MQ=255)
                                gATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAg                                                                                                                   >  2:3294545/1‑150 (MQ=255)
                                 aTCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGt                                                                                                                  <  1:4171234/150‑1 (MQ=255)
                                     caTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACg                                                                                                              <  1:3776080/150‑1 (MQ=255)
                                     caTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACg                                                                                                              >  2:2575388/1‑150 (MQ=255)
                                      aTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACgg                                                                                                             >  1:4941274/1‑150 (MQ=255)
                                           cAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGtta                                                                                                        >  2:1521171/1‑150 (MQ=255)
                                             aGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGttatt                                                                                                      >  1:2399165/1‑150 (MQ=255)
                                              gACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGttatta                                                                                                     >  1:74163/1‑150 (MQ=255)
                                              gACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGttatta                                                                                                     >  2:1592308/1‑150 (MQ=255)
                                                   tACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGcaca                                                                                                <  1:4037380/150‑1 (MQ=255)
                                                    aCTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACaa                                                                                               <  2:3723410/150‑1 (MQ=255)
                                                     cTTACAAATTACCTACAAGCTTACAACAGTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAAc                                                                                              <  2:5136272/150‑1 (MQ=255)
                                                          aaaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACg                                                                                         >  2:3040322/1‑150 (MQ=255)
                                                          aaaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACg                                                                                         >  2:735349/1‑150 (MQ=255)
                                                          aaaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACg                                                                                         >  2:4245211/1‑150 (MQ=255)
                                                           aaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAAGATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGt                                                                                        >  2:2844306/1‑150 (MQ=255)
                                                                 cTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCt                                                                                  >  1:5079682/1‑150 (MQ=255)
                                                                  tACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCtt                                                                                 <  1:3758179/150‑1 (MQ=255)
                                                                  gaCAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCtt                                                                                 <  2:4632932/149‑1 (MQ=255)
                                                                   aCAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTc                                                                                <  1:586275/150‑1 (MQ=255)
                                                                    cAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTcg                                                                               >  2:1699486/1‑150 (MQ=255)
                                                                        cTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCa                                                                           >  2:1158389/1‑150 (MQ=255)
                                                                         ttACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCAt                                                                          >  1:1581721/1‑150 (MQ=255)
                                                                         ttACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCAt                                                                          >  1:2046101/1‑150 (MQ=255)
                                                                         ttACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCAt                                                                          >  1:2282151/1‑150 (MQ=255)
                                                                           acaacaTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATcc                                                                        <  2:4661828/150‑1 (MQ=255)
                                                                                aTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCtttt                                                                   <  1:3136577/150‑1 (MQ=255)
                                                                                 ttAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCttttt                                                                  <  1:2492932/150‑1 (MQ=255)
                                                                                 ttAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCttttt                                                                  <  2:4718099/150‑1 (MQ=255)
                                                                                  tAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTc                                                                 <  1:3581029/150‑1 (MQ=255)
                                                                                    gCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCga                                                               <  2:2092547/150‑1 (MQ=255)
                                                                                      aGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCgaga                                                             >  1:481169/1‑150 (MQ=255)
                                                                                         cgcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATc                                                          <  1:4257651/150‑1 (MQ=255)
                                                                                         cgcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATc                                                          <  2:1206240/150‑1 (MQ=255)
                                                                                         cgcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATc                                                          >  1:303103/1‑150 (MQ=255)
                                                                                           cTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGc                                                        >  2:1574481/1‑150 (MQ=255)
                                                                                                  tgGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAgg                                                 >  1:1830811/1‑150 (MQ=255)
                                                                                                   ggCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGt                                                >  2:259662/1‑150 (MQ=255)
                                                                                                       ccTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCAtt                                            <  1:3421354/150‑1 (MQ=255)
                                                                                                         tGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAg                                          >  2:2419924/1‑150 (MQ=255)
                                                                                                             aaTGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAAc                                      >  2:4904812/1‑150 (MQ=255)
                                                                                                              aTGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAAcc                                     <  2:41904/150‑1 (MQ=255)
                                                                                                                gATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTg                                   <  2:4992200/150‑1 (MQ=255)
                                                                                                                gATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTg                                   <  1:4351826/150‑1 (MQ=255)
                                                                                                                 aTAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGc                                  >  1:2791552/1‑150 (MQ=255)
                                                                                                                 aTAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGc                                  >  2:1031259/1‑150 (MQ=255)
                                                                                                                  tAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCt                                 <  1:3891784/150‑1 (MQ=255)
                                                                                                                    aaGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCg                               >  2:4457777/1‑150 (MQ=255)
                                                                                                                     aGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCga                              <  2:1412903/150‑1 (MQ=255)
                                                                                                                      ggCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgac                             >  1:544763/1‑150 (MQ=255)
                                                                                                                       tcATATATGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgacc                            <  2:3516499/149‑1 (MQ=255)
                                                                                                                        cATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgacca                           >  2:3474741/1‑150 (MQ=255)
                                                                                                                        cATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgacca                           >  2:3762303/1‑150 (MQ=255)
                                                                                                                         atatAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccag                          >  1:79219/1‑150 (MQ=255)
                                                                                                                         atatAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccag                          <  2:504326/150‑1 (MQ=255)
                                                                                                                             aGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccagacca                      >  1:1128558/1‑150 (MQ=255)
                                                                                                                              ggATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGaccagaccag                     >  2:4409076/1‑150 (MQ=255)
                                                                                                                               gATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAgg                    >  1:866256/1‑150 (MQ=255)
                                                                                                                                    tgtAAATAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACtt               <  2:4818101/150‑1 (MQ=255)
                                                                                                                                      tAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCAGTAGCAACCTGCTCGACCAGACCAGGCACTTGc             >  1:2205730/1‑150 (MQ=255)
                                                                                                                                        aCTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGt           <  2:4182620/150‑1 (MQ=255)
                                                                                                                                        aCTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGt           <  1:2087095/150‑1 (MQ=255)
                                                                                                                                         cTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCATATTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTc          >  2:3384015/1‑150 (MQ=255)
                                                                                                                                         cTAAGATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTc          <  2:4780950/150‑1 (MQ=255)
                                                                                                                                          tAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCa         <  2:39101/150‑1 (MQ=255)
                                                                                                                                          tAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCa         <  1:4224010/150‑1 (MQ=255)
                                                                                                                                           aaTATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAg        >  2:895420/1‑150 (MQ=255)
                                                                                                                                           aaTATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAg        >  2:2133953/1‑150 (MQ=255)
                                                                                                                                              atGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTc     >  1:2719903/1‑150 (MQ=255)
                                                                                                                                              atGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTc     <  2:4115054/150‑1 (MQ=255)
                                                                                                                                               tGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTcc    <  1:1989472/150‑1 (MQ=255)
                                                                                                                                                ggTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCt   <  2:3249377/150‑1 (MQ=255)
                                                                                                                                                 gTGAAATACGACCCCCGGTCGAAACCGGGTTCAGACGTCACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCNg  <  1:3553273/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTG  >  CP009273/1603293‑1603587

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