Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,684,105 C→A 100% intergenic (+97/‑4) manA → / → ydgA mannose‑6‑phosphate isomerase/DUF945 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,684,1050CA99.1% 744.5 / ‑2.0 229intergenic (+97/‑4)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (124/103);  total (125/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 236 total aligned reads displayed.

GACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGAC  >  CP009273/1683960‑1684248
                                                                                                                                                 |                                                                                                                                               
gACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATaataat                                                                                                                                               <  2:1616938/148‑1 (MQ=255)
  cTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTa                                                                                                                                             <  2:2905868/148‑1 (MQ=255)
     tCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGa                                                                                                                                          <  2:3497992/148‑1 (MQ=255)
     tCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGa                                                                                                                                          <  2:1550723/148‑1 (MQ=255)
     tCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGa                                                                                                                                          <  2:3612948/148‑1 (MQ=255)
       aaaGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAAt                                                                                                                                        <  2:3405080/148‑1 (MQ=255)
         aGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATaa                                                                                                                                      <  2:526435/148‑1 (MQ=255)
          ggCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATaaa                                                                                                                                     <  1:3606337/148‑1 (MQ=255)
                ggCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTg                                                                                                                               <  2:753434/148‑1 (MQ=255)
                  ccGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGt                                                                                                                             <  2:1727883/148‑1 (MQ=255)
                   cGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTa                                                                                                                            <  2:722175/148‑1 (MQ=255)
                   cGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTa                                                                                                                            <  1:3010195/148‑1 (MQ=255)
                     tttAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGc                                                                                                                          >  2:1242576/1‑148 (MQ=255)
                      ttAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCg                                                                                                                         >  1:2148936/1‑148 (MQ=255)
                       tAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCgg                                                                                                                        <  1:1727285/148‑1 (MQ=255)
                       tAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCgg                                                                                                                        >  2:2204288/1‑148 (MQ=255)
                        aGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGt                                                                                                                       <  2:510306/148‑1 (MQ=255)
                        aGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGt                                                                                                                       <  2:3804218/148‑1 (MQ=255)
                         gcgcGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTa                                                                                                                      <  2:1620787/148‑1 (MQ=255)
                              tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGt                                                                                                                 <  1:610066/148‑1 (MQ=255)
                                 ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                                                                                              <  2:964849/148‑1 (MQ=255)
                                   acaacaAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGTCTTATAACACCTTCAGGCGGCCAGTCCGCATGATTACATTTTATGGATAATAATTATGAATAAATCGCTGGAAGCGGTAGGCGTCAttg                                                                                                            >  1:986787/1‑148 (MQ=255)
                                        aaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcg                                                                                                       <  1:1974995/148‑1 (MQ=255)
                                           cTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATTAATAAATAGCTGGTAGCGGTCGGCGTCATTGTTGCGCTa                                                                                                    >  1:2607114/1‑148 (MQ=255)
                                             gTAAGAGCTTACTGAAAAAATAAACATCTCTTGCTAAGCTTATTAAAGGCTAATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAGTTTCTGGATCATAATTATGAATAAATCGCGGGTACCGGGAGGACGCATTGTGGCGCTAgg                                                                                                  >  2:2309366/1‑148 (MQ=255)
                                              tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGc                                                                                                 >  1:1001024/1‑148 (MQ=255)
                                              tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGc                                                                                                 >  1:2309814/1‑148 (MQ=255)
                                                  agCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                                                                                             >  2:2556165/1‑148 (MQ=255)
                                                    ctaaCTCAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAATTTATGGATAATAAATAAGAATAAATCGCTGGAAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTc                                                                                           >  2:3470758/4‑148 (MQ=255)
                                                     ttACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCt                                                                                          <  1:1502189/148‑1 (MQ=255)
                                                       aCTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg                                                                                        >  1:2800312/1‑148 (MQ=255)
                                                        cTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                                                                                       <  1:607530/148‑1 (MQ=255)
                                                        cTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                                                                                       <  1:944829/148‑1 (MQ=255)
                                                        cTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                                                                                       <  1:3107786/148‑1 (MQ=255)
                                                          gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa                                                                                     <  2:3005593/148‑1 (MQ=255)
                                                           aaaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAg                                                                                    >  2:1228813/1‑148 (MQ=255)
                                                           aaaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAg                                                                                    >  1:229077/1‑148 (MQ=255)
                                                            aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAgg                                                                                   <  1:3487446/148‑1 (MQ=255)
                                                            aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAgg                                                                                   <  1:1452167/148‑1 (MQ=255)
                                                             aaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc                                                                                  <  1:1884525/148‑1 (MQ=255)
                                                             aaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCAGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc                                                                                  >  1:1333248/1‑148 (MQ=255)
                                                               aaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcgg                                                                                >  1:1137615/1‑148 (MQ=255)
                                                               aaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcgg                                                                                >  2:1789028/1‑148 (MQ=255)
                                                                 ttAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg                                                                              >  2:1351653/1‑148 (MQ=255)
                                                                  tAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGgcgc                                                                             >  1:1735329/1‑148 (MQ=255)
                                                                  tAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGgcgc                                                                             >  1:767768/1‑148 (MQ=255)
                                                                  tAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGgcgc                                                                             >  1:2422158/1‑148 (MQ=255)
                                                                  tAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGgcgc                                                                             >  2:2864265/1‑148 (MQ=255)
                                                                  tAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGgcgc                                                                             >  1:3097106/1‑148 (MQ=255)
                                                                    aCATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAt                                                                           <  2:1014089/148‑1 (MQ=255)
                                                                     cATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATg                                                                          <  2:2321144/148‑1 (MQ=255)
                                                                     cATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATg                                                                          <  2:2234297/148‑1 (MQ=255)
                                                                      aTCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATgg                                                                         >  1:1210932/1‑148 (MQ=255)
                                                                      aTCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATgg                                                                         >  2:2996572/1‑148 (MQ=255)
                                                                       tctcTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGt                                                                        <  1:124436/148‑1 (MQ=255)
                                                                       tctcTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGt                                                                        <  1:494043/148‑1 (MQ=255)
                                                                         tctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGtat                                                                      <  2:2521141/148‑1 (MQ=255)
                                                                         tctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGtat                                                                      <  2:2380855/148‑1 (MQ=255)
                                                                         tctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGtat                                                                      <  1:3445512/148‑1 (MQ=255)
                                                                          ctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGtata                                                                     >  2:2624416/1‑148 (MQ=255)
                                                                          ctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCATGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGtata                                                                     >  2:2624393/1‑148 (MQ=255)
                                                                              cTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAgg                                                                 <  1:2937393/148‑1 (MQ=255)
                                                                                aaGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCa                                                               >  1:3178313/1‑148 (MQ=255)
                                                                                aaGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCa                                                               >  1:728630/1‑148 (MQ=255)
                                                                                   cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaaga                                                            <  2:2200445/148‑1 (MQ=255)
                                                                                    ttattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaagaa                                                           <  2:3066692/148‑1 (MQ=255)
                                                                                    ttattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaagaa                                                           >  2:326709/1‑148 (MQ=255)
                                                                                      attaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAagaaga                                                         >  2:2800483/1‑148 (MQ=255)
                                                                                       ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAt                                                        >  2:348829/1‑148 (MQ=255)
                                                                                       ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAt                                                        <  1:1520306/148‑1 (MQ=255)
                                                                                        taAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAtt                                                       >  1:3379169/1‑148 (MQ=255)
                                                                                           aGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaa                                                    <  1:2759197/148‑1 (MQ=255)
                                                                                            ggCTTATAACACCTTCAGGCGGGCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaaa                                                   <  2:2003970/148‑1 (MQ=255)
                                                                                            ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaaa                                                   >  1:529778/1‑148 (MQ=255)
                                                                                            ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaaa                                                   <  2:2003459/148‑1 (MQ=255)
                                                                                              cTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAAcc                                                 >  1:2161602/1‑148 (MQ=255)
                                                                                                tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCa                                               >  2:1502970/1‑148 (MQ=255)
                                                                                                tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCa                                               <  2:907557/148‑1 (MQ=255)
                                                                                                tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCa                                               >  1:1841325/1‑148 (MQ=255)
                                                                                                tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCa                                               <  1:2030979/148‑1 (MQ=255)
                                                                                                  taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtc                                             <  2:1959984/148‑1 (MQ=255)
                                                                                                  taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtc                                             >  2:323181/1‑148 (MQ=255)
                                                                                                     cacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCg                                          >  1:2601873/1‑148 (MQ=255)
                                                                                                     cacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCg                                          >  1:2844573/1‑148 (MQ=255)
                                                                                                        cTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAg                                       <  1:3416799/148‑1 (MQ=255)
                                                                                                          tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGAc                                     >  2:1858874/1‑148 (MQ=255)
                                                                                                            aGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACAt                                   <  1:4062/148‑1 (MQ=255)
                                                                                                             ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATg                                  <  1:2072013/148‑1 (MQ=255)
                                                                                                             cccggcCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATg                                  >  2:2738254/3‑148 (MQ=255)
                                                                                                              gcggcCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATgg                                 <  2:1220716/148‑1 (MQ=255)
                                                                                                               cggcCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGt                                <  2:2149333/148‑1 (MQ=255)
                                                                                                                               tttCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAAcgcg                >  1:3018239/1‑148 (MQ=255)
                                                                                                                                ttCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACgcgc               >  1:3808671/1‑148 (MQ=255)
                                                                                                                                ttCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACgcgc               >  2:3427884/1‑148 (MQ=255)
                                                                                                                                 tCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCa              >  2:2107410/1‑148 (MQ=255)
                                                                                                                                 tCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCa              >  1:1777315/1‑148 (MQ=255)
                                                                                                                                  ctttttATGGATAATACTTATGAATATATCGCTGGTTGCGTTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCTAGAAGATTGAAACCCTTCTCGAAGACATGGTCGCTCAGGCGAACGCGCaa             <  2:3136192/146‑1 (MQ=255)
                                                                                                                                     tttATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTc          >  2:2144147/1‑148 (MQ=255)
                                                                                                                                             taataaTTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGAc  <  1:2204288/148‑1 (MQ=255)
                                                                                                                                             taataaTTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGAc  >  1:1625487/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGAC  >  CP009273/1683960‑1684248

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: