Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,724,092 T→C A245A (GCA→GCG)  U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,0920TC99.8% 1525.9 / NA 398A245A (GCA→GCG) U712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (194/202);  total (196/202)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 403 total aligned reads displayed.

TCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAAAACTCTGCCAGCTCAAGCGCACATTGATAGGCAAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGACTGACTGCGTATCTCAAATGCGTGACAAGGCGTTGATAGGAAAT  >  CP006881/3723945‑3724238
                                                                                                                                                   |                                                                                                                                                  
tCCTGATGAATAAAAAACCTTCCAGAGATCACTCACTGCAAGTCAAAGCGGTATATCCAAGAGTCCTTGATGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGc                                                                                                                                                  <  1:2059698/150‑1 (MQ=255)
 cctgATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGcc                                                                                                                                                 <  1:1247297/147‑1 (MQ=255)
   tgATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGGGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATGTTTCAAAAACTCCGcccg                                                                                                                                               >  2:1588497/2‑148 (MQ=255)
    gATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGATATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGc                                                                                                                                              >  2:684002/1‑150 (MQ=255)
        aataaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTcca                                                                                                                                          <  1:3000165/150‑3 (MQ=255)
         ataaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAg                                                                                                                                         >  1:1298743/1‑150 (MQ=255)
         ataaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAg                                                                                                                                         <  1:222638/150‑1 (MQ=255)
          taaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAgc                                                                                                                                        <  1:1982844/150‑1 (MQ=255)
              aaaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGcac                                                                                                                                    <  1:337126/150‑1 (MQ=255)
              aaaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGcac                                                                                                                                    >  1:3305690/1‑150 (MQ=255)
               aaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCGGCGTGATCTAGCCGGCCCCGGATTCCGGCAGATGTAATTGCTATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGcacc                                                                                                                                   >  2:145808/1‑149 (MQ=255)
               aaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGcaca                                                                                                                                   <  1:3640669/150‑1 (MQ=255)
                 aCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACAtt                                                                                                                                 <  1:3249784/150‑1 (MQ=255)
                 aCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAAGTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACAtt                                                                                                                                 >  1:863047/1‑150 (MQ=255)
                  cTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTg                                                                                                                                >  1:155221/1‑150 (MQ=255)
                   tcccAGAGATCACCCCCTGCTAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATCCCGGCAAATGTAAATGATATGCATTTTTTAAAAACTCCGCCAGCTCCAGCGCACATTGa                                                                                                                               <  1:516366/148‑1 (MQ=255)
                       agagATCCCTCCCTGCAAGCCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAgg                                                                                                                           <  2:3299508/150‑1 (MQ=255)
                       agagATCACCCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAgg                                                                                                                           <  1:2421814/150‑1 (MQ=255)
                           aTCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCgaa                                                                                                                       >  2:4176540/1‑150 (MQ=255)
                             cACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATg                                                                                                                     >  2:3201327/1‑150 (MQ=255)
                             cACGCCCTGGAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATAAAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATg                                                                                                                     <  2:1787756/150‑1 (MQ=255)
                              aCTCCCTGCAGGTCAAAGCGGTTATTCCGAGAGTTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGa                                                                                                                    <  1:1780304/150‑1 (MQ=255)
                              aCTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGGAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGa                                                                                                                    >  2:3928380/1‑150 (MQ=255)
                              aCTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCCCCAGCTCCAGCGCACATTGATAGGCGAATGa                                                                                                                    <  1:1464361/150‑1 (MQ=255)
                               cTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGaa                                                                                                                   >  1:308470/1‑150 (MQ=255)
                                tCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCgaatgaat                                                                                                                  <  1:2711049/150‑1 (MQ=255)
                                tCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATACAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCgaatgaat                                                                                                                  <  2:615828/150‑1 (MQ=255)
                                 cccTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATa                                                                                                                 >  2:2651733/1‑150 (MQ=255)
                                   cTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACt                                                                                                               <  1:4614071/150‑1 (MQ=255)
                                    tGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACtt                                                                                                              >  2:4488856/1‑150 (MQ=255)
                                    tGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACtt                                                                                                              >  2:238105/1‑150 (MQ=255)
                                    tGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACtt                                                                                                              <  2:2243121/150‑1 (MQ=255)
                                    tGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACtt                                                                                                              >  2:217193/1‑150 (MQ=255)
                                       aaGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACttttt                                                                                                           <  2:3128473/150‑1 (MQ=255)
                                       aaGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACttttt                                                                                                           >  1:828601/1‑150 (MQ=255)
                                         gTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGGTAGGCGAATGAATACttttttt                                                                                                         >  2:1667224/1‑150 (MQ=255)
                                         gTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACttttttt                                                                                                         <  1:4179322/150‑1 (MQ=255)
                                           cAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTgg                                                                                                       <  2:2957794/150‑1 (MQ=255)
                                           cAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTgg                                                                                                       >  2:3975025/1‑150 (MQ=255)
                                            aaaGCGGGTATTCCGAGAGATCTTGAAGACGCTGGACAAGCAGCGTGATATAGCCGGCCTCGGATCCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGa                                                                                                      <  1:3978525/150‑1 (MQ=255)
                                             aaGCGGTTATTCCGAGAGACCTTGAAGCCGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGAt                                                                                                     <  2:2598206/150‑1 (MQ=255)
                                               gCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGa                                                                                                   >  1:2323949/1‑150 (MQ=255)
                                               gCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCACGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGa                                                                                                   <  2:4574859/150‑1 (MQ=255)
                                                 ggTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTGTTTGGATGaaa                                                                                                 >  2:3249983/1‑150 (MQ=255)
                                                  gTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAAt                                                                                                <  1:1377090/150‑1 (MQ=255)
                                                  gTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAAt                                                                                                <  1:119626/150‑1 (MQ=255)
                                                  gTTATCCGGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAAt                                                                                                <  2:3843088/150‑1 (MQ=255)
                                                   ttATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATa                                                                                               <  2:3702938/150‑1 (MQ=255)
                                                   gtATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGAGTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATa                                                                                               <  2:2296344/149‑1 (MQ=255)
                                                     aTTCCGAGAGATCTTGAAGACGCTGGACACGCAGCGTGATATAGCCGGCATGGGAATCCGGCAAATGATATCGAAATCCACATATCAAAAAGACGGCCGGCTGCAGCGCACAGATGTAGGCGAATGAATACTTTAATGTGATGAAAGaat                                                                                             >  2:1286261/1‑149 (MQ=255)
                                                      ttCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAg                                                                                            >  2:3562955/1‑150 (MQ=255)
                                                      ttCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAg                                                                                            >  2:2083681/1‑150 (MQ=255)
                                                       tCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGc                                                                                           >  1:2570021/1‑150 (MQ=255)
                                                         cGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAt                                                                                         <  1:909753/150‑1 (MQ=255)
                                                         cGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAt                                                                                         <  1:31496/150‑1 (MQ=255)
                                                          gagagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGGTATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAtc                                                                                        >  1:657469/1‑150 (MQ=255)
                                                          gagagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAtc                                                                                        <  1:4362708/150‑1 (MQ=255)
                                                           agagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAtct                                                                                       <  2:3360635/150‑1 (MQ=255)
                                                             agaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATctcc                                                                                     <  2:1342808/150‑1 (MQ=255)
                                                              gaTCTTGAAGACGGTGGACATGCAGCGTGATATAGCCGGCCTCGGATCCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATctcct                                                                                    <  2:561705/150‑1 (MQ=255)
                                                                 cTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGt                                                                                 <  1:166869/150‑1 (MQ=255)
                                                                 cTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGt                                                                                 <  2:4629135/150‑1 (MQ=255)
                                                                   tggaGACGCTGGGCATGCAGGGTGATATAGCCGGCCCCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTcc                                                                               <  1:3548562/147‑1 (MQ=255)
                                                                   tGAAGACGCTGGCCATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAGTAAAGCATCTCCTCGTcc                                                                               <  2:4355602/150‑1 (MQ=255)
                                                                   tGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTcc                                                                               >  2:4068576/1‑150 (MQ=255)
                                                                       gACGCTGGACATGCAGCGTGATATAGCCGGCCCCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATAc                                                                           <  2:3081085/150‑1 (MQ=255)
                                                                            tGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGGTGAAATAAAGCATCTCCTCGTCCATACGATgg                                                                      >  2:3501820/1‑150 (MQ=255)
                                                                             ggACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGa                                                                     <  2:2027301/150‑1 (MQ=255)
                                                                                cATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAgcg                                                                  >  2:1728350/1‑150 (MQ=255)
                                                                                 aTGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAgcgc                                                                 <  1:3688758/150‑1 (MQ=255)
                                                                                    cAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTc                                                              >  2:2457007/1‑150 (MQ=255)
                                                                                     agagTGATATAACCGCCCTTGGATCCCGGCAAATGGAATTGATATTCATTTTTCAAAAACTCTGCCAGCTCCAGCTCACATTGAGTGGCGAATGAAGACTTTTTTTGGATGTAATAAAGCATCTCCTCTACCACACGATGGTTCGCTTct                                                             >  1:4080819/4‑149 (MQ=255)
                                                                                     aGCGTGATAAAGCCGGCCTCGGATTCCGGCAAATGTAATTCATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGATCGCTTCa                                                             <  1:993491/150‑1 (MQ=255)
                                                                                      gCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCAt                                                            <  2:1060984/150‑1 (MQ=255)
                                                                                          gATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGa                                                        >  2:1173933/1‑150 (MQ=255)
                                                                                          gATATAGCCGGCCTCGGATTCCGGCAAATGTAAGAGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCGCATTGATAGGCGAATTAATACTTGTTTTGGATGAAGTAAAGGATCTCCTCGTCCATACGATGGAGCGCTTCATTTGa                                                        <  2:3456079/150‑1 (MQ=255)
                                                                                           atatACCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGAt                                                       >  1:1159839/1‑150 (MQ=255)
                                                                                              taGCCGGCCTCGGATTCCGGCAAATGTGATTGACATCCATTTTTCAAAACCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATGCTTTTTTTGGATGAAATAAGGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTcc                                                    <  1:2544805/150‑1 (MQ=255)
                                                                                               aGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCa                                                   >  2:4279570/1‑150 (MQ=255)
                                                                                                gCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCaa                                                  <  2:1957939/150‑1 (MQ=255)
                                                                                                   ggCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCg                                               <  2:4768203/150‑1 (MQ=255)
                                                                                                         ggATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCggctgac                                         >  2:1110417/1‑150 (MQ=255)
                                                                                                         ggATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCggctgac                                         >  1:1102937/1‑150 (MQ=255)
                                                                                                                  cacaTGTAATTGCTACTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAtcg                                <  1:4398555/147‑2 (MQ=255)
                                                                                                                  cAAATGTAATTGATATTCATTTTCCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAtcg                                <  2:4132518/150‑2 (MQ=255)
                                                                                                                   aaaTGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAtcgc                               <  1:1455121/150‑3 (MQ=255)
                                                                                                                    aaTGTAATTGATATTCATTTTCCAAAAACTCCGCCCGCCCCAGGGCACAGTGATAGGCGAGTGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCa                              <  1:1019782/150‑1 (MQ=255)
                                                                                                                     aTGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGGGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCaa                             >  2:2508454/1‑150 (MQ=255)
                                                                                                                     aTGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCaa                             >  2:1021964/1‑150 (MQ=255)
                                                                                                                        tAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATg                          <  2:89781/150‑1 (MQ=255)
                                                                                                                        tAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATg                          >  1:1923606/1‑150 (MQ=255)
                                                                                                                         aaTTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGc                         >  1:4537027/1‑150 (MQ=255)
                                                                                                                                   catttttCAAAAACTCCGACAGCTCCAGCGCACGTCGATAGGCGAATCAATACTGTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGc               >  1:3731638/4‑150 (MQ=255)
                                                                                                                                        ttCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAGTAAAGCATCCCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTgg          >  2:2811428/1‑149 (MQ=255)
                                                                                                                                        ttCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGGGACAAGGCGTTGa          >  2:3620792/1‑150 (MQ=255)
                                                                                                                                         tCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGAt         >  2:2675341/1‑150 (MQ=255)
                                                                                                                                           aaaaaCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAg       >  2:1865764/1‑150 (MQ=255)
                                                                                                                                             aaaCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGa     >  2:3798112/1‑150 (MQ=255)
                                                                                                                                              aaCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGaa    >  2:1917189/1‑150 (MQ=255)
                                                                                                                                              aaCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGaa    <  2:141767/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAAAACTCTGCCAGCTCAAGCGCACATTGATAGGCAAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGACTGACTGCGTATCTCAAATGCGTGACAAGGCGTTGATAGGAAAT  >  CP006881/3723945‑3724238

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