Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,577,836 C→A 9.8% D283Y (GAC→TAC)  yjiA ← GTP‑binding protein, putative GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,577,8360CA 9.8% 990.2 / 53.8 275D283Y (GAC→TAC) yjiAGTP‑binding protein, putative GTPase
Reads supporting (aligned to +/- strand):  ref base C (118/130);  new base A (18/9);  total (136/139)
Fisher's exact test for biased strand distribution p-value = 6.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

GGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGT  >  CP009273/4577688‑4577985
                                                                                                                                                    |                                                                                                                                                     
 gCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCgg                                                                                                                                                     >  2:3282887/1‑150 (MQ=255)
        gCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTa                                                                                                                                              <  1:4377717/150‑1 (MQ=255)
           gCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAACCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTagag                                                                                                                                           <  1:1751490/150‑1 (MQ=255)
            cTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAgg                                                                                                                                          <  1:1224282/150‑1 (MQ=255)
            cTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAgg                                                                                                                                          >  2:4257873/1‑150 (MQ=255)
            cTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTAGGCGCTGTAGAgg                                                                                                                                          >  1:2332380/1‑150 (MQ=255)
                 ccTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTAGGCGCTGTAGAGGCGCTg                                                                                                                                     >  1:3373975/1‑150 (MQ=255)
                  cTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTgg                                                                                                                                    >  2:681992/1‑150 (MQ=255)
                   ttATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGa                                                                                                                                   >  2:3800397/1‑150 (MQ=255)
                     aTCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACg                                                                                                                                 >  2:2739529/1‑150 (MQ=255)
                        cGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGccc                                                                                                                              <  2:4532210/150‑1 (MQ=255)
                        cGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGccc                                                                                                                              <  2:3978264/150‑1 (MQ=255)
                         ggCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCAACGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCt                                                                                                                             >  2:539845/1‑150 (MQ=255)
                          gCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTg                                                                                                                            <  2:3085990/150‑1 (MQ=255)
                          gCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTg                                                                                                                            >  2:2569860/1‑150 (MQ=255)
                           cccATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTgg                                                                                                                           <  2:568538/150‑1 (MQ=255)
                            ccATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGa                                                                                                                          >  1:4023497/1‑150 (MQ=255)
                              aTAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAc                                                                                                                        <  1:1428283/150‑1 (MQ=255)
                              aTAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAc                                                                                                                        <  2:1935707/150‑1 (MQ=255)
                              aTAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTAGGCGCTGTAGAGGCGCTGGACGCCCTGGAAc                                                                                                                        >  2:2697344/1‑150 (MQ=255)
                                  ccccTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAcagca                                                                                                                    <  1:2735436/150‑1 (MQ=255)
                                  ccccTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTAGGCGCTGTAGAGGCGCTGGACGCCCTGGACcagca                                                                                                                    >  2:904710/1‑150 (MQ=255)
                                     cTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATCAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGc                                                                                                                 <  2:2070686/150‑1 (MQ=255)
                                      ttACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCg                                                                                                                >  1:117769/1‑150 (MQ=255)
                                       tACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTAGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCgg                                                                                                               >  1:1349609/1‑150 (MQ=255)
                                        aCTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGt                                                                                                              <  2:3301440/150‑1 (MQ=255)
                                        aCTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGt                                                                                                              <  2:3301503/150‑1 (MQ=255)
                                         cgtCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGtt                                                                                                             <  1:3744606/148‑1 (MQ=255)
                                         cTTCCTCAACCCCGCAACCGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGtt                                                                                                             <  1:2702547/150‑1 (MQ=255)
                                         cTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGtt                                                                                                             >  1:4222138/1‑150 (MQ=255)
                                         cTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGtt                                                                                                             <  1:1666646/150‑1 (MQ=255)
                                          ttCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTa                                                                                                            >  1:862301/1‑150 (MQ=255)
                                          ttCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTa                                                                                                            <  1:2497961/150‑1 (MQ=255)
                                           tCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAg                                                                                                           >  1:2490312/1‑150 (MQ=255)
                                           tCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAg                                                                                                           <  2:2642806/150‑1 (MQ=255)
                                            ccTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAgg                                                                                                          >  2:479235/1‑150 (MQ=255)
                                             cTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGt                                                                                                         <  2:1677962/150‑1 (MQ=255)
                                              tCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGtt                                                                                                        <  2:2552694/150‑1 (MQ=255)
                                                  ccccGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTAGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGcc                                                                                                    >  1:2030597/1‑150 (MQ=255)
                                                   cccGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCg                                                                                                   <  1:3812832/150‑1 (MQ=255)
                                                   cccGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCg                                                                                                   >  1:1703071/1‑150 (MQ=255)
                                                    ccGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGt                                                                                                  >  2:120002/1‑150 (MQ=255)
                                                     cGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTc                                                                                                 <  2:1053406/150‑1 (MQ=255)
                                                        aaaCGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAAt                                                                                              <  1:567471/150‑1 (MQ=255)
                                                           cGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCca                                                                                           >  1:2229216/1‑150 (MQ=255)
                                                            gCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCcac                                                                                          <  1:309472/150‑1 (MQ=255)
                                                            gCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCcac                                                                                          >  1:2530605/1‑150 (MQ=255)
                                                            gCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCcac                                                                                          >  1:2957606/1‑150 (MQ=255)
                                                             cAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCcaca                                                                                         <  1:494981/150‑1 (MQ=255)
                                                              aGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAg                                                                                        >  2:1742518/1‑150 (MQ=255)
                                                               gCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGc                                                                                       <  1:3013073/150‑1 (MQ=255)
                                                                  cGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATc                                                                                    <  1:3610710/150‑1 (MQ=255)
                                                                   gAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATcc                                                                                   >  2:3474713/1‑150 (MQ=255)
                                                                   gAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATcc                                                                                   >  2:2023535/1‑150 (MQ=255)
                                                                                  gCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCa                                                                    <  1:999592/150‑1 (MQ=255)
                                                                                    aGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGt                                                                  <  1:396062/150‑1 (MQ=255)
                                                                                     gCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGtt                                                                 <  1:1926118/150‑1 (MQ=255)
                                                                                     gCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGtt                                                                 >  2:1023722/1‑150 (MQ=255)
                                                                                       tGGATCCCGATAAGCACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTa                                                               >  2:50421/1‑150 (MQ=255)
                                                                                       tGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTa                                                               >  2:2501128/1‑150 (MQ=255)
                                                                                       tGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTa                                                               >  2:2501284/1‑150 (MQ=255)
                                                                                         gATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGGCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATc                                                             >  1:2939842/1‑150 (MQ=255)
                                                                                          aTCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGAAGCAGTTTATca                                                            >  1:4080687/1‑149 (MQ=255)
                                                                                             ccGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATccgcc                                                         <  2:2673148/150‑1 (MQ=255)
                                                                                              cGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATccgccg                                                        <  2:3477047/150‑1 (MQ=255)
                                                                                              cGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATccgccg                                                        <  2:1862597/150‑1 (MQ=255)
                                                                                                aTAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGAc                                                      >  1:2908557/1‑150 (MQ=255)
                                                                                                aTAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGAc                                                      >  1:2908518/1‑150 (MQ=255)
                                                                                                aTAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGAc                                                      >  2:4512531/1‑150 (MQ=255)
                                                                                                aTAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTAGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGAc                                                      >  2:4100928/1‑150 (MQ=255)
                                                                                                   aaCACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTcc                                                   >  1:68128/1‑150 (MQ=255)
                                                                                                    acacCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCca                                                  >  1:2348051/1‑150 (MQ=255)
                                                                                                     cacCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCcag                                                 >  1:1041887/1‑150 (MQ=255)
                                                                                                          tCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCcagcagca                                            <  2:2834262/150‑1 (MQ=255)
                                                                                                             tGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTAGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGt                                         >  1:4541006/1‑150 (MQ=255)
                                                                                                              gCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTAGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGtt                                        <  2:2332380/150‑1 (MQ=255)
                                                                                                                tATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGtttt                                      <  2:1062994/150‑1 (MQ=255)
                                                                                                                 aTGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTc                                     >  2:1916126/1‑150 (MQ=255)
                                                                                                                   gCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCa                                   >  1:3330070/1‑150 (MQ=255)
                                                                                                                   gCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCa                                   >  2:4032582/1‑150 (MQ=255)
                                                                                                                    cGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCAt                                  <  2:2298878/150‑1 (MQ=255)
                                                                                                                    cGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCAt                                  <  2:3735674/150‑1 (MQ=255)
                                                                                                                    cGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCAt                                  >  1:2345753/1‑150 (MQ=255)
                                                                                                                     ggTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATc                                 <  2:6360/150‑1 (MQ=255)
                                                                                                                      gTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCa                                >  2:1496499/1‑150 (MQ=255)
                                                                                                                           tCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCgg                           <  2:4524408/150‑1 (MQ=255)
                                                                                                                            cATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCggg                          <  2:215798/150‑1 (MQ=255)
                                                                                                                             aTCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGa                         <  2:1831054/150‑1 (MQ=255)
                                                                                                                                 ccccACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACt                     >  2:2174939/1‑150 (MQ=255)
                                                                                                                                 ccccACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACt                     <  1:3556360/150‑1 (MQ=255)
                                                                                                                                   ccACGGCCTGTCCCAGTAGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTc                   <  1:1436308/150‑1 (MQ=255)
                                                                                                                                    cACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCg                  >  2:1628351/1‑150 (MQ=255)
                                                                                                                                      cGGCCTGTCCCAGTCGGCGCTGTAGAGGCGGTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACACCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCt                <  2:2913025/150‑1 (MQ=255)
                                                                                                                                       ggCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTg               >  1:2571199/1‑150 (MQ=255)
                                                                                                                                           tGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGatat           >  1:2638560/1‑150 (MQ=255)
                                                                                                                                               ccAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTAc       <  1:1794329/150‑1 (MQ=255)
                                                                                                                                                cAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTAcc      <  2:2965611/150‑1 (MQ=255)
                                                                                                                                                 aGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCg     <  2:28871/150‑1 (MQ=255)
                                                                                                                                                 aGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCg     <  2:1333641/150‑1 (MQ=255)
                                                                                                                                                   tCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCggg   <  1:286183/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGT  >  CP009273/4577688‑4577985

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