Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 524,496 T→G 100% pseudogene (401/491 nt) ybbD → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273524,4960TG95.7% 1073.2 / 14.2 279pseudogene (401/491 nt)ybbDpseudogene
Reads supporting (aligned to +/- strand):  ref base T (9/3);  new base G (133/134);  total (142/137)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 281 total aligned reads displayed.

CCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGA  >  CP009273/524351‑524640
                                                                                                                                                 |                                                                                                                                                
   tAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTaaa                                                                                                                                             >  1:2188893/1‑148 (MQ=255)
       aaaTGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtt                                                                                                                                         <  2:2194432/148‑1 (MQ=255)
          tgatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCttttt                                                                                                                                      >  1:556445/1‑148 (MQ=255)
            atgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCttttttt                                                                                                                                    <  2:1965472/148‑1 (MQ=255)
             tgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtttttttt                                                                                                                                   >  2:1583615/1‑148 (MQ=255)
              gatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttg                                                                                                                                  >  2:3315401/1‑148 (MQ=255)
               atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttgt                                                                                                                                 <  2:3070505/148‑1 (MQ=255)
               atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttgt                                                                                                                                 >  2:149202/1‑148 (MQ=255)
               atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttgt                                                                                                                                 <  2:3399930/148‑1 (MQ=255)
               atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttgt                                                                                                                                 <  2:1295271/148‑1 (MQ=255)
                   gggTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGt                                                                                                                             <  2:1447895/148‑1 (MQ=255)
                       aGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGaaa                                                                                                                         <  2:689466/148‑1 (MQ=255)
                        ggTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAAt                                                                                                                        >  1:1467954/1‑148 (MQ=255)
                         gTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTTGTGAAAtt                                                                                                                       <  2:2813535/148‑1 (MQ=255)
                          tttAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTa                                                                                                                      >  1:213484/1‑148 (MQ=255)
                           taaataaaaCATTGACCTTGAGGGCCAGGTAACATGGGGCGAAAACCGTTCAAAAGAATCATCCAGATACAAAAAAAGTTTGTTTTAAATATGGTAAGATACCAGGGCAAGGGTTTGGGTTAAACCTTTTTTTTGTTGGTGAAATTaa                                                                                                                     <  2:744358/143‑1 (MQ=255)
                            tAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAAt                                                                                                                    <  1:2129642/148‑1 (MQ=255)
                             aaGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATg                                                                                                                   <  1:2934922/148‑1 (MQ=255)
                                aaaaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTa                                                                                                                <  2:1894398/148‑1 (MQ=255)
                                    cATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCa                                                                                                            <  2:3468144/148‑1 (MQ=255)
                                     aTTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCat                                                                                                           <  1:780721/148‑1 (MQ=255)
                                      ttGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCatt                                                                                                          >  1:89206/1‑148 (MQ=255)
                                      ttGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCatt                                                                                                          >  1:1891685/1‑148 (MQ=255)
                                        gACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCattat                                                                                                        <  1:2000181/148‑1 (MQ=255)
                                        gACAGTGAGGGCCAGGGAACAAGTGGAGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCattat                                                                                                        <  2:2293491/148‑1 (MQ=255)
                                          cAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTAttt                                                                                                      <  2:2902580/148‑1 (MQ=255)
                                             tGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTtat                                                                                                   >  2:921625/1‑148 (MQ=255)
                                              gAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTtata                                                                                                  <  2:1907153/148‑1 (MQ=255)
                                                 ggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACa                                                                                               <  1:1735257/148‑1 (MQ=255)
                                                  gCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACaa                                                                                              >  1:620795/1‑148 (MQ=255)
                                                   ccAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACaaa                                                                                             <  1:309269/148‑1 (MQ=255)
                                                    cAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAat                                                                                            <  2:867599/148‑1 (MQ=255)
                                                    cAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAat                                                                                            <  2:1462701/148‑1 (MQ=255)
                                                     aGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAata                                                                                           <  1:3369620/148‑1 (MQ=255)
                                                       ggAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATAtt                                                                                         <  1:634042/148‑1 (MQ=255)
                                                        gAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTg                                                                                        >  2:1687102/1‑148 (MQ=255)
                                                        gAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTg                                                                                        >  2:1066298/1‑148 (MQ=255)
                                                         aaCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTgg                                                                                       <  1:1331631/148‑1 (MQ=255)
                                                         aaCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTgg                                                                                       <  2:1741589/148‑1 (MQ=255)
                                                          aCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGa                                                                                      >  1:2340266/1‑148 (MQ=255)
                                                              gTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTg                                                                                  >  2:2984605/1‑148 (MQ=255)
                                                                ggCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGa                                                                                >  2:2324143/1‑148 (MQ=255)
                                                                ggCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGa                                                                                >  2:1207490/1‑148 (MQ=255)
                                                                 gCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAt                                                                               >  2:1846743/1‑148 (MQ=255)
                                                                  cGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAtt                                                                              <  1:1370118/148‑1 (MQ=255)
                                                                   gAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTg                                                                             >  2:2548330/1‑148 (MQ=255)
                                                                     aaaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGcc                                                                           >  1:3026860/1‑148 (MQ=255)
                                                                     aaaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGcc                                                                           <  2:283867/148‑1 (MQ=255)
                                                                     aaaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGcc                                                                           >  1:2478211/1‑148 (MQ=255)
                                                                      aaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                                                                          >  1:256264/1‑148 (MQ=255)
                                                                      aaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                                                                          >  1:1877427/1‑148 (MQ=255)
                                                                       aTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTg                                                                         <  1:940207/148‑1 (MQ=255)
                                                                          gTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATg                                                                      >  1:217407/1‑148 (MQ=255)
                                                                           tATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGc                                                                     >  1:3126887/1‑148 (MQ=255)
                                                                           tATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGc                                                                     >  1:168929/1‑148 (MQ=255)
                                                                           tATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGc                                                                     <  1:1021133/148‑1 (MQ=255)
                                                                              cAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtg                                                                  >  2:960952/1‑148 (MQ=255)
                                                                                  gAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCa                                                              >  2:3300953/1‑148 (MQ=255)
                                                                                   aaTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                                             >  1:1199007/1‑148 (MQ=255)
                                                                                   aaTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                                             >  2:3206251/1‑148 (MQ=255)
                                                                                    aTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCaga                                                            <  1:2985076/148‑1 (MQ=255)
                                                                                    aTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCaga                                                            <  1:321197/148‑1 (MQ=255)
                                                                                      gATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGt                                                          <  2:1855156/148‑1 (MQ=255)
                                                                                        tCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTaa                                                        <  2:680105/148‑1 (MQ=255)
                                                                                         ccAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAt                                                       <  2:1658000/148‑1 (MQ=255)
                                                                                         ccAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAt                                                       <  1:2615295/148‑1 (MQ=255)
                                                                                           aGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTg                                                     <  1:1116433/148‑1 (MQ=255)
                                                                                               aCAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTg                                                 >  2:2401430/1‑148 (MQ=255)
                                                                                                cAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGa                                                <  2:386146/148‑1 (MQ=255)
                                                                                                 aaaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGaa                                               <  2:3104364/148‑1 (MQ=255)
                                                                                                  aaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAAt                                              >  1:2885361/1‑148 (MQ=255)
                                                                                                   aaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAAtt                                             <  2:1000660/148‑1 (MQ=255)
                                                                                                   aaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAAtt                                             >  1:400373/1‑148 (MQ=255)
                                                                                                     aaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTaa                                           <  2:3497269/148‑1 (MQ=255)
                                                                                                     aaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTaa                                           >  2:2126044/1‑148 (MQ=255)
                                                                                                     aaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTaa                                           <  1:2380035/148‑1 (MQ=255)
                                                                                                         tgtgGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGt                                       <  2:544355/148‑1 (MQ=255)
                                                                                                         tgtgGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGt                                       <  2:489049/148‑1 (MQ=255)
                                                                                                           tggtgTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCt                                     <  2:615732/148‑1 (MQ=255)
                                                                                                            ggtgTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACACATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTa                                    >  2:1975863/1‑148 (MQ=255)
                                                                                                             gtgtAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTAt                                   >  2:2684266/1‑148 (MQ=255)
                                                                                                                tAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATccc                                <  1:2720226/148‑1 (MQ=255)
                                                                                                                  aaTATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCat                              <  1:3559129/148‑1 (MQ=255)
                                                                                                                   atatGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCata                             >  1:728157/1‑148 (MQ=255)
                                                                                                                     atGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATc                           >  2:2756100/1‑148 (MQ=255)
                                                                                                                     atGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATc                           >  2:2669177/1‑148 (MQ=255)
                                                                                                                       ggTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGa                         >  2:381171/1‑148 (MQ=255)
                                                                                                                        gTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGaa                        >  2:3521241/1‑148 (MQ=255)
                                                                                                                         tAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAg                       <  2:2600574/148‑1 (MQ=255)
                                                                                                                         tAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAg                       <  1:3315401/148‑1 (MQ=255)
                                                                                                                            gATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAgtcg                    >  1:2833088/1‑148 (MQ=255)
                                                                                                                             aTACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAgtcgt                   >  1:1437459/1‑148 (MQ=255)
                                                                                                                               aCCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCa                 >  1:3489107/1‑148 (MQ=255)
                                                                                                                                ccAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCaa                <  1:438065/148‑1 (MQ=255)
                                                                                                                                ccAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCaa                <  1:1583615/148‑1 (MQ=255)
                                                                                                                                ccAGGGCAAGGGGATGGGGTAAACCATTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTACCCCATATCGAAGTCGTCaa                <  2:1774639/148‑1 (MQ=255)
                                                                                                                                   gggCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACtt             <  2:1955764/148‑1 (MQ=255)
                                                                                                                                      cAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGt          <  2:1165977/148‑1 (MQ=255)
                                                                                                                                       aaGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTa         <  1:2299768/148‑1 (MQ=255)
                                                                                                                                         ggggATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGt       <  1:3310260/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGA  >  CP009273/524351‑524640

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: