Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,512,114 C→A 6.5% E104* (GAA→TAA)  aroK ← shikimate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,512,1140CA 6.5% 634.5 / 2.1 243E104* (GAA→TAA) aroKshikimate kinase I
Reads supporting (aligned to +/- strand):  ref base C (87/125);  new base A (14/14);  total (104/139)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

ACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCT  >  CP009273/3511965‑3512262
                                                                                                                                                     |                                                                                                                                                    
aCGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTc                                                                                                                                                      <  2:1500017/150‑1 (MQ=255)
aCGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTc                                                                                                                                                      <  2:3865581/150‑1 (MQ=255)
  gAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCaa                                                                                                                                                    <  2:2910384/150‑1 (MQ=255)
  gAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGAGGGTCGTTTCaa                                                                                                                                                    <  2:3625131/150‑1 (MQ=255)
       gTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTTTAACGCTGCGTGGGTGCAAGTTGCTTTTAGTTTGTCGTTTGAATATaa                                                                                                                                               >  1:1534665/1‑150 (MQ=255)
       gTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAACGCTGCGTGCGTGCAAGTTTCTTTTCGATGGTCGTTTCAAAATaa                                                                                                                                               >  2:1133654/1‑150 (MQ=255)
        tCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATaaa                                                                                                                                              >  2:363831/1‑150 (MQ=255)
          aCGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACg                                                                                                                                            <  1:697284/150‑1 (MQ=255)
          aCGTCGGCAATCTCGTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACg                                                                                                                                            <  1:3263021/150‑1 (MQ=255)
           cGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGa                                                                                                                                           >  2:2603894/1‑150 (MQ=255)
           cGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCAGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAATATaaaaca                                                                                                                                           >  2:3054672/1‑147 (MQ=255)
            gTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGAc                                                                                                                                          >  1:2174537/1‑150 (MQ=255)
            gTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGAc                                                                                                                                          <  2:3597940/150‑1 (MQ=255)
            gTCGGCAATCTCTTCATACAGCGGATTGCGTTCATGGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGAc                                                                                                                                          <  1:3505349/150‑1 (MQ=255)
             tCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACa                                                                                                                                         <  2:1567388/150‑1 (MQ=255)
              cGGCAATCTCGTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACaa                                                                                                                                        <  1:1014988/150‑1 (MQ=255)
                         tCATACAGCGGATTGCGTTCTTTGGCCAACGATTCCAGAACTTCACGGGGCGGTGTTTCAACGTGCAGAAACGGGCGTTTTTTATCACGCTGCGTGCGTGAAAGTTGCTTTTCGATGGTCGTTTCAAGATACACGACAACGCCACGAGCg                                                                                                                             <  1:3011824/150‑1 (MQ=255)
                         tCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAACCCTGCGGGCGTGCAAGTTGCTTTTAGAAGGTCGTTTAAAGATAAAACACAAAGCCACCAGcc                                                                                                                             >  1:1959078/1‑149 (MQ=255)
                          cATACAGCGGTTTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCgg                                                                                                                            <  2:4153396/150‑1 (MQ=255)
                              ccgCGGATTGCTTTCATGGGCAAACGATACCACAACTTCCCGCGGCGGTTTTTCAACGTGCAGCAACGGGCTTTTTTTATCCCGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTACAAGATAAACGACAACGCCACGAGCGGAAAg                                                                                                                        <  1:3237674/148‑1 (MQ=255)
                              cAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCCAGTTGCTTTTAGATGGTCGTTTCAAGATAAAAGACAACGCCACGAGCGGAAAg                                                                                                                        >  2:793388/1‑150 (MQ=255)
                                gggttttGGGTTAATTGGCCAACGCTTCCAGAACTTCACGGGGGGGTTTTACAACTTGCAGAAAGGGGGGTTTTTTATCACGCTGGGTGCGTGCATTTTGTTTTTCGTTGTTCTTTTAAAGAAAAACGACAACGCCACGAGCGGAAAGAc                                                                                                                      <  1:4350530/145‑1 (MQ=255)
                                gCGGTTTGGGTTCTTTGGCCAACGCTTCCAGAACTTCACGCGGGGGTGTTTCAACGTGCAGCAAGGGGCGTTTTTTATCACGCTGCGTGCGTGCATGTTGCTTTTCGTGGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGAc                                                                                                                      <  1:4115643/150‑1 (MQ=255)
                                 cGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACg                                                                                                                     >  1:762353/1‑150 (MQ=255)
                                 cGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAACGCTGCGTGCGTGCAAGTTGCTTTTAGATGGTCGTTTAAATATAAAAAAAAAACCCACCAGAGGAAAGaag                                                                                                                     >  1:564198/1‑148 (MQ=255)
                                  ggATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACTCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAACGCTGAGTGCGTGCAAGTTGCTTTTAGATGGTCGTTTGAAGATAAAAAACAAAGCCACGATCGGAAAGAAgg                                                                                                                    >  1:2970058/1‑150 (MQ=255)
                                  ggATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTTCTTTTAGATGGTAGTTTGAAGATAAAAGAAAACCCCCCCAGCGGAAACACgg                                                                                                                    >  1:4237656/1‑150 (MQ=255)
                                  ggATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAAAGACAAAGCCACCAGCGGAAAGACgg                                                                                                                    >  1:2971003/1‑150 (MQ=255)
                                  ggATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGCCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACgg                                                                                                                    >  2:1854030/1‑150 (MQ=255)
                                   gATTGCGTTCATTGGCCACCGCTTCCAGAACTGCCCGCGGCGGTGTTTCAACGGGCAGCAACGGGCGTTTTTTATCCCGCGGCGTGCGTGCATGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGt                                                                                                                   <  1:263054/150‑1 (MQ=255)
                                    aTTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGtt                                                                                                                  <  2:1747016/150‑1 (MQ=255)
                                    aTTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGtt                                                                                                                  <  2:2635362/150‑1 (MQ=255)
                                          ttCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGt                                                                                                            >  2:366847/1‑150 (MQ=255)
                                            cATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAAAGACAACGCCACGAGCGGAAAGACGGTTACGCGttt                                                                                                          >  2:4162217/1‑150 (MQ=255)
                                              ttGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCa                                                                                                        >  2:1038464/1‑150 (MQ=255)
                                                ggCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTAGATGGTCGTTTAAATATAAAAGACAACGCCACCAGCGGAAACACGGGTAAGCGTTTAAcc                                                                                                      >  1:3929476/1‑149 (MQ=255)
                                                       gCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAAAGACAACGCCCCGAGCGGAAAGACGGTTACGCGTTTCACGGGAttta                                                                                               >  1:3739437/1‑149 (MQ=255)
                                                       gCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTGTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGGCGTTTCAAGATAAACGACAACGCCACCAGCGGAAAGACGGTTACGCGGTTCACGGGCTTTc                                                                                               >  2:2695239/1‑150 (MQ=255)
                                                        cTTCCAGAACTTCACGCGGCGGTGTTTCACCGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTca                                                                                              <  2:1675357/150‑1 (MQ=255)
                                                        cTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTca                                                                                              <  2:2060524/150‑1 (MQ=255)
                                                          tCCAGAACTTCACGCGGCGGTGTTTCACCGTGCAGCAACGGGCGTTTTTTATCCCGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTcaca                                                                                            <  2:3116035/150‑1 (MQ=255)
                                                             aGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGAAGGTTACGCGTTTCACGGGATTTCACagag                                                                                         >  1:1487154/1‑150 (MQ=255)
                                                               aaCTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGTTGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAgcc                                                                                       <  2:598481/150‑1 (MQ=255)
                                                               aaCTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAgcc                                                                                       <  2:2470153/150‑1 (MQ=255)
                                                                 cTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAgccgc                                                                                     >  2:2387855/1‑150 (MQ=255)
                                                                 cTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAgccgc                                                                                     >  2:892411/1‑150 (MQ=255)
                                                                 cTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAgccgc                                                                                     >  2:3001869/1‑150 (MQ=255)
                                                                    cACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAAAGCCACCAGCGGAAAGACGGTTAAGCGTTTAACGGGATTTCACAGAgccgccgc                                                                                  >  2:1268598/1‑150 (MQ=255)
                                                                          gtggTTTTTCAACGTTCTGCACGGGGCTTTTTTTATCTCGGTGCGTGGGTTAATTTTGCTTTGCGTTGGTCGTTTAAATATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAg                                                                            <  1:4384591/148‑1 (MQ=255)
                                                                           gggTGTTTCACGGTGCAGACAGGGGGTTTTTTTATCCCGGTGCGGGGGTGATTTTTGTTTTTCGTTGTTGTTTAAAATAAACACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGc                                                                           <  1:1133654/149‑1 (MQ=255)
                                                                           cggTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGc                                                                           <  2:2266170/150‑1 (MQ=255)
                                                                           cggTGTTTAAACGTGAAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGAAAGTTGATTTTCGATGGTCTTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGc                                                                           <  2:3484869/150‑1 (MQ=255)
                                                                              tgtTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAg                                                                        <  2:4372255/150‑1 (MQ=255)
                                                                                 ttCAACGTGCAGAAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTAAAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGcac                                                                     <  2:42038/150‑1 (MQ=255)
                                                                                    aaCGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAAt                                                                  <  2:4437093/150‑1 (MQ=255)
                                                                                       gggcagcaAGGGGGTTTTTTTATCCCGGTGGGTGGGTGCATTTTGCTTTTCTTTGGTGTTTAAAATAAAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATAcc                                                               <  1:2688281/148‑1 (MQ=255)
                                                                                         gcagcaACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCt                                                             <  2:861582/150‑1 (MQ=255)
                                                                                          cagcaACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTg                                                            <  2:4064324/150‑1 (MQ=255)
                                                                                            gcaACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGtt                                                          <  1:1844324/150‑1 (MQ=255)
                                                                                             caACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTTCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCCCAAAACCCTGttt                                                         >  2:2484121/1‑150 (MQ=255)
                                                                                             caACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAACCGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGttt                                                         <  2:514262/150‑1 (MQ=255)
                                                                                                gggggGTTTTTTATCCCGCTGGGTGCGTGCATGTTGGTTTTCGTTGGTGGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTtctc                                                      <  1:4299338/149‑1 (MQ=255)
                                                                                                    cGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTc                                                  >  1:2572842/1‑150 (MQ=255)
                                                                                                     gTTTTTTATCCCGCTGCGTGCGTGAACGTTGCTTTTCGATGGTGTTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCa                                                 <  1:3493928/150‑1 (MQ=255)
                                                                                                      ttttttATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCaa                                                <  1:2540524/150‑1 (MQ=255)
                                                                                                        ttttATCACGCTGCGTGCGTGCATGTTGCTTTTCGATGTTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACt                                              <  1:1287887/150‑1 (MQ=255)
                                                                                                         tttATCACGGTGGGTGGGTGAATGTTGCTTTTCGTGGGTCGTTTAAATATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTc                                             <  1:3914648/150‑1 (MQ=255)
                                                                                                           tATCACGCTGCGTGCGTGAAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCAt                                           <  1:1306643/150‑1 (MQ=255)
                                                                                                             tCACGCTGCGTGCGTGCATGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATtg                                         <  2:3725145/150‑1 (MQ=255)
                                                                                                             tCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATtg                                         <  1:1076551/150‑1 (MQ=255)
                                                                                                             tCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATtg                                         >  2:1842290/1‑150 (MQ=255)
                                                                                                             tCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATAACCTGTTTCTCGGTCAACTCATtg                                         >  2:1842160/1‑150 (MQ=255)
                                                                                                              cACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGATCCGCCGCCAGTATCCAGCACAAAACCCTGTTTATCGGTCAACTCATtga                                        >  1:2806963/1‑150 (MQ=255)
                                                                                                               aCGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCCTtgat                                       >  1:3692595/1‑150 (MQ=255)
                                                                                                                   tgcgtgcgTGCATGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGAcc                                   <  2:2345515/150‑1 (MQ=255)
                                                                                                                    gcgtgcgtGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCt                                  >  2:4262091/1‑150 (MQ=255)
                                                                                                                     cgtgggtgCTAGTTGTTTTTCGTTGTTCTTTTCAATATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCtt                                 <  1:781279/150‑1 (MQ=255)
                                                                                                                     cgtgcgtgCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCtt                                 >  2:2924937/1‑150 (MQ=255)
                                                                                                                      gtgcgtgcAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCttt                                <  2:335638/150‑1 (MQ=255)
                                                                                                                         cgtgcAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTttct                             >  2:4462174/1‑150 (MQ=255)
                                                                                                                         cgtgcAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTATCGGTCAACTCATTGATGACCTTttct                             >  1:2920142/1‑150 (MQ=255)
                                                                                                                         cgtgcAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCAGTTTCTCGGTCAACTCATTGATGACCttttat                             >  1:1903162/1‑148 (MQ=255)
                                                                                                                           tgcAAGTTGCTTTTCGAGGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTttcttc                           <  1:2527789/150‑1 (MQ=255)
                                                                                                                           tgaatgTTGATTTTCGATGGTCTTTTAAATATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTttcttc                           <  1:3064864/145‑1 (MQ=255)
                                                                                                                            gcAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGACAGTAGCCAGCTCAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTcg                          <  2:4016236/150‑1 (MQ=255)
                                                                                                                              aaGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCATTTCGTcgcg                        >  1:2370608/1‑150 (MQ=255)
                                                                                                                                      tttCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACGGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGa                <  1:4499257/150‑1 (MQ=255)
                                                                                                                                       ttCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGaa               <  2:3336877/150‑1 (MQ=255)
                                                                                                                                         cGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGc             >  2:2115459/1‑150 (MQ=255)
                                                                                                                                          taTGGTCGTTTCATGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGcc            <  1:2234637/149‑1 (MQ=255)
                                                                                                                                          gATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACATAGCCGCCGCCAGTAGCCAGCACAAAACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGcc            >  1:370452/1‑150 (MQ=255)
                                                                                                                                           agggTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCct           <  2:3797437/148‑1 (MQ=255)
                                                                                                                                           aTGGTCGTTTCAAGATAAACGACATCGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCct           <  2:1917407/150‑1 (MQ=255)
                                                                                                                                            tGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCctt          <  2:3539701/150‑1 (MQ=255)
                                                                                                                                             ggTCTTTTAAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCcttc         <  1:548538/150‑1 (MQ=255)
                                                                                                                                              gTCGTTTCAAGATACACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCcttct        <  2:2779960/150‑1 (MQ=255)
                                                                                                                                              gTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCcttct        <  2:518100/150‑1 (MQ=255)
                                                                                                                                               tCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCcttctt       <  2:1393172/150‑1 (MQ=255)
                                                                                                                                                cGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCttcttc      <  1:971284/150‑1 (MQ=255)
                                                                                                                                                   ttCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGcc   >  1:2188032/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCT  >  CP009273/3511965‑3512262

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑