Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 51,352 A→T 17.6% G85G (GGT→GGA apaG ← protein associated with Co2+ and Mg2+ efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927351,3520AT17.6% 921.5 / 181.2 347G85G (GGT→GGAapaGprotein associated with Co2+ and Mg2+ efflux
Reads supporting (aligned to +/- strand):  ref base A (147/138);  new base T (38/23);  total (185/162)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 348 total aligned reads displayed.

CAATAAGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCAGATTGCG  >  CP009273/51204‑51501
                                                                                                                                                    |                                                                                                                                                     
cAATAAGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCAc                                                                                                                                                      >  1:32023/1‑150 (MQ=255)
cAATAAGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCAc                                                                                                                                                      >  1:312656/1‑150 (MQ=255)
 aaTAAGGTATGTCGCCATATTCTTTTAATGATTGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCAcc                                                                                                                                                     <  1:1510709/150‑1 (MQ=255)
 aaTAAGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTcc                                                                                                                                                     <  2:3985539/150‑1 (MQ=255)
    aaGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCt                                                                                                                                                  >  1:3873010/1‑150 (MQ=255)
     aGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGGCCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTg                                                                                                                                                 <  2:59130/150‑1 (MQ=255)
       gTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGgt                                                                                                                                               >  1:4229070/1‑150 (MQ=255)
        tATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGgtg                                                                                                                                              <  2:2399258/150‑1 (MQ=255)
             cGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTg                                                                                                                                         >  2:4094787/1‑150 (MQ=255)
              gCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTgg                                                                                                                                        >  1:3808471/1‑150 (MQ=255)
              gCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTgg                                                                                                                                        >  1:1654310/1‑150 (MQ=255)
               ccATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGt                                                                                                                                       <  2:2057754/150‑1 (MQ=255)
                cATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTa                                                                                                                                      >  2:479469/1‑150 (MQ=255)
                 atatTCTTTTAATGTATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGCACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTAc                                                                                                                                     <  1:4439838/150‑1 (MQ=255)
                 atatTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTAc                                                                                                                                     >  2:2329376/1‑150 (MQ=255)
                  tatTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTAct                                                                                                                                    >  2:2759757/1‑150 (MQ=255)
                   atTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTActc                                                                                                                                   >  2:4420073/1‑150 (MQ=255)
                         ttAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGcc                                                                                                                             >  1:450549/1‑150 (MQ=255)
                          tAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCa                                                                                                                            <  2:171721/150‑1 (MQ=255)
                           aatgaatgAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAg                                                                                                                           >  1:471462/1‑150 (MQ=255)
                           aatgaatgAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAg                                                                                                                           >  1:2021007/1‑150 (MQ=255)
                                atgaGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGcgcg                                                                                                                      >  1:2915934/1‑150 (MQ=255)
                                 tgaGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGa                                                                                                                     >  2:355159/1‑150 (MQ=255)
                                 tgaGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGa                                                                                                                     <  1:2325009/150‑1 (MQ=255)
                                  gaGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGAt                                                                                                                    >  1:4737053/1‑150 (MQ=255)
                                   aGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATa                                                                                                                   <  2:216917/150‑1 (MQ=255)
                                     tgtgGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAg                                                                                                                 >  2:4721164/1‑150 (MQ=255)
                                      gtgGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGt                                                                                                                >  1:4962904/1‑150 (MQ=255)
                                       tgGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTg                                                                                                               <  1:4178935/150‑1 (MQ=255)
                                        gggAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTgg                                                                                                              >  1:1030806/1‑150 (MQ=255)
                                        gggAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTgg                                                                                                              <  2:1300503/150‑1 (MQ=255)
                                           aaCGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTg                                                                                                           <  2:28049/150‑1 (MQ=255)
                                              ggCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGAc                                                                                                        >  2:2303674/1‑150 (MQ=255)
                                              ggCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCAAGGCGCGATAAGTGGCTGGAc                                                                                                        <  2:3412150/150‑1 (MQ=255)
                                               gCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACg                                                                                                       >  2:164774/1‑150 (MQ=255)
                                                 gAGTCGGAAGACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGnc                                                                                                     <  1:1408472/150‑1 (MQ=255)
                                                    tCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAAc                                                                                                  >  2:3067756/1‑150 (MQ=255)
                                                     cGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTGCGTAGTGACCCAGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGAACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAAcc                                                                                                 >  2:2075922/1‑150 (MQ=255)
                                                     cGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATATTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAAcc                                                                                                 >  1:530146/1‑150 (MQ=255)
                                                       gAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCAc                                                                                               >  1:3538913/1‑150 (MQ=255)
                                                         aTACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTc                                                                                             >  2:3880523/1‑150 (MQ=255)
                                                             gggAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTc                                                                                         <  2:1913834/150‑1 (MQ=255)
                                                             gggAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTc                                                                                         <  1:2666581/150‑1 (MQ=255)
                                                              ggAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCg                                                                                        <  1:4563879/150‑1 (MQ=255)
                                                                aaTGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGcc                                                                                      >  2:4509219/1‑150 (MQ=255)
                                                                aaTGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGcc                                                                                      <  1:3963110/150‑1 (MQ=255)
                                                                aaTGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGcc                                                                                      >  1:3000999/1‑150 (MQ=255)
                                                                aaTGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGcc                                                                                      <  1:3282375/150‑1 (MQ=255)
                                                                  tGTAGATGCTGAAAGGGACGCCATTTTCATCGATCATTTGGTTGTGTCCCTGCATGGTGCACAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCt                                                                                    <  1:3978833/150‑1 (MQ=255)
                                                                  ggTCGAGGCTGACCGGGACGCCAGGTTCATCGATCAGGTCGTAGTGACCCTGCATGGGGCCCAGCGGGGTTTCACTGCTTGCACCGCGGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCCCTCCGTCGCCCt                                                                                    <  2:3523470/149‑1 (MQ=255)
                                                                   gTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTg                                                                                   <  1:1270989/150‑1 (MQ=255)
                                                                         gCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTc                                                                             <  1:3440735/150‑1 (MQ=255)
                                                                         gCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTc                                                                             <  1:2862876/150‑1 (MQ=255)
                                                                             aaaGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACAACTCCTTCGCCCTGGACTTCGGtt                                                                         <  2:746089/150‑1 (MQ=255)
                                                                               aGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTc                                                                       >  1:834007/1‑150 (MQ=255)
                                                                               aGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTc                                                                       <  2:4306819/150‑1 (MQ=255)
                                                                                  gACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACg                                                                    <  2:2562007/150‑1 (MQ=255)
                                                                                     gCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgcc                                                                 <  1:4822480/150‑1 (MQ=255)
                                                                                      ccATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgcca                                                                >  2:2703255/1‑150 (MQ=255)
                                                                                          tttCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgc                                                            >  2:2614045/1‑150 (MQ=255)
                                                                                          tttCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgc                                                            >  1:1101032/1‑150 (MQ=255)
                                                                                            tCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgcca                                                          >  1:2882187/1‑150 (MQ=255)
                                                                                            tCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgcca                                                          <  1:69620/150‑1 (MQ=255)
                                                                                             cATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgccat                                                         <  1:3426167/150‑1 (MQ=255)
                                                                                              aTCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgccatt                                                        >  2:4248306/1‑150 (MQ=255)
                                                                                                cGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTgg                                                      >  2:804087/1‑150 (MQ=255)
                                                                                                cGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTgg                                                      >  2:3168429/1‑150 (MQ=255)
                                                                                                 gATCATTTCGGAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGt                                                     <  2:3332868/150‑1 (MQ=255)
                                                                                                          gTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCa                                            >  2:1705091/1‑150 (MQ=255)
                                                                                                               gACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAAc                                       <  2:1438896/150‑1 (MQ=255)
                                                                                                               gACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTTCTGGTTCTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAAc                                       <  1:2429230/150‑1 (MQ=255)
                                                                                                                  ccTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGcc                                    <  1:3636020/150‑1 (MQ=255)
                                                                                                                   gtGCGTGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCGTTGGTGATCAGCCAGTAACGccc                                   >  1:2487894/2‑150 (MQ=255)
                                                                                                                      cATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaa                                >  2:1103201/1‑150 (MQ=255)
                                                                                                                      cATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaa                                <  1:2078537/150‑1 (MQ=255)
                                                                                                                       aTGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaac                               >  1:2939244/1‑150 (MQ=255)
                                                                                                                       aTGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaac                               <  1:2373893/150‑1 (MQ=255)
                                                                                                                         ggTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaacaa                             >  2:4624640/1‑150 (MQ=255)
                                                                                                                          gTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCaacaac                            >  2:3745667/1‑150 (MQ=255)
                                                                                                                             cccAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGc                         >  2:2964687/1‑150 (MQ=255)
                                                                                                                             cccAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGc                         >  2:510292/1‑150 (MQ=255)
                                                                                                                             cccAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGc                         <  2:1153016/150‑1 (MQ=255)
                                                                                                                                  cGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTgg                    <  1:4868551/150‑1 (MQ=255)
                                                                                                                                    gggTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGgcg                  >  2:4941330/1‑150 (MQ=255)
                                                                                                                                    gggTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGgcg                  <  1:4598314/150‑1 (MQ=255)
                                                                                                                                     ggTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGgcgc                 >  2:5125897/1‑150 (MQ=255)
                                                                                                                                     ggTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGgcgc                 >  2:1627092/1‑150 (MQ=255)
                                                                                                                                        ttCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCg              <  2:2915676/150‑1 (MQ=255)
                                                                                                                                         tCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGc             >  2:3735544/1‑150 (MQ=255)
                                                                                                                                           aaTGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGccc           <  1:159612/150‑1 (MQ=255)
                                                                                                                                           aaTGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGccc           <  1:2249168/150‑1 (MQ=255)
                                                                                                                                            aTGATTGCTCCGCTGGTGTACTGGTACTCGTCGCCAGGCGCGATAAGTGGCTGGTCGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCATCCAGTAACGCCCCAACAACTGCACTGGCGCTCGcccc          <  1:5039701/150‑1 (MQ=255)
                                                                                                                                            aTGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCTCTGGCGCTCGcccc          <  2:861507/150‑1 (MQ=255)
                                                                                                                                            aTGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGcccc          <  1:136345/150‑1 (MQ=255)
                                                                                                                                            aTGAGTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGcccc          <  2:4885732/150‑1 (MQ=255)
                                                                                                                                             tGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCa         >  1:3578369/1‑150 (MQ=255)
                                                                                                                                             tGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCa         >  2:2417270/1‑150 (MQ=255)
                                                                                                                                              gATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCAg        >  2:1373544/1‑150 (MQ=255)
                                                                                                                                                ttGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCAGAt      >  1:630592/1‑150 (MQ=255)
                                                                                                                                                 tGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCAGAtt     <  1:1753587/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAATAAGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCAGATTGCG  >  CP009273/51204‑51501

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