Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,724,078 C→T S250N (AGT→AAT)  U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,0780CT99.2% 1394.4 / NA 401S250N (AGT→AAT) U712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (191/205);  total (195/206)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 407 total aligned reads displayed.

TATTACACGCTCAATCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAAAACTCTGCCAGCTCAAGCGCACATTGATAGGCAAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGACTGACTGCGTATCTCAAATGCGTGACAAGG  >  CP006881/3723931‑3724224
                                                                                                                                                   |                                                                                                                                                  
tactACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAttt                                                                                                                                                  <  2:3392560/147‑1 (MQ=255)
  ctACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATACCGGCAAATGTAATTGATATTCATtttt                                                                                                                                                <  1:4717394/149‑1 (MQ=255)
     cacaCTGAATCCTGAACAATAAAAAACCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaa                                                                                                                                             <  1:2232928/150‑1 (MQ=255)
      acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGTCGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGGAAATGAAATTCAGTTTTCaaa                                                                                                                                            >  2:2258226/7‑150 (MQ=255)
      acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaa                                                                                                                                            >  1:2726462/7‑150 (MQ=255)
      acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaa                                                                                                                                            >  2:1114192/7‑150 (MQ=255)
       cactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaaa                                                                                                                                           >  2:2457174/6‑150 (MQ=255)
        actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGCCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaaaa                                                                                                                                          <  2:3328424/146‑1 (MQ=255)
         ctgaaTCCTGATCAATAAACAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAAc                                                                                                                                         <  2:3294739/147‑1 (MQ=255)
          tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAAct                                                                                                                                        <  2:1469522/148‑1 (MQ=255)
           gaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAActc                                                                                                                                       >  2:2547172/2‑150 (MQ=255)
           gaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAActc                                                                                                                                       >  1:315898/2‑150 (MQ=255)
            aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAActcc                                                                                                                                      <  1:3056024/150‑2 (MQ=255)
             aTCCTGATCAATAAACAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAActccg                                                                                                                                     <  2:250348/150‑3 (MQ=255)
              tCCTGATGAATAAAAAACCTTCCAGAGATCACTCACTGCAAGTCAAAGCGGTATATCCAAGAGTCCTTGATGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGc                                                                                                                                    <  1:2059698/150‑1 (MQ=255)
               cctgATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGcc                                                                                                                                   >  2:3457740/4‑150 (MQ=255)
                ctgATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCa                                                                                                                                  <  2:550464/148‑1 (MQ=255)
                 tgATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGGGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATGTTTCAAAAACTCCGcccg                                                                                                                                 >  2:1588497/2‑148 (MQ=255)
                 tgATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAg                                                                                                                                 >  2:3536798/2‑150 (MQ=255)
                  gATCAATAAAAACACTTCCAGAGATCACTCCCTGCAAGTCAAGGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGc                                                                                                                                <  1:1137514/150‑1 (MQ=255)
                  gATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGc                                                                                                                                >  1:1920196/1‑150 (MQ=255)
                    tcaataaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTc                                                                                                                              >  2:1436527/3‑150 (MQ=255)
                     caataaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTcc                                                                                                                             >  1:4442357/2‑149 (MQ=255)
                      aataaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTcca                                                                                                                            <  1:3000165/150‑3 (MQ=255)
                      aataaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTcca                                                                                                                            <  2:1160323/150‑3 (MQ=255)
                       ataaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAg                                                                                                                           >  1:1298743/1‑150 (MQ=255)
                        taaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAgc                                                                                                                          <  1:1965384/150‑1 (MQ=255)
                        taaaAAAACTTCCAGAGATCACTCCCTGCAAGACAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAgc                                                                                                                          <  1:981180/150‑1 (MQ=255)
                          aaaaaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAgcgc                                                                                                                        >  2:1299568/1‑150 (MQ=255)
                           aaaaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGca                                                                                                                       >  2:2787810/1‑150 (MQ=255)
                           aaaaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGca                                                                                                                       <  2:3048880/150‑1 (MQ=255)
                            aaaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGcac                                                                                                                      >  1:3305690/1‑150 (MQ=255)
                            aaaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGcac                                                                                                                      <  1:337126/150‑1 (MQ=255)
                                cTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTg                                                                                                                  <  1:1594262/150‑1 (MQ=255)
                                   ccAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGGTTCCGGCAAATGTAATTGATATTCATTTTCCAAAAACTCCGCCAGCTCCAGCGCACATTGATa                                                                                                               >  2:1073149/1‑150 (MQ=255)
                                     agagATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAgg                                                                                                             >  2:2197596/1‑150 (MQ=255)
                                     agagATCACCCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAgg                                                                                                             <  1:2421814/150‑1 (MQ=255)
                                       agaTCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTGTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGcg                                                                                                           >  2:2905580/1‑149 (MQ=255)
                                         aTCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCgaa                                                                                                         >  2:4176540/1‑150 (MQ=255)
                                          tCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCATAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAAt                                                                                                        >  1:903112/1‑150 (MQ=255)
                                           cACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATg                                                                                                       >  2:3201327/1‑150 (MQ=255)
                                           cACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATg                                                                                                       <  2:677085/150‑1 (MQ=255)
                                            aCTCCCTGCAGGTCAAAGCGGTTATTCCGAGAGTTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGa                                                                                                      <  1:1780304/150‑1 (MQ=255)
                                            aCTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGa                                                                                                      <  1:758495/150‑1 (MQ=255)
                                              tCCCTGCAGGTCAAGGCGGTCATTCGGAGAGATCTTGGAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCgaatgaat                                                                                                    <  1:4241136/150‑1 (MQ=255)
                                              tCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCgaatgaat                                                                                                    <  2:2037585/150‑1 (MQ=255)
                                              tCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCgaatgaat                                                                                                    >  2:389467/1‑150 (MQ=255)
                                                ccTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATAc                                                                                                  <  1:4377462/150‑1 (MQ=255)
                                                 cTGCAAGACAATGCGGTTATTCCGAGAGATCTTGACGACGCTGGACATGCAGCGAGATACAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGGTCCAGCGCACATTGATAGGCGAATGAATACt                                                                                                 <  1:3709444/150‑1 (MQ=255)
                                                  tGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACtt                                                                                                >  2:4488856/1‑150 (MQ=255)
                                                  tGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACtt                                                                                                >  2:238105/1‑150 (MQ=255)
                                                  tGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACtt                                                                                                <  1:3572281/150‑1 (MQ=255)
                                                        tCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTg                                                                                          >  1:4037076/1‑150 (MQ=255)
                                                         cAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTgg                                                                                         >  2:3975025/1‑150 (MQ=255)
                                                             gCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGa                                                                                     >  2:2663542/1‑150 (MQ=255)
                                                             gCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGa                                                                                     >  1:2323949/1‑150 (MQ=255)
                                                             gCGGTTATCCCGAGAGATCTTGGAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGa                                                                                     <  1:849176/150‑1 (MQ=255)
                                                               ggTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGaaa                                                                                   >  2:1726807/1‑150 (MQ=255)
                                                               ggTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGaaa                                                                                   >  1:3255674/1‑150 (MQ=255)
                                                               ggTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGaaa                                                                                   <  2:4308233/150‑1 (MQ=255)
                                                               ggTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGaaa                                                                                   >  2:438918/1‑150 (MQ=255)
                                                               ggTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGaaa                                                                                   >  1:4558248/1‑150 (MQ=255)
                                                               ggTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTGTTTGGATGaaa                                                                                   >  2:3249983/1‑150 (MQ=255)
                                                                gTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAAt                                                                                  <  1:1377090/150‑1 (MQ=255)
                                                                 gtATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGAGTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATa                                                                                 <  2:2296344/149‑1 (MQ=255)
                                                                   aTTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATaaa                                                                               <  2:4631822/150‑1 (MQ=255)
                                                                     tCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGc                                                                             >  1:2570021/1‑150 (MQ=255)
                                                                       cGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAt                                                                           >  2:1765603/1‑150 (MQ=255)
                                                                       cGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAt                                                                           <  2:3692631/150‑1 (MQ=255)
                                                                       cGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAt                                                                           >  2:1238631/1‑150 (MQ=255)
                                                                        gagagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAtc                                                                          <  1:1101444/150‑1 (MQ=255)
                                                                        gagagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAtc                                                                          <  2:1988127/150‑1 (MQ=255)
                                                                           agaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATctcc                                                                       <  1:4148697/150‑1 (MQ=255)
                                                                              tCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCg                                                                    >  2:3833019/1‑150 (MQ=255)
                                                                               cTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGt                                                                   <  1:166869/150‑1 (MQ=255)
                                                                                 tGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTcc                                                                 >  2:4068576/1‑150 (MQ=255)
                                                                                    aGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATa                                                              <  1:3643168/150‑1 (MQ=255)
                                                                                            gACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAg                                                      <  2:1345699/150‑1 (MQ=255)
                                                                                                      gTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGACATTCATTTTTCAACAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAAAACAGCATCTCCTCCTCCACACGACGGAGCGCTTCAttt                                            >  2:1921543/1‑150 (MQ=255)
                                                                                                         atatAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCGTCATTTGAt                                         >  2:612798/1‑150 (MQ=255)
                                                                                                             aGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCa                                     <  2:1803056/150‑1 (MQ=255)
                                                                                                               ccGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAAc                                   <  2:3730472/150‑1 (MQ=255)
                                                                                                               ccGGCCACGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAAc                                   <  1:1187284/150‑1 (MQ=255)
                                                                                                                 ggCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCg                                 <  2:4768203/150‑1 (MQ=255)
                                                                                                                  gCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCgg                                >  2:4134642/1‑149 (MQ=255)
                                                                                                                         attCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCggctgaccg                         <  1:1490919/148‑3 (MQ=255)
                                                                                                                         attCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCggctgaccg                         >  2:3125890/3‑148 (MQ=255)
                                                                                                                                 aaaTGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAtcgc                 <  1:1455121/150‑3 (MQ=255)
                                                                                                                                 aaaTGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAtcgc                 <  1:1217267/150‑3 (MQ=255)
                                                                                                                                  aaTGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCa                <  1:3125890/150‑1 (MQ=255)
                                                                                                                                  aaTGTAATTGATATTCATTTTCCAAAAACTCCGCCCGCCCCAGGGCACAGTGATAGGCGAGTGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCa                <  1:1019782/150‑1 (MQ=255)
                                                                                                                                   aTGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCaa               >  2:1021964/1‑150 (MQ=255)
                                                                                                                                   aTGTAATTGATATTCATTGTTCAAAAACTCCGCCAGCTCCAGCACACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAGAGCATCTCCTCTACCATACGGTGGAGTGCTTCATTTGGTTCCATCCGGCAGACCGCGTAGCGCaa               >  1:2542535/1‑150 (MQ=255)
                                                                                                                                    tGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCaaa              >  2:2141268/1‑150 (MQ=255)
                                                                                                                                    tGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCaaa              <  2:4593475/150‑1 (MQ=255)
                                                                                                                                      tAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATg            <  2:89781/150‑1 (MQ=255)
                                                                                                                                       aaTTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGc           >  1:3956030/1‑150 (MQ=255)
                                                                                                                                        aTTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCg          <  1:2647731/150‑1 (MQ=255)
                                                                                                                                           gATATCCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGa       <  1:533678/150‑1 (MQ=255)
                                                                                                                                            atatTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGGGAc      >  2:2633446/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TATTACACGCTCAATCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAAAACTCTGCCAGCTCAAGCGCACATTGATAGGCAAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGACTGACTGCGTATCTCAAATGCGTGACAAGG  >  CP006881/3723931‑3724224

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑