Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 828,304 A→T 6.2% W35R (TGG→AGG)  ybiA ← DUF1768 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273828,3040AT 6.2% 499.8 / 19.1 160W35R (TGG→AGG) ybiADUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base A (72/78);  new base T (4/6);  total (76/84)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 161 total aligned reads displayed.

TTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAA  >  CP009273/828159‑828451
                                                                                                                                                 |                                                                                                                                                   
 tAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAggc                                                                                                                                                  >  2:1993886/1‑148 (MQ=255)
 tAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAggc                                                                                                                                                  >  1:2079268/1‑148 (MQ=255)
 tAGAACGATCGCGCCCCATGCGTGCGGCGACCAGAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAggc                                                                                                                                                  >  1:144680/1‑148 (MQ=255)
    aaCGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAggcggc                                                                                                                                               >  1:761540/1‑148 (MQ=255)
    aaCGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAggcggc                                                                                                                                               >  1:2874642/1‑148 (MQ=255)
          cgcgCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGtt                                                                                                                                         <  2:1987720/148‑1 (MQ=255)
                   tgcgtgcgGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATc                                                                                                                                >  2:2968652/1‑148 (MQ=255)
                   tgcgtgcgGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATc                                                                                                                                >  2:316776/1‑148 (MQ=255)
                    gcgtgcgGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCa                                                                                                                               <  2:3093557/148‑1 (MQ=255)
                    gcgtgcgGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCa                                                                                                                               <  1:792150/148‑1 (MQ=255)
                    gcgtgcgGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCa                                                                                                                               >  2:1088068/1‑148 (MQ=255)
                    gcgtgcgGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCa                                                                                                                               <  2:1678429/148‑1 (MQ=255)
                      gtgcgGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCAcc                                                                                                                             <  2:706768/148‑1 (MQ=255)
                          ggcgACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAg                                                                                                                         >  2:2156918/1‑148 (MQ=255)
                           gcgACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGt                                                                                                                        >  2:255186/1‑148 (MQ=255)
                           gcgACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGt                                                                                                                        <  2:2453713/148‑1 (MQ=255)
                           gcgACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGt                                                                                                                        <  2:432779/148‑1 (MQ=255)
                            cgACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCTGGCGGCAAAGTTAGAAAAATCACCGTAGTc                                                                                                                       <  1:1818042/148‑1 (MQ=255)
                            cgACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCACAGGTAGAAAAATCACCGTAGTc                                                                                                                       >  2:1946864/1‑148 (MQ=255)
                              aCCATAGGGGAAGAAACCCGACTGATCTCTTCGCGGTATTTTTCATAAAGGATTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCAt                                                                                                                     <  1:1785343/148‑1 (MQ=255)
                                cATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCg                                                                                                                   <  2:2125501/148‑1 (MQ=255)
                                    ggggAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTgg                                                                                                               >  1:1858285/1‑148 (MQ=255)
                                      ggAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTg                                                                                                             <  1:2968652/148‑1 (MQ=255)
                                      ggAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTg                                                                                                             <  1:316776/148‑1 (MQ=255)
                                       gaagaaACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGc                                                                                                            >  2:2331898/1‑148 (MQ=255)
                                           aaaCCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTa                                                                                                        >  1:2521481/1‑148 (MQ=255)
                                                cGACAGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAAt                                                                                                   >  2:1971896/1‑148 (MQ=255)
                                                  aCGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAAttt                                                                                                 <  2:2598462/148‑1 (MQ=255)
                                                  aCGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAAttt                                                                                                 <  2:2598819/148‑1 (MQ=255)
                                                   cGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTa                                                                                                <  1:2078989/148‑1 (MQ=255)
                                                     gATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTatg                                                                                              >  2:338488/1‑148 (MQ=255)
                                                     gATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTatg                                                                                              <  2:1103616/148‑1 (MQ=255)
                                                     gATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTatg                                                                                              >  1:1541716/1‑148 (MQ=255)
                                                      aTCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTatga                                                                                             >  2:1360317/1‑148 (MQ=255)
                                                      aTCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTatga                                                                                             >  2:1360170/1‑148 (MQ=255)
                                                       tctcTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTatgat                                                                                            >  1:1049399/1‑148 (MQ=255)
                                                         tctTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTGGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCTTAGTCATCGCTGGTGCTGTAGAAATTTATGATgg                                                                                          <  2:672652/148‑1 (MQ=255)
                                                         tctTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTGGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCTTAGTCATCGCTGGTGCTGTAGAAATTTATGATgg                                                                                          <  2:2332877/148‑1 (MQ=255)
                                                          ctTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGgt                                                                                         >  2:1937404/1‑148 (MQ=255)
                                                            tCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGgtgt                                                                                       >  2:2082255/1‑148 (MQ=255)
                                                             cgcgGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTc                                                                                      <  1:673876/148‑1 (MQ=255)
                                                               cgGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCtt                                                                                    >  1:1784049/1‑148 (MQ=255)
                                                                  tATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCTGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCa                                                                                 <  2:1073157/148‑1 (MQ=255)
                                                                  tATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCa                                                                                 >  1:1461880/1‑148 (MQ=255)
                                                                  tATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCa                                                                                 <  1:984010/148‑1 (MQ=255)
                                                                   aTTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCAt                                                                                >  1:3311870/1‑148 (MQ=255)
                                                                   aTTTTTCATCAAGGAATTTTTGTGCCCGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCAt                                                                                >  2:2939703/1‑148 (MQ=255)
                                                                   aTTTTACATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCAt                                                                                <  1:3163582/148‑1 (MQ=255)
                                                                    tttttCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCAtt                                                                               >  1:1041062/1‑147 (MQ=255)
                                                                     ttttCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCTGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCAtta                                                                              >  1:2694466/1‑146 (MQ=255)
                                                                       ttCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCTGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACg                                                                            <  1:1156642/148‑1 (MQ=255)
                                                                        tcatcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACgt                                                                           >  1:2246800/1‑148 (MQ=255)
                                                                         catcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACgtg                                                                          >  2:2325324/1‑148 (MQ=255)
                                                                         catcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACgtg                                                                          <  2:3244352/148‑1 (MQ=255)
                                                                         catcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACgtg                                                                          <  2:3245796/148‑1 (MQ=255)
                                                                          atcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCGTTACgtgt                                                                         <  2:1252585/148‑1 (MQ=255)
                                                                             aaGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTgg                                                                      <  1:970293/148‑1 (MQ=255)
                                                                             aaGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTgg                                                                      <  1:2943585/148‑1 (MQ=255)
                                                                              aGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGa                                                                     >  2:1503500/1‑148 (MQ=255)
                                                                                 aaTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATcc                                                                  <  1:3224193/148‑1 (MQ=255)
                                                                                 aaTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATcc                                                                  <  1:996180/148‑1 (MQ=255)
                                                                                 aaTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATcc                                                                  <  1:3054637/148‑1 (MQ=255)
                                                                                   tttttGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCtt                                                                >  2:2789979/1‑148 (MQ=255)
                                                                                   tttttGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCtt                                                                >  2:2789739/1‑148 (MQ=255)
                                                                                   tttttGGGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCtt                                                                <  2:655472/148‑1 (MQ=255)
                                                                                    ttttGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCttt                                                               <  2:544714/148‑1 (MQ=255)
                                                                                        gtgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGc                                                           <  2:3184217/148‑1 (MQ=255)
                                                                                         tgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCt                                                          <  2:1924457/148‑1 (MQ=255)
                                                                                            cTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGa                                                       <  2:1041062/148‑1 (MQ=255)
                                                                                             tGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGaa                                                      <  1:547211/148‑1 (MQ=255)
                                                                                                aaaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCTGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACgg                                                   >  2:27166/1‑148 (MQ=255)
                                                                                                   tAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCa                                                <  1:2262293/148‑1 (MQ=255)
                                                                                                   tAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCa                                                >  1:2307276/1‑148 (MQ=255)
                                                                                                    aGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCAt                                               <  2:660573/148‑1 (MQ=255)
                                                                                                     gtgtTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCAtt                                              >  1:1190752/1‑148 (MQ=255)
                                                                                                        ttCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTaaa                                           <  2:542732/148‑1 (MQ=255)
                                                                                                            gAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCTGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGc                                       <  1:1934161/148‑1 (MQ=255)
                                                                                                             aGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCTGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCa                                      <  1:2441550/148‑1 (MQ=255)
                                                                                                             aGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCa                                      <  2:2340167/148‑1 (MQ=255)
                                                                                                             aGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCa                                      <  2:2521481/148‑1 (MQ=255)
                                                                                                                 aGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTAt                                  <  2:1086905/148‑1 (MQ=255)
                                                                                                                  ggCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATg                                 >  1:1982532/1‑148 (MQ=255)
                                                                                                                  ggCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATg                                 <  2:1310795/148‑1 (MQ=255)
                                                                                                                      aGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtgg                             >  2:1736843/1‑148 (MQ=255)
                                                                                                                       ggTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggt                            >  2:299124/1‑148 (MQ=255)
                                                                                                                         ttttCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTAATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtgg                          >  2:713569/1‑148 (MQ=255)
                                                                                                                          tttCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTggtggt                         <  1:2224620/148‑1 (MQ=255)
                                                                                                                            tCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTa                       <  1:2497236/148‑1 (MQ=255)
                                                                                                                             cccGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAt                      <  1:3293321/148‑1 (MQ=255)
                                                                                                                             cccGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAt                      <  1:3293260/148‑1 (MQ=255)
                                                                                                                               cGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTg                    <  2:382558/148‑1 (MQ=255)
                                                                                                                                   aaCCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAg                <  1:1782711/148‑1 (MQ=255)
                                                                                                                                      cTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCg             <  1:2592229/148‑1 (MQ=255)
                                                                                                                                      cTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCg             <  2:1541716/148‑1 (MQ=255)
                                                                                                                                      cTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCg             <  2:1541695/148‑1 (MQ=255)
                                                                                                                                      cTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCg             <  1:280253/148‑1 (MQ=255)
                                                                                                                                       ttGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGc            <  1:195549/148‑1 (MQ=255)
                                                                                                                                        tGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCt           >  1:1398174/1‑148 (MQ=255)
                                                                                                                                          aTGGGCCTGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCttt         >  1:1597661/1‑148 (MQ=255)
                                                                                                                                           tGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTc        >  2:770108/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAA  >  CP009273/828159‑828451

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: