Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,603,440 G→A 91.9% 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,4400GA91.9% 907.3 / 55.0 285intergenic (‑161/+46)yneE/uxaBbestrophin family putative inner membrane protein/altronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (8/15);  new base A (147/114);  total (156/129)
Fisher's exact test for biased strand distribution p-value = 5.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 288 total aligned reads displayed.

GCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCT  >  CP009273/1603292‑1603589
                                                                                                                                                    |                                                                                                                                                     
  gTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTaaaa                                                                                                                                                    <  1:392823/150‑1 (MQ=255)
     gAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATaa                                                                                                                                                 >  1:1506633/1‑150 (MQ=255)
     gAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATaa                                                                                                                                                 <  1:1651592/150‑1 (MQ=255)
         aGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGaaa                                                                                                                                             >  2:3449846/1‑150 (MQ=255)
         aGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGaaa                                                                                                                                             >  2:3040432/1‑150 (MQ=255)
            cGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAAccc                                                                                                                                          <  1:3011137/150‑1 (MQ=255)
            cGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAAccc                                                                                                                                          <  2:3280790/150‑1 (MQ=255)
              cccTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCgg                                                                                                                                        >  1:2070918/1‑150 (MQ=255)
               ccTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGt                                                                                                                                       <  2:3936956/150‑1 (MQ=255)
                cTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTc                                                                                                                                      <  1:2100642/150‑1 (MQ=255)
                  aTTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGa                                                                                                                                    >  2:3248669/1‑150 (MQ=255)
                   ttAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGaa                                                                                                                                   >  2:660846/1‑150 (MQ=255)
                     aaaCGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAAc                                                                                                                                 <  1:2303628/150‑1 (MQ=255)
                      aaCGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAAcc                                                                                                                                >  1:2526574/1‑150 (MQ=255)
                            aTTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTc                                                                                                                          >  2:3706815/1‑150 (MQ=255)
                             tgaCGATCACATTCTCCCGACGCTACTTACACAGTCCCGACCCGCTGACCCCATGAGCCGGCGCGGCATGGGGCACCTGTTAATGCTAAAGGCATATAGGATGTTGTAACTACTCTGGTGAAAGAAGAAACCCGGTCGACACCGGGTTCa                                                                                                                         <  2:2142434/148‑1 (MQ=255)
                              tACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAg                                                                                                                        >  1:4082678/1‑150 (MQ=255)
                                cGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGaa                                                                                                                      <  2:4444734/150‑1 (MQ=255)
                                 gATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCAGGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCCGAAg                                                                                                                     >  2:3828043/1‑150 (MQ=255)
                                  aTCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGt                                                                                                                    >  2:287630/1‑150 (MQ=255)
                                   tCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTa                                                                                                                   <  2:4475054/150‑1 (MQ=255)
                                     acaTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAAc                                                                                                                 >  1:2836505/1‑150 (MQ=255)
                                      caTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACg                                                                                                                >  2:3063675/1‑150 (MQ=255)
                                       aTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACgg                                                                                                               >  1:4770392/1‑150 (MQ=255)
                                       aTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACgg                                                                                                               >  2:113816/1‑150 (MQ=255)
                                        ttCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGgt                                                                                                              >  2:4340345/1‑150 (MQ=255)
                                         tctcAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGgtg                                                                                                             >  2:2425335/1‑150 (MQ=255)
                                          ctcAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGgtgt                                                                                                            >  2:3830685/1‑150 (MQ=255)
                                          ctcAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGgtgt                                                                                                            <  2:2511286/150‑1 (MQ=255)
                                           tcAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGtt                                                                                                           <  2:2636331/150‑1 (MQ=255)
                                               gACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGttatta                                                                                                       <  1:111110/150‑1 (MQ=255)
                                                  gCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCCGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGca                                                                                                    >  2:2402559/1‑150 (MQ=255)
                                                  gCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCCGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGca                                                                                                    >  2:900666/1‑150 (MQ=255)
                                                   cTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGcac                                                                                                   >  2:1744927/1‑150 (MQ=255)
                                                    tACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGcaca                                                                                                  >  1:1857048/1‑150 (MQ=255)
                                                     aCTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACaa                                                                                                 >  2:3128660/1‑150 (MQ=255)
                                                     aCTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACaa                                                                                                 >  2:1259042/1‑150 (MQ=255)
                                                     aCTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACaa                                                                                                 >  1:115541/1‑150 (MQ=255)
                                                      cTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAAc                                                                                                >  2:4792910/1‑150 (MQ=255)
                                                       ttACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACg                                                                                               <  2:2004668/150‑1 (MQ=255)
                                                         aCAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACACCGGa                                                                                             >  2:1349693/1‑150 (MQ=255)
                                                         aCAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGa                                                                                             >  2:2174347/1‑150 (MQ=255)
                                                           aaaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACg                                                                                           <  2:4076311/150‑1 (MQ=255)
                                                            aaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGt                                                                                          >  2:1460736/1‑150 (MQ=255)
                                                                aCCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAg                                                                                      <  2:342664/150‑1 (MQ=255)
                                                                     cAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTcg                                                                                 >  2:1591051/1‑150 (MQ=255)
                                                                        gCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCgcgc                                                                              <  2:2162114/150‑1 (MQ=255)
                                                                        gCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCgcgc                                                                              <  2:4561381/150‑1 (MQ=255)
                                                                          ttACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCAt                                                                            >  1:2641820/1‑150 (MQ=255)
                                                                          ttACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCAt                                                                            <  1:3374553/150‑1 (MQ=255)
                                                                               acaTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCtt                                                                       <  2:1346526/150‑1 (MQ=255)
                                                                                      cAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCgag                                                                >  1:1620728/1‑150 (MQ=255)
                                                                                       aGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCgaga                                                               >  2:187101/1‑150 (MQ=255)
                                                                                          cgcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATc                                                            <  1:3811805/150‑1 (MQ=255)
                                                                                          cgcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATc                                                            <  1:3529840/150‑1 (MQ=255)
                                                                                            cTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGc                                                          >  1:3624254/1‑150 (MQ=255)
                                                                                             tGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCa                                                         >  1:2682286/1‑150 (MQ=255)
                                                                                                aTGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATcc                                                      >  1:846297/1‑150 (MQ=255)
                                                                                                  gtgGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAg                                                    <  2:1856863/150‑1 (MQ=255)
                                                                                                  gtgGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAg                                                    <  1:4150806/150‑1 (MQ=255)
                                                                                                   tgGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAgg                                                   >  2:1414226/1‑150 (MQ=255)
                                                                                                   tgGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAgg                                                   >  2:4017027/1‑150 (MQ=255)
                                                                                                      cACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCa                                                <  1:300597/150‑1 (MQ=255)
                                                                                                       aCCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCAt                                               <  1:1259042/150‑1 (MQ=255)
                                                                                                       aCCTGTTAATGATAAAGGAATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCAt                                               <  1:687627/150‑1 (MQ=255)
                                                                                                           gTTAATGATAAAGGCATATAGGATGTTGTAACTACTATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGGTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGc                                           >  1:4502106/1‑150 (MQ=255)
                                                                                                             tAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCaa                                         >  1:1383382/1‑150 (MQ=255)
                                                                                                                tGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCt                                      >  2:4432584/1‑150 (MQ=255)
                                                                                                                  aTAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGc                                    >  1:1026225/1‑150 (MQ=255)
                                                                                                                   tAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCt                                   <  1:2916336/150‑1 (MQ=255)
                                                                                                                    aaaGGCATATAGGATGTTGTAACTAATCTGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTc                                  <  1:2764713/150‑1 (MQ=255)
                                                                                                                    aaaGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTc                                  >  1:4565186/1‑150 (MQ=255)
                                                                                                                    aaaGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTc                                  <  1:4690799/150‑1 (MQ=255)
                                                                                                                    aaaGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTc                                  <  1:4219300/150‑1 (MQ=255)
                                                                                                                    aaaGGCATATAGGATGTTGTAACTAATAGGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTc                                  <  2:1941216/150‑1 (MQ=255)
                                                                                                                      aGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCga                                <  1:328768/150‑1 (MQ=255)
                                                                                                                       ggCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgac                               >  2:3815525/1‑150 (MQ=255)
                                                                                                                         cATATAGGATGTTGTAACTACTATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGACGTAACGGTGTTATGAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgacca                             <  1:4634887/150‑1 (MQ=255)
                                                                                                                          atatAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccag                            >  1:476646/1‑150 (MQ=255)
                                                                                                                            ataGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccagac                          <  2:1102791/150‑1 (MQ=255)
                                                                                                                            ataGGATCTTGTAACTAATATGGTAAAATAATAAACCCGGTCGAAACCGCGTTCAGAAGTAAAGTTGTAATTAGCACAACGGACGTACAGCTTCTCTCATCACTTATTCGAGAATCGCATCCACGACATTAGCAACCTGCTCtacaagac                          >  1:760296/1‑150 (MQ=255)
                                                                                                                              aGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccagacca                        >  2:2309208/1‑150 (MQ=255)
                                                                                                                               ggATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGaccagaccag                       >  2:2761066/1‑150 (MQ=255)
                                                                                                                                   gttgtAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGctc                   >  2:4177637/1‑148 (MQ=255)
                                                                                                                                   gttgtAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCAc                   >  2:4177969/1‑150 (MQ=255)
                                                                                                                                     tgtAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACtt                 <  2:2123707/150‑1 (MQ=255)
                                                                                                                                          cTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCCTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTc            <  1:1005867/150‑1 (MQ=255)
                                                                                                                                          cTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTc            >  2:4432497/1‑150 (MQ=255)
                                                                                                                                           tAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCa           >  2:838950/1‑150 (MQ=255)
                                                                                                                                            aaTATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAg          <  1:608867/150‑1 (MQ=255)
                                                                                                                                              tatGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGt        <  1:4439137/150‑1 (MQ=255)
                                                                                                                                              tatGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGt        >  1:3197538/1‑150 (MQ=255)
                                                                                                                                              tatGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGt        >  1:4558419/1‑150 (MQ=255)
                                                                                                                                               atGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTc       <  1:899017/150‑1 (MQ=255)
                                                                                                                                                tGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTcc      >  2:3596074/1‑150 (MQ=255)
                                                                                                                                                tGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTcc      >  1:3681104/1‑150 (MQ=255)
                                                                                                                                                tGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTcc      >  2:1523991/1‑150 (MQ=255)
                                                                                                                                                 ggTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCt     <  2:2101496/150‑1 (MQ=255)
                                                                                                                                                 ggTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCt     <  1:3656263/150‑1 (MQ=255)
                                                                                                                                                 ggTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCt     <  1:4271707/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCT  >  CP009273/1603292‑1603589

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