Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,603,440 G→A 86.5% 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,4400GA86.5% 779.4 / 96.2 275intergenic (‑161/+46)yneE/uxaBbestrophin family putative inner membrane protein/altronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (10/27);  new base A (109/128);  total (120/155)
Fisher's exact test for biased strand distribution p-value = 3.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 280 total aligned reads displayed.

GTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCT  >  CP009273/1603294‑1603586
                                                                                                                                                  |                                                                                                                                                  
gTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTaaaa                                                                                                                                                 <  2:4199365/150‑1 (MQ=255)
gTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTaaaa                                                                                                                                                 <  2:1952411/150‑1 (MQ=255)
  ggAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATa                                                                                                                                               <  2:1885053/150‑1 (MQ=255)
   gAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATaa                                                                                                                                              <  1:3944636/150‑1 (MQ=255)
     aaaGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGa                                                                                                                                            <  1:113728/150‑1 (MQ=255)
     aaaGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGa                                                                                                                                            >  2:321575/1‑150 (MQ=255)
      aaGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTTGAAGCTTAAAACATTAGCAGGCGCTGCATGTGGCAACTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGaa                                                                                                                                           <  1:3693347/150‑1 (MQ=255)
      aaGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGaa                                                                                                                                           <  1:173970/150‑1 (MQ=255)
        gTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGTATGTTGTAACTAATATGGTAAAATAAGAAAc                                                                                                                                         <  2:1932568/150‑1 (MQ=255)
        gTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAAc                                                                                                                                         <  1:2784584/150‑1 (MQ=255)
         tCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAAcc                                                                                                                                        <  1:1313358/150‑1 (MQ=255)
         tCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAAcc                                                                                                                                        >  1:633078/1‑150 (MQ=255)
           gCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCg                                                                                                                                      <  2:625888/150‑1 (MQ=255)
             ccTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGt                                                                                                                                    <  2:1443585/150‑1 (MQ=255)
               tATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCg                                                                                                                                  <  1:4419068/150‑1 (MQ=255)
                aTTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGa                                                                                                                                 <  2:2243978/150‑1 (MQ=255)
                 ttAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGaa                                                                                                                                <  2:3088187/150‑1 (MQ=255)
                 ttAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGaa                                                                                                                                <  2:3877694/150‑1 (MQ=255)
                   aaaCGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAAc                                                                                                                              <  1:2436149/150‑1 (MQ=255)
                   aaaCGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAAc                                                                                                                              <  1:304436/150‑1 (MQ=255)
                        tGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGt                                                                                                                         <  2:3639857/150‑1 (MQ=255)
                        tGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGt                                                                                                                         <  1:7942/150‑1 (MQ=255)
                          aTTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTc                                                                                                                       >  1:3469337/1‑150 (MQ=255)
                           ttACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCa                                                                                                                      >  1:3963197/1‑150 (MQ=255)
                            tACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAg                                                                                                                     >  2:123959/1‑150 (MQ=255)
                              cGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGaa                                                                                                                   >  2:863570/1‑150 (MQ=255)
                                aTCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGt                                                                                                                 >  1:3320261/1‑150 (MQ=255)
                                aTCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGt                                                                                                                 <  1:253416/150‑1 (MQ=255)
                                  cacaTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTaa                                                                                                               <  1:3738517/150‑1 (MQ=255)
                                     aTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACgg                                                                                                            >  2:2567134/1‑150 (MQ=255)
                                      ttCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGgt                                                                                                           <  2:3135152/150‑1 (MQ=255)
                                      ttCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGgt                                                                                                           <  1:106015/150‑1 (MQ=255)
                                       tctcAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGgtg                                                                                                          >  1:2074348/1‑150 (MQ=255)
                                        ctcAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGgtgt                                                                                                         >  1:1386796/1‑150 (MQ=255)
                                         tcAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGtt                                                                                                        >  2:3229124/1‑150 (MQ=255)
                                           aaGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGttat                                                                                                      >  1:4003982/1‑150 (MQ=255)
                                            aGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGttatt                                                                                                     <  1:2251051/150‑1 (MQ=255)
                                             gACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGttatta                                                                                                    <  1:4053697/150‑1 (MQ=255)
                                              aCGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACAGGGTTCAGAAGTAACGGTGTTATTAg                                                                                                   >  2:1976979/1‑150 (MQ=255)
                                                gCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGca                                                                                                 >  2:3308884/1‑150 (MQ=255)
                                                gCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGca                                                                                                 <  2:263697/150‑1 (MQ=255)
                                                   aCTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACaa                                                                                              <  1:1903294/150‑1 (MQ=255)
                                                    cTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAAc                                                                                             <  1:3893796/150‑1 (MQ=255)
                                                    cTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAAc                                                                                             >  2:81173/1‑150 (MQ=255)
                                                       aCAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGa                                                                                          >  2:743277/1‑150 (MQ=255)
                                                        cAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGAc                                                                                         >  2:3666026/1‑150 (MQ=255)
                                                         aaaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGGAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACg                                                                                        >  2:2995173/1‑150 (MQ=255)
                                                         aaaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACg                                                                                        >  1:3264698/1‑150 (MQ=255)
                                                         aaaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACg                                                                                        >  2:1646829/1‑150 (MQ=255)
                                                         aaaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACg                                                                                        >  1:986042/1‑150 (MQ=255)
                                                                 tACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCtt                                                                                >  2:2156106/1‑150 (MQ=255)
                                                                 tACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCtt                                                                                <  1:2374004/150‑1 (MQ=255)
                                                                  aCAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTc                                                                               >  1:223296/1‑150 (MQ=255)
                                                                  aCAAGCTTACAACATGAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTc                                                                               <  2:1838567/150‑1 (MQ=255)
                                                                   cAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTcg                                                                              <  1:3966632/150‑1 (MQ=255)
                                                                     aGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTcgcg                                                                            >  1:1984770/1‑150 (MQ=255)
                                                                       cTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCa                                                                          >  1:3175795/1‑150 (MQ=255)
                                                                       cTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCa                                                                          >  1:3175629/1‑150 (MQ=255)
                                                                       cTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCa                                                                          <  2:2798721/150‑1 (MQ=255)
                                                                       cTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCa                                                                          >  1:341832/1‑150 (MQ=255)
                                                                       cTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCa                                                                          >  1:2644882/1‑150 (MQ=255)
                                                                            aacaTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCt                                                                     <  1:3885402/150‑1 (MQ=255)
                                                                              caTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATGCGCACAACGGCCGTACAGCTTCGCGCATCCCttt                                                                   >  2:4410312/1‑150 (MQ=255)
                                                                               aTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCtttt                                                                  >  1:2021635/1‑150 (MQ=255)
                                                                               aTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCtttt                                                                  >  1:4269109/1‑150 (MQ=255)
                                                                                ttAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGGTATTAGCACAACGGACGTACAGCTTCGCGCATCCCttttt                                                                 >  2:3079832/1‑150 (MQ=255)
                                                                                  aGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCg                                                               >  1:1126865/1‑150 (MQ=255)
                                                                                     aGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTGTCgaga                                                            >  1:922650/1‑150 (MQ=255)
                                                                                      ggCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGaa                                                           >  2:2578032/1‑150 (MQ=255)
                                                                                       gcgcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAAt                                                          <  2:3443298/150‑1 (MQ=255)
                                                                                        cgcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATc                                                         >  1:753471/1‑150 (MQ=255)
                                                                                         gcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCg                                                        <  1:542282/150‑1 (MQ=255)
                                                                                         gcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCg                                                        >  2:1395812/1‑150 (MQ=255)
                                                                                         gcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCg                                                        >  2:3470999/1‑150 (MQ=255)
                                                                                            gCATGTGGCACCTGTTAATGATAAAGGCATATAGGCTGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCAt                                                     <  2:1570867/150‑1 (MQ=255)
                                                                                            gCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCAt                                                     <  2:523071/150‑1 (MQ=255)
                                                                                             cATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATc                                                    <  2:3031445/150‑1 (MQ=255)
                                                                                               tgtgGCACCTGTTAATGATAAAGGCATATAGGATTTTGTAACTAATATGGTAAAATAAGAAACCAGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCa                                                  <  2:2522906/150‑1 (MQ=255)
                                                                                               tgtgGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCa                                                  <  2:2666669/150‑1 (MQ=255)
                                                                                                           tAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCaa                                      >  1:499211/1‑150 (MQ=255)
                                                                                                             aTGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAAcc                                    >  2:1676299/1‑150 (MQ=255)
                                                                                                              tGATAAAGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCt                                   >  2:3160518/1‑150 (MQ=255)
                                                                                                                 tAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCt                                <  1:4147926/150‑1 (MQ=255)
                                                                                                                  aaaGGCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTc                               <  2:1386796/150‑1 (MQ=255)
                                                                                                                   aaGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCg                              <  1:3484448/150‑1 (MQ=255)
                                                                                                                     ggCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgac                            <  2:580987/150‑1 (MQ=255)
                                                                                                                      gCATATAGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgacc                           >  2:4115052/1‑150 (MQ=255)
                                                                                                                          ataGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccagac                       <  1:2308764/150‑1 (MQ=255)
                                                                                                                          ataGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccagac                       <  1:2308835/150‑1 (MQ=255)
                                                                                                                           taGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccagacc                      >  1:1098652/1‑150 (MQ=255)
                                                                                                                           taGGATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCCAGAAACCCAAACAGGTCAATAGCCACCAGCaccacccaacc                      >  2:248145/1‑145 (MQ=255)
                                                                                                                             ggATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGaccagaccag                    <  1:3800349/150‑1 (MQ=255)
                                                                                                                              gATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAgg                   <  1:13561/150‑1 (MQ=255)
                                                                                                                              gATGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAgg                   >  2:1406148/1‑150 (MQ=255)
                                                                                                                               aTGTTGTAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGc                  >  2:4081103/1‑150 (MQ=255)
                                                                                                                                tgttgtAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGGTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCa                 >  2:1062413/1‑150 (MQ=255)
                                                                                                                                  ttgtAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACt               <  1:1646796/150‑1 (MQ=255)
                                                                                                                                  ttgtAACTAATATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACt               <  1:2995148/150‑1 (MQ=255)
                                                                                                                                      aaCTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCg           >  1:2911500/1‑150 (MQ=255)
                                                                                                                                          aaTATGGTAAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAg       <  1:70581/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCT  >  CP009273/1603294‑1603586

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: