Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,724,200 T→C A209A (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,2000TC100.0% 1127.3 / NA 288A209A (GCA→GCGU712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (127/161);  total (127/161)

Only 100 of 295 total aligned reads displayed.

CTCCGGCAAATGTAATTGATATTCACTTTTCAAAAACTCTGCCAGCTCAAGCGCACATTGATAGGCAAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGACTGACTGCGTATCTCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTCTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAATACATGCTCTCCGCATCCATCT  >  CP006881/3724053‑3724345
                                                                                                                                                   |                                                                                                                                                 
ttCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCACTTGATTCCCACCGGCTGACCGc                                                                                                                                                 >  2:2857017/2‑150 (MQ=255)
 tCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGGGCGCCTCATTTGATTCCAACCGGCTGGCcgcg                                                                                                                                                >  2:2824339/1‑150 (MQ=255)
  ccGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGt                                                                                                                                               <  1:2005015/150‑1 (MQ=255)
   cGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTa                                                                                                                                              <  1:3334887/150‑1 (MQ=255)
   cGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTa                                                                                                                                              >  2:1248996/1‑150 (MQ=255)
    ggCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAt                                                                                                                                             >  1:1126634/1‑150 (MQ=255)
     gCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAtc                                                                                                                                            >  2:2600888/1‑150 (MQ=255)
     gCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAtc                                                                                                                                            <  2:2045071/150‑1 (MQ=255)
       aaaTGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAtcgc                                                                                                                                          >  2:1010134/1‑148 (MQ=255)
         aTGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCaa                                                                                                                                        <  2:1053037/150‑1 (MQ=255)
          tGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCaaa                                                                                                                                       >  2:3911032/1‑150 (MQ=255)
          tGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCCCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCaga                                                                                                                                       >  1:2078091/1‑148 (MQ=255)
            tAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATg                                                                                                                                     <  1:2777460/150‑1 (MQ=255)
             aaTTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGc                                                                                                                                    <  2:446545/150‑1 (MQ=255)
              aTTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCg                                                                                                                                   <  2:2940698/150‑1 (MQ=255)
               ttGATATTCATGTTTCAAAAACTCCGCCAGCGCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGt                                                                                                                                  <  2:3506408/150‑1 (MQ=255)
                 gATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGa                                                                                                                                >  1:2894204/1‑150 (MQ=255)
                  atatTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAGGCATCTCCTCGTCCATACGATGGGGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGGGAc                                                                                                                               >  2:968423/1‑150 (MQ=255)
                      tCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAgg                                                                                                                           <  1:3789250/150‑1 (MQ=255)
                       catttttCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGc                                                                                                                          <  2:1791450/147‑1 (MQ=255)
                        atttttCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCCTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGGCCAGGCg                                                                                                                         >  2:1324898/3‑150 (MQ=255)
                          ttttCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTAGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGtt                                                                                                                       <  2:2427812/150‑1 (MQ=255)
                            ttCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGa                                                                                                                     >  2:2732928/1‑150 (MQ=255)
                             tCAAAAACTCCGCCAGCTCCAGGGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGAt                                                                                                                    <  1:3492913/150‑1 (MQ=255)
                              cAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATa                                                                                                                   >  1:1664701/1‑150 (MQ=255)
                              cAAAAACTCCGCCACCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATa                                                                                                                   <  2:167890/150‑1 (MQ=255)
                                aaaaCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAgg                                                                                                                 <  2:2567057/150‑1 (MQ=255)
                                aaaaCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGGGGCAAGGCGTTGATAgg                                                                                                                 >  2:1592561/1‑150 (MQ=255)
                                 aaaCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACCATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGa                                                                                                                >  2:363169/1‑150 (MQ=255)
                                  aaCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGTCCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGaa                                                                                                               >  1:3800880/1‑150 (MQ=255)
                                  aaCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGGCAAGGCGTTGATAGGaa                                                                                                               >  2:2351254/1‑150 (MQ=255)
                                   aCTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGaaa                                                                                                              >  2:2074464/1‑150 (MQ=255)
                                     tccgCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATg                                                                                                            <  2:3800334/147‑1 (MQ=255)
                                     cccgCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTGTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATg                                                                                                            <  2:1122447/147‑1 (MQ=255)
                                      ccgCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGc                                                                                                           <  1:2496121/148‑1 (MQ=255)
                                        gCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTg                                                                                                         >  2:1731582/1‑150 (MQ=255)
                                        gCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTGGGTAGGAAATGCTg                                                                                                         >  2:2468565/1‑150 (MQ=255)
                                         ccAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGt                                                                                                        <  2:78422/150‑1 (MQ=255)
                                            gCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGtttt                                                                                                     <  1:3403444/150‑1 (MQ=255)
                                            gCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGtttt                                                                                                     <  1:2543784/150‑1 (MQ=255)
                                             cTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGGTAGGGAATGCTGTTTTc                                                                                                    >  2:3599293/1‑150 (MQ=255)
                                             cTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTc                                                                                                    <  2:3066880/150‑1 (MQ=255)
                                             cTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTc                                                                                                    <  1:724244/150‑1 (MQ=255)
                                              cccaGCGCCCATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGAATCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCa                                                                                                   <  1:2221944/147‑1 (MQ=255)
                                                   cgcACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGGGCTTCATTTGATACCAACCGGCTGACCGCGTATCGCAAATGTGTGAGTAGTCGATCATAGCAAATGGTGTTTTCATCCat                                                                                              <  1:1736343/150‑2 (MQ=255)
                                                   cgcACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCAc                                                                                              <  1:2953493/150‑1 (MQ=255)
                                                   cgcACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCAc                                                                                              >  2:917278/1‑150 (MQ=255)
                                                     cacaTTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCt                                                                                            <  2:423883/150‑1 (MQ=255)
                                                      acaTTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCtt                                                                                           >  1:1009595/1‑150 (MQ=255)
                                                        aTTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTc                                                                                         >  2:405654/1‑150 (MQ=255)
                                                         ttGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCg                                                                                        <  1:3740026/150‑1 (MQ=255)
                                                            aTAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGAtt                                                                                     <  2:3816095/150‑1 (MQ=255)
                                                             tAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTa                                                                                    <  2:3490382/150‑1 (MQ=255)
                                                             tAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTa                                                                                    <  2:1367938/150‑1 (MQ=255)
                                                                  gaaTGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTa                                                                               >  2:3859884/2‑150 (MQ=255)
                                                                  gaaTGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTa                                                                               >  2:120930/2‑150 (MQ=255)
                                                                   aaTGAATACTTTTTTTTGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTAt                                                                              <  1:1478325/150‑1 (MQ=255)
                                                                    aTGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGAAtt                                                                             >  1:794538/1‑150 (MQ=255)
                                                                      gAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGt                                                                           <  2:828996/150‑1 (MQ=255)
                                                                          aCTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGACAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTa                                                                       <  1:116563/150‑1 (MQ=255)
                                                                           ctgtttttGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTAt                                                                      <  1:1929181/147‑1 (MQ=255)
                                                                           cTTTTTTTGGATGAAATAAAGCACCTCCTCGTCCATCCGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTAt                                                                      <  1:3946715/150‑1 (MQ=255)
                                                                              tttttGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTAtttt                                                                   >  1:3235220/1‑150 (MQ=255)
                                                                                tttGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTtct                                                                 <  2:1510924/150‑1 (MQ=255)
                                                                                 ttGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTtctc                                                                <  1:580280/150‑1 (MQ=255)
                                                                                  tGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTctct                                                               >  2:671914/1‑150 (MQ=255)
                                                                                  tGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATCCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTctct                                                               <  2:2289396/150‑1 (MQ=255)
                                                                                    gATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTAt                                                             >  2:1482359/1‑150 (MQ=255)
                                                                                     aTGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATc                                                            <  2:1893119/150‑1 (MQ=255)
                                                                                        aaaTAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCAtt                                                         <  1:1270643/150‑1 (MQ=255)
                                                                                        aaaTAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCAtt                                                         >  1:2575476/1‑150 (MQ=255)
                                                                                         aaTAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATCCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCAttt                                                        <  2:3289578/150‑1 (MQ=255)
                                                                                             aaGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTAtcatttcatt                                                    <  1:923136/150‑1 (MQ=255)
                                                                                              aGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGGCCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTg                                                   >  2:2765751/1‑150 (MQ=255)
                                                                                               gCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGa                                                  <  1:1041119/150‑1 (MQ=255)
                                                                                                   ctcctcGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCAt                                              <  1:3152742/150‑1 (MQ=255)
                                                                                                   ctcctcGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCAt                                              <  2:621301/150‑1 (MQ=255)
                                                                                                    tcctcGTCCATCCGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATg                                             <  1:2360826/150‑1 (MQ=255)
                                                                                                      ctcGACCATACGATGGAGTGCTTCATTTGATTCCACCGGGCTGAGCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATgaa                                           <  1:2151595/150‑3 (MQ=255)
                                                                                                        cGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGTTCATGGTCg                                         >  2:3680334/1‑150 (MQ=255)
                                                                                                        cGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCg                                         <  2:373314/150‑1 (MQ=255)
                                                                                                        cGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCg                                         >  2:557438/1‑150 (MQ=255)
                                                                                                         gTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGt                                        >  2:2113685/1‑150 (MQ=255)
                                                                                                                 gATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCa                                <  2:3133509/150‑1 (MQ=255)
                                                                                                                  aTGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAg                               <  1:2435686/150‑1 (MQ=255)
                                                                                                                      aGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCATATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGACGTATGCTTCAGCGCt                           <  1:2611990/150‑1 (MQ=255)
                                                                                                                         gcTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCTTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTGTTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGaa                        <  1:1855421/150‑1 (MQ=255)
                                                                                                                         gcTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGaa                        <  1:363169/150‑1 (MQ=255)
                                                                                                                          cTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAAt                       >  2:1271124/1‑150 (MQ=255)
                                                                                                                          cTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAAt                       <  1:14898/150‑1 (MQ=255)
                                                                                                                          cTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAAt                       >  2:3932313/1‑150 (MQ=255)
                                                                                                                            tCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAATAc                     <  1:3042557/150‑1 (MQ=255)
                                                                                                                                  gATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAATACATGctc               <  1:405654/150‑1 (MQ=255)
                                                                                                                                   aTTCCAACCGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAATACATGctct              >  2:740599/1‑150 (MQ=255)
                                                                                                                                          ccGGCTGACCGCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAATACATGCTCTCCGcatc       <  2:1856757/150‑1 (MQ=255)
                                                                                                                                           cggctgaccgCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAATACATGCTCTCCGcatcc      <  2:2169489/147‑1 (MQ=255)
                                                                                                                                              ctgaccgCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAATACATGCTCTCCGcatccatc   <  1:248509/150‑1 (MQ=255)
                                                                                                                                              ctgaccgCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGGCGTATGCTTCAGCGCTGCAGACATGCTCTCCGcatccatc   >  1:201915/1‑150 (MQ=255)
                                                                                                                                              ctgaccgCGTACCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTATCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAATACATGCTCTCCGcatccatc   <  1:1310956/150‑1 (MQ=255)
                                                                                                                                               tgaccgCGTATCGCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTTTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAATACATGCTCTCCGCATCCATCt  <  1:1248996/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTCCGGCAAATGTAATTGATATTCACTTTTCAAAAACTCTGCCAGCTCAAGCGCACATTGATAGGCAAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGACTGACTGCGTATCTCAAATGCGTGACAAGGCGTTGATAGGAAATGCTGTTTTCATCCACCTTTCGATTAAAGTATTGTTCTATTTTCTCTATCATTTCTTTGATCATGGTCGTATGCTTCAGCGCTGAATACATGCTCTCCGCATCCATCT  >  CP006881/3724053‑3724345

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