Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,574,523 T→A 100% L175F (TTA→TTT)  ydeM ← putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,574,5230TA96.1% 906.9 / 19.3 282L175F (TTA→TTT) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base T (9/2);  new base A (144/127);  total (153/129)
Fisher's exact test for biased strand distribution p-value = 7.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 285 total aligned reads displayed.

TGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTG  >  CP009273/1574379‑1574667
                                                                                                                                                |                                                                                                                                                
     tAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATg                                                                                                                                          <  1:1926236/148‑1 (MQ=255)
      aaaCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGa                                                                                                                                         <  1:464497/148‑1 (MQ=255)
      aaaCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGa                                                                                                                                         >  1:1706774/1‑148 (MQ=255)
       aaCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGAt                                                                                                                                        >  1:127627/1‑148 (MQ=255)
       aaCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGAt                                                                                                                                        >  1:2533346/1‑148 (MQ=255)
        aCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATa                                                                                                                                       >  1:198689/1‑148 (MQ=255)
         cTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATaa                                                                                                                                      <  1:1943899/148‑1 (MQ=255)
          ttGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATaaa                                                                                                                                     >  1:3331917/1‑148 (MQ=255)
           tGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAAc                                                                                                                                    <  1:1147860/148‑1 (MQ=255)
           tGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAAc                                                                                                                                    <  1:3592167/148‑1 (MQ=255)
               aTAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCa                                                                                                                                <  2:3759934/148‑1 (MQ=255)
               aTAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCa                                                                                                                                <  1:499073/148‑1 (MQ=255)
                 aaGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAag                                                                                                                              >  2:2657416/1‑148 (MQ=255)
                  aGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAaga                                                                                                                             <  1:1841526/148‑1 (MQ=255)
                  aGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAaga                                                                                                                             >  2:3595515/1‑148 (MQ=255)
                       tGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTa                                                                                                                        <  1:3606008/148‑1 (MQ=255)
                       tGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTa                                                                                                                        <  1:530989/148‑1 (MQ=255)
                       tGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTa                                                                                                                        <  1:1517174/148‑1 (MQ=255)
                       tGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTa                                                                                                                        <  2:376850/148‑1 (MQ=255)
                         ggAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAAt                                                                                                                      >  1:2078813/1‑148 (MQ=255)
                            ggCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGa                                                                                                                   <  1:680363/148‑1 (MQ=255)
                             gCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGAc                                                                                                                  <  2:1396801/148‑1 (MQ=255)
                              cacaGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACa                                                                                                                 >  2:794569/1‑148 (MQ=255)
                                caGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACAtt                                                                                                               >  2:2972861/1‑148 (MQ=255)
                                 aGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTa                                                                                                              >  1:2971458/1‑148 (MQ=255)
                                 aGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTa                                                                                                              <  1:1793483/148‑1 (MQ=255)
                                     aaaTCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGt                                                                                                          <  2:41316/148‑1 (MQ=255)
                                     aaaTCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGt                                                                                                          >  1:3539662/1‑148 (MQ=255)
                                     aaaTCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGt                                                                                                          >  2:2481421/1‑148 (MQ=255)
                                      aaTCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGtt                                                                                                         <  2:2649561/148‑1 (MQ=255)
                                       aTCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGTATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGtta                                                                                                        >  1:159271/1‑148 (MQ=255)
                                        tCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGttat                                                                                                       >  1:754249/1‑148 (MQ=255)
                                          aaTGATACGGAATTTGTTCTAACTATTACCACTGAAATCAATATTCGGCGTCCAGGTTTCTAGCAATTCGATAAATTGCTTATGTTTACTGCCTATAGATTTAATAAAATCATAAACCTCAAGAGGGTAATTGACATTAACGttatta                                                                                                     <  2:3126843/148‑1 (MQ=255)
                                          aaTGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGttatta                                                                                                     <  2:877928/148‑1 (MQ=255)
                                          aaTGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGttatta                                                                                                     <  2:102825/148‑1 (MQ=255)
                                            tGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAAt                                                                                                   <  2:225590/148‑1 (MQ=255)
                                                aCGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACg                                                                                               >  1:1578271/1‑148 (MQ=255)
                                                 cGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACgg                                                                                              <  2:488362/148‑1 (MQ=255)
                                                    aaTGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTCATGGACATTAACGTGATTATTGACGGtga                                                                                           >  1:2614230/1‑146 (MQ=255)
                                                    aaTGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTa                                                                                           <  2:323647/148‑1 (MQ=255)
                                                       gtgtTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACg                                                                                        <  2:2358580/148‑1 (MQ=255)
                                                        tgtTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGt                                                                                       <  1:1462251/148‑1 (MQ=255)
                                                         gtTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTa                                                                                      <  1:362547/148‑1 (MQ=255)
                                                              cACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAACCCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTaaa                                                                                 <  2:739773/148‑1 (MQ=255)
                                                              cACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTaaa                                                                                 >  2:1519527/1‑148 (MQ=255)
                                                              cACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTaaa                                                                                 >  1:3271182/1‑148 (MQ=255)
                                                                cTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAAct                                                                               >  2:1459578/1‑148 (MQ=255)
                                                                 tATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAActc                                                                              <  2:2310259/148‑1 (MQ=255)
                                                                  aTGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAActct                                                                             <  1:1001712/148‑1 (MQ=255)
                                                                  aTGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAActct                                                                             >  2:1919562/1‑148 (MQ=255)
                                                                   tGACCACTGAACTCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTa                                                                            <  2:2498651/148‑1 (MQ=255)
                                                                    gACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTAc                                                                           >  2:1816773/1‑148 (MQ=255)
                                                                    gACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTAc                                                                           >  2:3543565/1‑148 (MQ=255)
                                                                       cACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTg                                                                        >  2:3404024/1‑148 (MQ=255)
                                                                        aCTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGa                                                                       >  2:3093881/1‑148 (MQ=255)
                                                                          tGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGata                                                                     >  2:2585293/1‑148 (MQ=255)
                                                                          tGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGata                                                                     >  2:2419499/1‑148 (MQ=255)
                                                                             aaTCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGa                                                                  >  2:1119386/1‑148 (MQ=255)
                                                                              aTCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGAt                                                                 <  2:2134953/148‑1 (MQ=255)
                                                                               tCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATTACGGTTAACGTATTAAACTCTACTTGATATGAtt                                                                >  1:3492692/1‑148 (MQ=255)
                                                                               tCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGAtt                                                                >  2:502279/1‑148 (MQ=255)
                                                                                cAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGAttt                                                               <  1:348617/148‑1 (MQ=255)
                                                                                cAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGAttt                                                               >  2:3500004/1‑148 (MQ=255)
                                                                                cAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGAttt                                                               >  2:851666/1‑148 (MQ=255)
                                                                                cAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGAAGGTTAACGTATTAAACTCTACTTGATATGAttt                                                               <  1:1063345/148‑1 (MQ=255)
                                                                                cAATATTCGGCGTCCCGGTTTCTAGCAATACGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGATATTAATGACGGTTAACGTATTAAACTCTACTTGATATGAttt                                                               >  2:1965219/1‑148 (MQ=255)
                                                                                 aaTATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTc                                                              <  1:1545924/148‑1 (MQ=255)
                                                                                  atatTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCa                                                             >  2:1709159/1‑148 (MQ=255)
                                                                                   tatTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAg                                                            <  1:1909868/148‑1 (MQ=255)
                                                                                   tatTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAg                                                            <  1:3789877/148‑1 (MQ=255)
                                                                                    atTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGa                                                           >  2:1587639/1‑148 (MQ=255)
                                                                                    atTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGa                                                           >  1:3408887/1‑148 (MQ=255)
                                                                                     ttCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGAc                                                          >  1:800682/1‑148 (MQ=255)
                                                                                      tCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACg                                                         >  2:3664523/1‑148 (MQ=255)
                                                                                          cGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATATAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTATACTTGATATCATTTCAGACGCTCg                                                     >  2:2088581/1‑148 (MQ=255)
                                                                                           gTCCCGGTTTCTAGCAATTCGGTAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGa                                                    <  2:271061/148‑1 (MQ=255)
                                                                                            tCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGAt                                                   <  1:2892982/148‑1 (MQ=255)
                                                                                            tCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGAt                                                   >  2:2967888/1‑148 (MQ=255)
                                                                                             cccGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATg                                                  >  2:2069400/1‑148 (MQ=255)
                                                                                                 gTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGgctg                                              >  1:2800519/1‑148 (MQ=255)
                                                                                                  tttCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGgctgc                                             >  2:1227440/1‑148 (MQ=255)
                                                                                                   ttCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGgctgct                                            <  1:1643949/148‑1 (MQ=255)
                                                                                                    tCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTa                                           <  1:713836/148‑1 (MQ=255)
                                                                                                    tCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGCACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTa                                           <  1:2248896/148‑1 (MQ=255)
                                                                                                         cAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACt                                      >  1:3001236/1‑148 (MQ=255)
                                                                                                           aTTCGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACttt                                    <  1:3141866/148‑1 (MQ=255)
                                                                                                              cGATAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGc                                 >  1:2721782/1‑148 (MQ=255)
                                                                                                                 tAAATTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCaaa                              >  1:3229154/1‑148 (MQ=255)
                                                                                                                    aTTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGt                           >  1:876020/1‑148 (MQ=255)
                                                                                                                    aTTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGt                           >  2:3175362/1‑148 (MQ=255)
                                                                                                                    aTTGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGt                           >  1:2321195/1‑148 (MQ=255)
                                                                                                                     ttGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTa                          >  2:136755/1‑148 (MQ=255)
                                                                                                                     ttGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTa                          >  2:3483084/1‑148 (MQ=255)
                                                                                                                      tGCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTAc                         <  2:3703517/148‑1 (MQ=255)
                                                                                                                      ggCATATGTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTAc                         <  2:666717/147‑1 (MQ=255)
                                                                                                                             gTTTACTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTAcc                  <  1:2992650/148‑1 (MQ=255)
                                                                                                                                 aCTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGaa              <  2:3043021/148‑1 (MQ=255)
                                                                                                                                 aCTGCCGATAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGaa              <  2:3343126/148‑1 (MQ=255)
                                                                                                                                         tAGATTTAAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCg      >  2:3189731/1‑148 (MQ=255)
                                                                                                                                             tttaAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTg  >  2:2905422/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTG  >  CP009273/1574379‑1574667

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: ‑