Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,919,922 G→T 100% D86Y (GAT→TAT)  ygdH → UPF0717 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,919,9220GT99.2% 799.6 / ‑3.5 247D86Y (GAT→TAT) ygdHUPF0717 family protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (112/133);  total (114/133)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-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 248 total aligned reads displayed.

TTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTT  >  CP009273/2919777‑2920066
                                                                                                                                                 |                                                                                                                                                
ttCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTAt                                                                                                                                                <  1:3231492/148‑1 (MQ=255)
 tCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATg                                                                                                                                               <  2:2446096/148‑1 (MQ=255)
     ttGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCg                                                                                                                                           >  2:1256792/1‑148 (MQ=255)
      tGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGa                                                                                                                                          <  1:3336382/148‑1 (MQ=255)
       gCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGaa                                                                                                                                         <  2:2637613/148‑1 (MQ=255)
        ccGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAat                                                                                                                                        <  2:1089788/148‑1 (MQ=255)
        ccGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAat                                                                                                                                        <  2:2068953/148‑1 (MQ=255)
         cGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAatt                                                                                                                                       <  1:1386348/148‑1 (MQ=255)
         cGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAatt                                                                                                                                       >  2:2054497/1‑148 (MQ=255)
          gTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAatta                                                                                                                                      <  2:10595/148‑1 (MQ=255)
            aCTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAattatt                                                                                                                                    >  2:2981019/1‑148 (MQ=255)
            aCTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAattatt                                                                                                                                    >  2:2138417/1‑148 (MQ=255)
              tGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTcg                                                                                                                                  <  2:3474870/148‑1 (MQ=255)
               gAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTcgc                                                                                                                                 >  2:1462139/1‑148 (MQ=255)
                  cTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCt                                                                                                                              <  2:2041508/148‑1 (MQ=255)
                  cTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCt                                                                                                                              >  2:2104169/1‑148 (MQ=255)
                          gTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGc                                                                                                                      <  2:1182318/148‑1 (MQ=255)
                            ttGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCa                                                                                                                    >  2:816894/1‑148 (MQ=255)
                              gACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAAc                                                                                                                  <  1:2579390/148‑1 (MQ=255)
                               aCCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACt                                                                                                                 >  1:1553212/1‑148 (MQ=255)
                                 cGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTg                                                                                                               >  2:2674987/1‑148 (MQ=255)
                                 cGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTg                                                                                                               <  1:655638/148‑1 (MQ=255)
                                    tAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTc                                                                                                            >  2:2725488/1‑148 (MQ=255)
                                    tAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTc                                                                                                            >  1:1253425/1‑148 (MQ=255)
                                     aaCAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTcg                                                                                                           >  1:3145528/1‑148 (MQ=255)
                                     aaCAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTcg                                                                                                           <  1:3207211/148‑1 (MQ=255)
                                        aGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCgg                                                                                                        <  2:2941195/148‑1 (MQ=255)
                                         gCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGt                                                                                                       >  1:2312702/1‑148 (MQ=255)
                                         gCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGt                                                                                                       >  1:492927/1‑148 (MQ=255)
                                          cAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTc                                                                                                      >  2:2365880/1‑148 (MQ=255)
                                          cAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTc                                                                                                      >  2:248956/1‑148 (MQ=255)
                                          cAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTc                                                                                                      >  1:731382/1‑148 (MQ=255)
                                          cAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTc                                                                                                      <  2:546042/148‑1 (MQ=255)
                                          cAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTc                                                                                                      <  1:1064315/148‑1 (MQ=255)
                                          cAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTc                                                                                                      >  1:3221410/1‑148 (MQ=255)
                                           aaaGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTcc                                                                                                     <  1:180692/148‑1 (MQ=255)
                                            aaGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCTCGCTTTGCAGGCCAACTTGTTCGCGGTCCt                                                                                                    >  1:1088906/1‑148 (MQ=255)
                                            aaGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCt                                                                                                    >  1:2034989/1‑148 (MQ=255)
                                            aaGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCt                                                                                                    >  1:2797558/1‑148 (MQ=255)
                                            aaGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCt                                                                                                    >  1:1044887/1‑148 (MQ=255)
                                             aGAATTGCTGTCTCGTTTTGAAAATTTTGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTg                                                                                                   >  2:3159710/1‑148 (MQ=255)
                                             aGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTg                                                                                                   <  2:131210/148‑1 (MQ=255)
                                               aagtGCTGTCTCGTTTTGAAAATTTCGATCTTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCg                                                                                                 <  1:3098961/145‑1 (MQ=255)
                                               aaTTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCg                                                                                                 <  1:3098979/148‑1 (MQ=255)
                                                 ttGCTGTCTCGTTTTGAAAATTTCGATATTCACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTg                                                                                               <  1:912892/148‑1 (MQ=255)
                                                 gtGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTg                                                                                               <  2:2928031/147‑1 (MQ=255)
                                                  tGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGa                                                                                              >  2:235492/1‑148 (MQ=255)
                                                  tGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGa                                                                                              >  2:236225/1‑148 (MQ=255)
                                                     tGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACAt                                                                                           <  2:2509062/148‑1 (MQ=255)
                                                      gTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACAtt                                                                                          >  1:2658330/1‑148 (MQ=255)
                                                        ctcGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATtct                                                                                        <  1:2781179/148‑1 (MQ=255)
                                                        ctcGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATtct                                                                                        <  2:1932189/148‑1 (MQ=255)
                                                          cGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTctct                                                                                      <  2:431593/148‑1 (MQ=255)
                                                          cGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTctct                                                                                      <  1:2007234/148‑1 (MQ=255)
                                                             tttGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCg                                                                                   >  1:1751918/1‑148 (MQ=255)
                                                             tttGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCg                                                                                   >  2:612224/1‑148 (MQ=255)
                                                             tttGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCg                                                                                   >  1:2431995/1‑148 (MQ=255)
                                                               tGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGtt                                                                                 <  1:2644903/148‑1 (MQ=255)
                                                               tGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGtt                                                                                 <  2:2339724/148‑1 (MQ=255)
                                                                 aaaaTTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTa                                                                               >  1:594223/1‑148 (MQ=255)
                                                                 aaaaTTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTa                                                                               >  2:1863089/1‑148 (MQ=255)
                                                                      ttCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGc                                                                          >  2:3391090/1‑148 (MQ=255)
                                                                        cGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCaa                                                                        >  2:3484963/1‑148 (MQ=255)
                                                                         gATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCaaa                                                                       <  2:3078619/148‑1 (MQ=255)
                                                                         gATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCaaa                                                                       <  1:546405/148‑1 (MQ=255)
                                                                          atatTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAAt                                                                      <  2:65129/148‑1 (MQ=255)
                                                                           tatTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATc                                                                     <  2:921736/148‑1 (MQ=255)
                                                                             ttAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCa                                                                   >  2:1140331/1‑148 (MQ=255)
                                                                                  gTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAaca                                                              >  2:1695883/1‑148 (MQ=255)
                                                                                     ttGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACcg                                                           <  1:1195694/148‑1 (MQ=255)
                                                                                      tGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACcgt                                                          <  2:2682137/148‑1 (MQ=255)
                                                                                         gcCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACcgttcg                                                       <  2:150686/148‑1 (MQ=255)
                                                                                          ccGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACcgttcgt                                                      <  1:2582839/148‑1 (MQ=255)
                                                                                          ccGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACcgttcgt                                                      >  1:2175127/1‑148 (MQ=255)
                                                                                          ccGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACcgttcgt                                                      <  1:959806/148‑1 (MQ=255)
                                                                                          ccGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTGTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACcgttcgt                                                      >  2:3196688/1‑148 (MQ=255)
                                                                                              gAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATACGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGAGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACAGTACGTTTTc                                                  >  1:1474470/1‑148 (MQ=255)
                                                                                                 cgcgGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCa                                               >  1:389776/1‑148 (MQ=255)
                                                                                                     gcgTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCt                                           <  1:78459/148‑1 (MQ=255)
                                                                                                       gTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGa                                         <  2:444917/148‑1 (MQ=255)
                                                                                                         aaaGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAAt                                       >  1:2209467/1‑148 (MQ=255)
                                                                                                           aGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAAtct                                     >  1:449222/1‑148 (MQ=255)
                                                                                                           aGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAAtct                                     >  1:2279473/1‑148 (MQ=255)
                                                                                                             cTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCg                                   >  2:1392234/1‑148 (MQ=255)
                                                                                                             cTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCg                                   >  2:598988/1‑148 (MQ=255)
                                                                                                              tGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGa                                  <  1:2995139/148‑1 (MQ=255)
                                                                                                                    tGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTc                            <  1:2154084/148‑1 (MQ=255)
                                                                                                                    tGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTc                            <  1:573607/148‑1 (MQ=255)
                                                                                                                        tAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTc                        <  2:2189131/148‑1 (MQ=255)
                                                                                                                           tCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCcac                     >  1:3487831/1‑148 (MQ=255)
                                                                                                                           tCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCcac                     >  2:2650781/1‑148 (MQ=255)
                                                                                                                           tCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCcac                     <  1:2984115/148‑1 (MQ=255)
                                                                                                                            cccccGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCcaca                    >  2:2841414/1‑148 (MQ=255)
                                                                                                                                 ggAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACt               >  2:3369191/1‑148 (MQ=255)
                                                                                                                                 ggAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACt               >  1:413955/1‑148 (MQ=255)
                                                                                                                                    agagGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACAtaactaac            <  1:612224/148‑1 (MQ=255)
                                                                                                                                        gCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTgg        >  1:3013045/1‑148 (MQ=255)
                                                                                                                                         cTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGt       >  2:2186278/1‑148 (MQ=255)
                                                                                                                                             tGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCttt   >  1:141011/1‑148 (MQ=255)
                                                                                                                                             tGTCTATGGGCGAATTAGTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCttt   <  2:1802399/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTT  >  CP009273/2919777‑2920066

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: