Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,122,994 C→A 100% R245S (CGC→AGC)  yceM → putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,122,9940CA96.0% 532.6 / 12.3 176R245S (CGC→AGC) yceMputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base A (80/89);  total (82/94)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 181 total aligned reads displayed.

GTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTG  >  CP009273/1122849‑1123137
                                                                                                                                                 |                                                                                                                                               
gTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAgagc                                                                                                                                               <  1:1333532/148‑1 (MQ=255)
gTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAgagc                                                                                                                                               >  2:891063/1‑148 (MQ=255)
  ttGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGACTGGCGTGAGGAGAGCgg                                                                                                                                             >  2:2553495/1‑148 (MQ=255)
   tGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCggg                                                                                                                                            <  1:1912207/148‑1 (MQ=255)
   tGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCggg                                                                                                                                            <  1:2737055/148‑1 (MQ=255)
     ccGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCa                                                                                                                                          >  2:3174983/1‑148 (MQ=255)
        aGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAggg                                                                                                                                       <  1:2501024/148‑1 (MQ=255)
         gCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGc                                                                                                                                      >  2:2068863/1‑148 (MQ=255)
            ccATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTa                                                                                                                                   <  2:3427205/148‑1 (MQ=255)
             cATTTTTCGGCTGGTTCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAg                                                                                                                                  <  1:17055/148‑1 (MQ=255)
              aTTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATCTGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGt                                                                                                                                 >  2:2598561/1‑148 (MQ=255)
                 tttCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCa                                                                                                                              >  2:685415/1‑148 (MQ=255)
                 tttCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCa                                                                                                                              >  2:2968183/1‑148 (MQ=255)
                    cGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATaa                                                                                                                           <  2:3330862/148‑1 (MQ=255)
                        tGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCg                                                                                                                       <  1:104643/148‑1 (MQ=255)
                         ggTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGa                                                                                                                      >  1:273225/1‑148 (MQ=255)
                          gTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGAt                                                                                                                     <  1:2886306/148‑1 (MQ=255)
                          gTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGAt                                                                                                                     <  1:857113/148‑1 (MQ=255)
                           tCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGAtt                                                                                                                    <  1:3117474/148‑1 (MQ=255)
                           tCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGAtt                                                                                                                    <  1:1159857/148‑1 (MQ=255)
                           tCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGAtt                                                                                                                    <  1:1472539/148‑1 (MQ=255)
                            ccTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTc                                                                                                                   <  2:1107641/148‑1 (MQ=255)
                              tttGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCt                                                                                                                 <  2:3256414/148‑1 (MQ=255)
                               ttGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTg                                                                                                                >  1:922719/1‑148 (MQ=255)
                               ttGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTg                                                                                                                >  2:51519/1‑148 (MQ=255)
                                 gCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGt                                                                                                              <  2:191918/148‑1 (MQ=255)
                                          caccTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGt                                                                                                     >  1:9446/1‑148 (MQ=255)
                                          caccTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGt                                                                                                     >  1:2935312/1‑148 (MQ=255)
                                           accTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTa                                                                                                    >  2:3129744/1‑148 (MQ=255)
                                           accTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAAACGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTa                                                                                                    >  2:3129823/1‑148 (MQ=255)
                                            ccTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTAc                                                                                                   >  2:199087/1‑148 (MQ=255)
                                            ccTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGGAGTGCATAAACCGATTCCTGGTTGGCAGAGTAc                                                                                                   >  1:3326432/1‑148 (MQ=255)
                                             cTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGATCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACg                                                                                                  <  1:2582839/148‑1 (MQ=255)
                                             cTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGATCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACg                                                                                                  <  1:2583036/148‑1 (MQ=255)
                                              tGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGc                                                                                                 <  1:919348/148‑1 (MQ=255)
                                              tGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGc                                                                                                 <  1:583754/148‑1 (MQ=255)
                                               gTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCt                                                                                                >  2:1501384/1‑148 (MQ=255)
                                               gTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCt                                                                                                >  2:3361480/1‑148 (MQ=255)
                                                 aTGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTg                                                                                              >  2:3147303/1‑148 (MQ=255)
                                                   gCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAg                                                                                            <  1:2038873/148‑1 (MQ=255)
                                                    cATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGTCATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGc                                                                                           >  1:1858717/1‑148 (MQ=255)
                                                    cATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGc                                                                                           >  1:2745217/1‑148 (MQ=255)
                                                     aTCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGCTATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCa                                                                                          <  2:2272054/148‑1 (MQ=255)
                                                     aTCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCa                                                                                          <  2:1380431/148‑1 (MQ=255)
                                                     aTCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCa                                                                                          <  2:1378725/148‑1 (MQ=255)
                                                      tCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCaa                                                                                         >  2:2956100/1‑148 (MQ=255)
                                                      tCGCCGTGACGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCaa                                                                                         >  2:2359023/1‑148 (MQ=255)
                                                          cGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTg                                                                                     <  1:1058809/148‑1 (MQ=255)
                                                           gTGCCGGATGTCATGGTGATACCGTGCAGGCCGTGACTGACGGTGAGCTCATCGAGATTTCGGTTATGCGCGATTGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGTTTCGTGGTTGGCAGATTACGCTTGTGCAACGTgg                                                                                    <  1:2989822/148‑1 (MQ=255)
                                                              ccGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGtt                                                                                 <  1:1114608/148‑1 (MQ=255)
                                                              ccGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGtt                                                                                 <  1:1110574/148‑1 (MQ=255)
                                                                 gAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGt                                                                              >  1:2875139/1‑148 (MQ=255)
                                                                   aGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCg                                                                            <  2:1694300/148‑1 (MQ=255)
                                                                    gTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCgg                                                                           <  1:3363154/148‑1 (MQ=255)
                                                                    gTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCgg                                                                           <  2:579893/148‑1 (MQ=255)
                                                                     tCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGc                                                                          >  2:2756237/1‑148 (MQ=255)
                                                                         cGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCtgtg                                                                      <  1:3012786/148‑1 (MQ=255)
                                                                             aaaCCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGcgcg                                                                  >  1:2455799/1‑148 (MQ=255)
                                                                              aaCCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCgg                                                                 >  1:2330797/1‑148 (MQ=255)
                                                                                ccGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCa                                                               <  1:1334785/148‑1 (MQ=255)
                                                                                 cGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCAc                                                              <  1:3282284/148‑1 (MQ=255)
                                                                                 cGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCAc                                                              <  1:1519233/148‑1 (MQ=255)
                                                                                 cGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCAc                                                              <  1:3214777/148‑1 (MQ=255)
                                                                                   tGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACtt                                                            >  1:2431989/1‑148 (MQ=255)
                                                                                     cAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCa                                                          <  2:2347460/148‑1 (MQ=255)
                                                                                     cAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCa                                                          <  2:465998/148‑1 (MQ=255)
                                                                                     cAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCa                                                          <  2:465560/148‑1 (MQ=255)
                                                                                      aGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCAt                                                         <  1:754523/148‑1 (MQ=255)
                                                                                      aGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGAGTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCAt                                                         <  1:276937/148‑1 (MQ=255)
                                                                                        gCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTg                                                       >  1:1789685/1‑148 (MQ=255)
                                                                                           gTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAAt                                                    >  1:647110/1‑148 (MQ=255)
                                                                                             gactgacGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAAtgt                                                  >  2:817669/1‑148 (MQ=255)
                                                                                                  acGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCa                                             >  2:3267656/1‑148 (MQ=255)
                                                                                                     gTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCaaaa                                          <  2:585827/148‑1 (MQ=255)
                                                                                                          cTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGa                                     >  1:3509435/1‑148 (MQ=255)
                                                                                                            cATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACa                                   <  2:1181749/148‑1 (MQ=255)
                                                                                                            cATCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACa                                   >  2:3557253/1‑148 (MQ=255)
                                                                                                              tCGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGt                                 <  1:2954413/148‑1 (MQ=255)
                                                                                                               cGACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGtt                                <  1:2196665/148‑1 (MQ=255)
                                                                                                                gACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTc                               <  1:1449288/148‑1 (MQ=255)
                                                                                                                gACATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTc                               <  2:2942271/148‑1 (MQ=255)
                                                                                                                 aCATTACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTcc                              >  1:1674937/1‑148 (MQ=255)
                                                                                                                    ttACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCa                           >  1:2061592/1‑148 (MQ=255)
                                                                                                                     tACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCaa                          >  2:2113497/1‑148 (MQ=255)
                                                                                                                     tACGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCaa                          >  1:3434675/1‑148 (MQ=255)
                                                                                                                      aCGGATATGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCaaa                         <  1:2840475/148‑1 (MQ=255)
                                                                                                                          atatGCGCGAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAAccgc                     >  1:2979926/1‑148 (MQ=255)
                                                                                                                                 cgAATGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGaa              <  1:2708475/148‑1 (MQ=255)
                                                                                                                                   aaTGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACa            >  2:658597/1‑148 (MQ=255)
                                                                                                                                    aTGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAg           >  2:646099/1‑148 (MQ=255)
                                                                                                                                     tGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAgg          >  2:943511/1‑148 (MQ=255)
                                                                                                                                     tGGCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAgg          <  2:1465349/148‑1 (MQ=255)
                                                                                                                                      ggCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGc         >  1:729047/1‑148 (MQ=255)
                                                                                                                                       gCGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGcc        >  1:27906/1‑148 (MQ=255)
                                                                                                                                        cGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCg       <  1:2128029/148‑1 (MQ=255)
                                                                                                                                        cGTGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCg       >  2:3361553/1‑148 (MQ=255)
                                                                                                                                          tGAGGAGAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTg     >  2:2871399/1‑148 (MQ=255)
                                                                                                                                           gaggagAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGGCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGc    <  1:51519/148‑1 (MQ=255)
                                                                                                                                            aggagAGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCt   >  2:90976/1‑148 (MQ=255)
                                                                                                                                             ggagAGCAGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTg  >  2:1570909/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTG  >  CP009273/1122849‑1123137

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: