Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,453,434 C→A 6.2% Q87K (CAA→AAA)  gspF → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,453,4340CA 6.2% 911.3 / 2.9 261Q87K (CAA→AAA) gspFgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (159/74);  new base A (2/26);  total (161/100)
Fisher's exact test for biased strand distribution p-value = 2.76e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

TCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCTGGTAAAAGCGGGCGAAAAGAGC  >  CP009273/3453285‑3453583
                                                                                                                                                       |                                                                                                                                                     
tCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTc                                                                                                                                                       >  1:1313559/1‑150 (MQ=255)
 ccTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTca                                                                                                                                                      >  2:3971809/1‑150 (MQ=255)
   tgctgGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTcaac                                                                                                                                                    <  1:165722/150‑1 (MQ=255)
       ggATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAg                                                                                                                                                >  1:793393/1‑150 (MQ=255)
       ggATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAg                                                                                                                                                >  2:5001812/1‑150 (MQ=255)
        gATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGc                                                                                                                                               <  2:515274/150‑1 (MQ=255)
        gATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGc                                                                                                                                               >  1:4969207/1‑150 (MQ=255)
         atatTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCa                                                                                                                                              >  2:4292273/1‑150 (MQ=255)
               gCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGtaata                                                                                                                                        <  1:1675074/150‑1 (MQ=255)
                cccccAAAACAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGtaataa                                                                                                                                       <  2:5387791/150‑1 (MQ=255)
                 ccccAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGCACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATaaa                                                                                                                                      <  1:4676917/150‑1 (MQ=255)
                   ccAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACg                                                                                                                                    <  2:2727986/150‑1 (MQ=255)
                     aaaaaaGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGAc                                                                                                                                  >  1:3574845/1‑150 (MQ=255)
                     aaaaaaGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGAc                                                                                                                                  <  1:2130080/150‑1 (MQ=255)
                     aaaaaaGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGAc                                                                                                                                  >  2:3299353/1‑150 (MQ=255)
                     aaaaaaGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGAc                                                                                                                                  <  2:4341235/150‑1 (MQ=255)
                      aaaaaGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACt                                                                                                                                 >  1:3998745/1‑150 (MQ=255)
                      aaaaaGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACt                                                                                                                                 >  1:380806/1‑150 (MQ=255)
                       aaaaGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTg                                                                                                                                <  2:2856572/150‑1 (MQ=255)
                        aaaGTTCGGGAGTAAAAACACGTCGCCTGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTgg                                                                                                                               >  1:1837446/1‑150 (MQ=255)
                             tCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGAc                                                                                                                          >  1:5239525/1‑150 (MQ=255)
                              cGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACg                                                                                                                         >  1:4327418/1‑150 (MQ=255)
                                ttataaaaaaCACGTCGCCCTAGTACCAGCATTTTTTAACTGCCGTTTATCACCCGCATTTTGTAAAATTTAAGCGCAGCGTATTTACCCCTTTAATAGAGCTTTGCCTTAATGGTTAAAAAAAGAATTAATAAACGACTGGGTGACgtg                                                                                                                       <  2:775277/145‑1 (MQ=255)
                                 gAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACgtgt                                                                                                                      >  2:106230/1‑150 (MQ=255)
                                  attAAAAACCCGGCGCCCGTGTTTCAGCCTTTGTTAACGGCCGCTTTTCCCCGGGCTGTTGGCAACCTTACGCGCCGCGGCATTACCCGTGTATGTGAGCGTTGCCTTATTGGGTCAAAAACGCTTTAATAAACGACTGGGTGACGTGtt                                                                                                                     <  2:2483770/148‑1 (MQ=255)
                                    tAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTaa                                                                                                                   >  1:3173042/1‑150 (MQ=255)
                                    tAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTaa                                                                                                                   >  1:5427403/1‑150 (MQ=255)
                                     aaaaaCACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTaaa                                                                                                                  >  1:2629929/1‑150 (MQ=255)
                                      aaaaCACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAAt                                                                                                                 >  1:3556679/1‑150 (MQ=255)
                                             gTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTAc                                                                                                          >  1:4191011/1‑150 (MQ=255)
                                                gCCCGAGCAACAGCCATTGTGAACGGCCGCTTTCCACCCGGCTGTTGGAAACCTTACGCGCAGCGGAATTACCCGTGGAAGAGAGCCTTGCCTTAATCGGAAAACAAAGCAGTAATACACGACTGGGTGACGTGTTAAATCAGGTACGCa                                                                                                       <  2:1964939/150‑1 (MQ=255)
                                                    gAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAgcgc                                                                                                   >  1:975535/1‑150 (MQ=255)
                                                         tCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATcc                                                                                              >  2:1606268/1‑150 (MQ=255)
                                                         tCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATcc                                                                                              >  1:4505033/1‑150 (MQ=255)
                                                         tCAGCAATAGTGAACTGACGCTTTTCACCCGGATGTGGGCAACATTAAGCGCAGGGGATTTACCCCTGTAAGATAGCCTTGCGTTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATcc                                                                                              <  2:2698597/150‑1 (MQ=255)
                                                           aGCCATTGTGAACTGACGCTTTTCACCGGGCTTTTGGAAACCTTAAGCGCTGCGGAATTCCCCCTGTAAGAGAGCCTTGCCGTAATCGGTCAAAAAAGCATTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCtt                                                                                            <  2:1798503/150‑1 (MQ=255)
                                                           aGCAATTGTGAACTTACGATTTTCACCCGGCAGTTGGCAACCTTAAGCGCCGCGGAATTACCCCTGTAAGAGAGCCTTGCCTTAATCGTTCAACAAAGCATTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCtt                                                                                            <  2:4563239/150‑1 (MQ=255)
                                                             ccATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGa                                                                                          >  1:646727/1‑150 (MQ=255)
                                                              ctttgTGACCTGACGCTTTTCACCCGGCTGTTGGCACCCTTACGCGCAGCGGCATTCCCCCTGGAAGAGAGCCTTGCCGTAATGGGTCAACAAAGCAGTAATACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGaa                                                                                         <  2:367413/146‑1 (MQ=255)
                                                              cATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGaa                                                                                         >  1:5155529/1‑150 (MQ=255)
                                                               aTAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAg                                                                                        >  2:3023656/1‑150 (MQ=255)
                                                               aTAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAg                                                                                        >  1:2645142/1‑150 (MQ=255)
                                                                  gTGAACGGCCGCTTTTCCCCGGGCAGTGGGCAACCTTACGCGCCGCGGCATTCCCCCTGTACGAGCGCCTTGCGTTACTGGGTCAACAACGCTGAAATACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGc                                                                                     <  2:3410658/150‑1 (MQ=255)
                                                                    gAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGcct                                                                                   >  2:4453889/1‑148 (MQ=255)
                                                                    gAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCAt                                                                                   >  2:4151307/1‑150 (MQ=255)
                                                                     aaCTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATc                                                                                  >  1:2703865/1‑150 (MQ=255)
                                                                     aaCTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATc                                                                                  >  1:2881999/1‑150 (MQ=255)
                                                                         gccGCTTTACACCCGCCTGTTGGCAACCTAACGCGCCGCGGCATTCCCCCTGTAATAGCGCCTTGCCGAAATGGTACACAAAAGCAGTAATACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCt                                                                              <  2:3165360/148‑1 (MQ=255)
                                                                         gACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCt                                                                              >  2:4373470/1‑150 (MQ=255)
                                                                          aCGCTTTTCCCCCGGCAGTTGGCAACCTTAAGCGCCGCGGAATCACCCCGGGAAGATAGCCTTGCCTTAAGCGTTCAAAAAAGCATTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCtt                                                                             <  2:3450603/150‑1 (MQ=255)
                                                                             gttttCCCCCGGCTGTTGCCAACCTTACCCGCCGCGGCTTTCCCCCGGTAAGAGAGCCTTGCCTTAATCGGACAACAAACAATAAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTcc                                                                          <  2:4244616/149‑1 (MQ=255)
                                                                             cTTTTCACCCGGCTGTTGGCAACCTTACGCGCCGCGGCTTTACCCCTGGAAGAGAGCCTGGCCGAAATCGGAAAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTcc                                                                          <  2:2189425/150‑1 (MQ=255)
                                                                             cTTTACCCCCGGCTGTTGGCAACCTTAAGCGCCGCGGAATTCCCCCTGTATGAGCGCCTGGCCTTAATCGGAAAAAAACGCATTAAAACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTcc                                                                          <  2:4644154/150‑1 (MQ=255)
                                                                             attttCCCCCGGGGTTTGGCACCCTTCCGCGCGGCGGCATCCCCCCGGGACGAGCGCCTGGCCGTAAGCGGACAACAACGCAGTCATACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTcc                                                                          <  1:3103334/149‑1 (MQ=255)
                                                                                ttCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGAt                                                                       >  2:2098283/1‑150 (MQ=255)
                                                                                  ccccgggCTGTTGTCAACCTTACGCGCCGCGGAATACCCCCTGTAATTGAGCTTTCCCGTAATCGGACAAAAACGAATAAATACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGc                                                                     <  1:3365533/145‑1 (MQ=255)
                                                                                   aCCGGGCTGTTGCCAACCTTACGCGCCGCGGAATACCCCCTGCATGACAGCTTTGCCTTAATGGGAAAAAAACGCTTTAAAACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCa                                                                    <  1:3226213/150‑1 (MQ=255)
                                                                                   aCCCGGCTTTGGGCAACATTACGCGCCGCGGATTAACCCCTGTATGAGAGCCTTGCCTAAATCGGAAAAAAAAGCAGTAAAAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCa                                                                    <  1:459815/150‑1 (MQ=255)
                                                                                    cccGGCTGTTGGCAACCTTACGCGCAGCGGCATTACCCCTGTAAGAGAGCAGTGCCGTAATGGTTAAACAAAGCATTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCAt                                                                   <  2:1335112/150‑1 (MQ=255)
                                                                                    cccGGCAGTTGGCAACCTTAAGCGCAGCGTCATTCCCCCTGGAAGAGAGCGTGGCCTTAATCGTACAACAAAGCATAAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCAt                                                                   <  2:4190728/150‑1 (MQ=255)
                                                                                        cgaGTTGGCACCATTCCGCGCCGCGGCATACCCCCTGCACGAGCGCCTGGCCGTAATGGGACAACACAGCTGTAATACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACa                                                               <  1:387746/148‑1 (MQ=255)
                                                                                           gTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGTTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCTATGCATTACAGCa                                                            >  2:2316127/1‑150 (MQ=255)
                                                                                           gTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCCTCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCa                                                            >  2:2235961/1‑150 (MQ=255)
                                                                                            ttGCCAACTTTACGCGCCGCGGATTACCCCCTGTATGTGAGCCTTGCGTTATTCGGAAAACAAACCATTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCAt                                                           <  1:4995125/150‑1 (MQ=255)
                                                                                            gtGGAAACATTAAGCGCAGCGGCTTTCCCCCTGGAATAGAGCCTTGCCGTATTGGGACAAAAAAGCATAAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCAt                                                           <  2:4022059/149‑1 (MQ=255)
                                                                                               gaaaCCTTCCGCGCCGCGGCATACCCCGGGTAAGAGCGCCTGGCCTTAAGCGGACAACAAAGCAGTAAAACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCAtttt                                                        <  1:4053431/148‑1 (MQ=255)
                                                                                               gCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCAtttt                                                        >  2:2324741/1‑150 (MQ=255)
                                                                                               gCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCAtttt                                                        >  2:4394620/1‑150 (MQ=255)
                                                                                               ccAACCTTAAGCGCCGCGGCATAACCCGTGGAAGAGCGCCTGGCCTTAATCGGTCAACAAAGCATAAATACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCAtttt                                                        <  1:1289359/149‑1 (MQ=255)
                                                                                                cAACCTTAAGCGCAGCGGCATTAACCCTGGAATAGAGCCTTGACGTAATCTGTAAACAAAGAAGTAATAAACGACTGGGTGACGTGTTAAATAAGGTACGCAGCGCCATCCTTTAAGGGCATCCCATTTCCGATGCATTAAAGCATTTTc                                                       >  2:1934636/1‑150 (MQ=255)
                                                                                                 aaCATTACGCGCCGCGGATTCCCCCCTGGATGAGAGCCTTCCCTTAATCGTAAAAAAAACCATAAAAACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTcc                                                      <  1:2226471/150‑1 (MQ=255)
                                                                                                    cTTATGCGCCGCGGCATTACCCCTGGAATAGAGCCTGGCGGTAATGGTTCAACAAAGCAGTAATACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCAc                                                   <  2:4531521/150‑1 (MQ=255)
                                                                                                    cTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCAc                                                   >  1:187956/1‑150 (MQ=255)
                                                                                                      tAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGc                                                 >  2:5241046/1‑150 (MQ=255)
                                                                                                          cgccggggCATTACCCGTGGATGTGCGCCTGGCCTTACTGGGACAACAAAGCTTTAATACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCtttt                                             <  2:608169/144‑1 (MQ=255)
                                                                                                             aGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCCTCCCCTGTCCGATGCATTACAGCATTTTCCCACGCTTTTCGa                                          >  1:1112068/1‑150 (MQ=255)
                                                                                                             aGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGa                                          >  1:5329459/1‑150 (MQ=255)
                                                                                                             aGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGa                                          >  1:1111945/1‑150 (MQ=255)
                                                                                                               cGGAATTACCCCTGTAAGATAGCCTTGCCTTAATCGTTCAAAAAAGCATTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGAtt                                        <  2:1730339/150‑1 (MQ=255)
                                                                                                                ggCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTc                                       >  2:3712354/1‑150 (MQ=255)
                                                                                                                ggCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTc                                       >  1:2038655/1‑150 (MQ=255)
                                                                                                                 gCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCg                                      >  2:1456415/1‑150 (MQ=255)
                                                                                                                  cATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGc                                     >  2:2611164/1‑150 (MQ=255)
                                                                                                                  cATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGc                                     >  1:4143195/1‑150 (MQ=255)
                                                                                                                        ccccTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTAt                                >  1:2077246/1‑150 (MQ=255)
                                                                                                                         cccTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATc                               >  2:2643374/1‑150 (MQ=255)
                                                                                                                          ccgggAAGAGAGCGTGGCCGTAATCGGACAACAACGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCg                              <  2:3075363/147‑1 (MQ=255)
                                                                                                                            tGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTCTCGTa                            >  2:86014/1‑150 (MQ=255)
                                                                                                                            tGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTa                            >  2:3172182/1‑150 (MQ=255)
                                                                                                                              gtaaGAGAGCCTTGCCTTAATCGTTAAACAAAGCAGTAATACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTAc                           <  2:3605548/148‑1 (MQ=255)
                                                                                                                                aaGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTAccc                         >  1:553548/1‑150 (MQ=255)
                                                                                                                                 agagagCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCt                        >  2:190753/1‑150 (MQ=255)
                                                                                                                                  gagagCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCTg                       >  2:481373/1‑150 (MQ=255)
                                                                                                                                   agagCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCTgg                      >  1:4546138/1‑150 (MQ=255)
                                                                                                                                      gCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCTGGTaa                   >  2:3636758/1‑150 (MQ=255)
                                                                                                                                      gCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCTGGTaa                   >  2:3101882/1‑150 (MQ=255)
                                                                                                                                      gCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCTGGTaa                   >  2:5253766/1‑150 (MQ=255)
                                                                                                                                                aaTCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCTGGTAAAAGCGGGCGa         >  2:4715757/1‑150 (MQ=255)
                                                                                                                                                   cGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCTGGTAAAAGCGGGCGaaaa      >  2:539335/1‑150 (MQ=255)
                                                                                                                                                    gcacaacaaACAATTAAAACACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCTGGTAAAAGCGGGCGAAAag     <  1:2125828/147‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCAGGTACGCAGCGCCATCCTTGAAGGGCATCCCCTTTCCGATGCATTACAGCATTTTCCCACGCTTTTCGATTCGCTCTATCGTACCCTGGTAAAAGCGGGCGAAAAGAGC  >  CP009273/3453285‑3453583

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