Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,939,313 G→T 94.6% intergenic (‑112/+127) tcdA ← / ← mltA tRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA/membrane‑bound lytic murein transglycosylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,939,3130GT94.6% 695.8 / 29.2 242intergenic (‑112/+127)tcdA/mltAtRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA/membrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base G (3/10);  new base T (110/119);  total (113/129)
Fisher's exact test for biased strand distribution p-value = 9.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 242 total aligned reads displayed.

GCTGACGCCATGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGT  >  CP009273/2939168‑2939456
                                                                                                                                                 |                                                                                                                                               
     cGCCATGCATCACTAATTACCACAGACAGAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCAACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTa                                                                                                                                          <  2:2647756/148‑1 (MQ=255)
       ccATGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAgg                                                                                                                                        <  1:3508918/148‑1 (MQ=255)
           gCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCgg                                                                                                                                    >  1:3024032/1‑148 (MQ=255)
           gCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCgg                                                                                                                                    <  1:2202314/148‑1 (MQ=255)
            cATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTACTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGa                                                                                                                                   <  1:530763/148‑1 (MQ=255)
              tCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATa                                                                                                                                 >  1:2733882/1‑148 (MQ=255)
              tCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATa                                                                                                                                 >  1:1419806/1‑148 (MQ=255)
              tCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATa                                                                                                                                 >  2:2333051/1‑148 (MQ=255)
              tCACTAAGTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATa                                                                                                                                 <  1:420051/148‑1 (MQ=255)
                 cTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAgg                                                                                                                              >  2:890504/1‑148 (MQ=255)
                    aTTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGt                                                                                                                           <  1:371146/148‑1 (MQ=255)
                    aTTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGt                                                                                                                           >  1:900222/1‑148 (MQ=255)
                      tACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTc                                                                                                                         >  2:1294964/1‑148 (MQ=255)
                       aCCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCa                                                                                                                        >  2:1112235/1‑148 (MQ=255)
                             gACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                  >  1:1522072/1‑148 (MQ=255)
                             gACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                  >  1:1521516/1‑148 (MQ=255)
                               cATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                <  1:1156617/148‑1 (MQ=255)
                               cATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                >  2:1433552/1‑148 (MQ=255)
                                aTAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                               >  1:2129956/1‑148 (MQ=255)
                                aTAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                               >  1:540123/1‑148 (MQ=255)
                                  aGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                             >  2:179399/1‑148 (MQ=255)
                                    cACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATACCTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                           <  1:110923/148‑1 (MQ=255)
                                     aCCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                          <  1:2186891/148‑1 (MQ=255)
                                        tCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                                                                                                       <  1:126323/148‑1 (MQ=255)
                                           gAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaaca                                                                                                    <  1:31597/148‑1 (MQ=255)
                                           gAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaaca                                                                                                    <  1:33089/148‑1 (MQ=255)
                                           gAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaaca                                                                                                    >  1:1992654/1‑148 (MQ=255)
                                            aaTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATACATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa                                                                                                   <  1:590454/148‑1 (MQ=255)
                                            aaTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa                                                                                                   >  2:3282084/1‑148 (MQ=255)
                                            aaTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa                                                                                                   <  1:945745/148‑1 (MQ=255)
                                            aaTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa                                                                                                   >  1:998796/1‑148 (MQ=255)
                                              ttAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg                                                                                                 >  2:1460931/1‑148 (MQ=255)
                                              ttAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg                                                                                                 >  2:976504/1‑148 (MQ=255)
                                                   aaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCAc                                                                                            <  2:2024822/148‑1 (MQ=255)
                                                      ggCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCATGAACAACGCGCACTTa                                                                                         >  2:3459795/1‑148 (MQ=255)
                                                      ggCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTa                                                                                         >  1:703917/1‑148 (MQ=255)
                                                        cGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGt                                                                                       <  1:2806450/148‑1 (MQ=255)
                                                        cGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGt                                                                                       <  1:2469339/148‑1 (MQ=255)
                                                        cGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGt                                                                                       <  2:2733882/148‑1 (MQ=255)
                                                           ggTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaa                                                                                    >  2:2885169/1‑148 (MQ=255)
                                                            gTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaat                                                                                   >  1:2544554/1‑148 (MQ=255)
                                                              tATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaataa                                                                                 >  1:2924433/1‑148 (MQ=255)
                                                              tATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaataa                                                                                 <  2:1286191/148‑1 (MQ=255)
                                                                tCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATc                                                                               >  2:1167321/1‑148 (MQ=255)
                                                                tCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATc                                                                               >  2:3233305/1‑148 (MQ=255)
                                                                  ccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAg                                                                             >  1:209159/1‑148 (MQ=255)
                                                                  ccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAg                                                                             >  2:822644/1‑148 (MQ=255)
                                                                   cTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAgg                                                                            <  2:1656612/148‑1 (MQ=255)
                                                                   cTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAgg                                                                            <  2:2130521/148‑1 (MQ=255)
                                                                    tCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGa                                                                           <  2:3145617/148‑1 (MQ=255)
                                                                    tCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGa                                                                           >  2:255009/1‑148 (MQ=255)
                                                                      aCCCTAATCGCATAAATGCCTGATGAGCTACTATTATCAGGCCTACATAACTCCTTCAATATTTTGGATTTGAATTATTTGTAGGCCGTATAAGTCTTTAACGCCGAATCCGGAATGAACATCGCGCACTTAGTCAATAATCAGGAtg                                                                         <  1:1241040/148‑1 (MQ=255)
                                                                          tAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAg                                                                     >  1:876674/1‑148 (MQ=255)
                                                                          tAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAg                                                                     <  2:2770432/148‑1 (MQ=255)
                                                                          tAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAg                                                                     >  1:2692353/1‑148 (MQ=255)
                                                                            aaCGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGt                                                                   <  2:1823321/148‑1 (MQ=255)
                                                                             aCGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGtt                                                                  <  2:586019/148‑1 (MQ=255)
                                                                                   aaaTGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCt                                                            <  1:1016697/148‑1 (MQ=255)
                                                                                         cTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCt                                                      >  1:84879/1‑148 (MQ=255)
                                                                                          tGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTg                                                     <  1:1217738/148‑1 (MQ=255)
                                                                                           gATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGt                                                    <  1:2331686/148‑1 (MQ=255)
                                                                                             tGCGCTACGTTTATCAGGCCTACATAACTCCTGCAATCTTTGGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCa                                                  <  2:1185527/148‑1 (MQ=255)
                                                                                             tGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCa                                                  <  2:484823/148‑1 (MQ=255)
                                                                                              gcgcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCat                                                 <  1:344495/148‑1 (MQ=255)
                                                                                              gcgcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCat                                                 >  2:1357942/1‑148 (MQ=255)
                                                                                                   aCGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCg                                            >  2:1163434/1‑148 (MQ=255)
                                                                                                   aCGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCg                                            >  1:1736672/1‑148 (MQ=255)
                                                                                                    cGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGt                                           >  2:2248187/1‑148 (MQ=255)
                                                                                                       ttATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaaa                                        <  1:1816782/148‑1 (MQ=255)
                                                                                                        tATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaaaa                                       <  1:3254250/148‑1 (MQ=255)
                                                                                                         aTCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaaaaa                                      <  2:1532505/148‑1 (MQ=255)
                                                                                                          tCAGGCCTACATAACTCCTGCAATATTTTGGATTTTCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAAc                                     >  1:1225007/1‑148 (MQ=255)
                                                                                                          tCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAAc                                     >  1:106264/1‑148 (MQ=255)
                                                                                                           cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACg                                    <  1:2333051/148‑1 (MQ=255)
                                                                                                             ggCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                  <  2:1506386/148‑1 (MQ=255)
                                                                                                             ggCCTACATAACTCATGCAATATTTTGGATAAGCATTATTAGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCAGAAAAAACGGc                                  >  2:442939/1‑148 (MQ=37)
                                                                                                              gcatACATAACTCGTGCAATATTTTGTATTTGGATTTTTTGTAGGACGTATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGct                                 <  2:3330955/145‑2 (MQ=37)
                                                                                                                cTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATa                               <  1:2792523/148‑1 (MQ=255)
                                                                                                                 tACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATAc                              <  1:1952744/148‑1 (MQ=255)
                                                                                                                 tACATAACTCCTGCAATAGTGTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATAc                              <  2:1949933/148‑1 (MQ=255)
                                                                                                                    aTAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGa                           <  1:179399/148‑1 (MQ=255)
                                                                                                                    aTAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGa                           <  2:1798988/148‑1 (MQ=255)
                                                                                                                    aTAACTCCTGCAATATTTTGGATGTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGa                           <  2:397714/148‑1 (MQ=255)
                                                                                                                     tAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGaa                          >  2:1866773/1‑148 (MQ=255)
                                                                                                                       aCTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAata                        >  1:1033930/1‑148 (MQ=255)
                                                                                                                        cTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAatat                       >  2:1105218/1‑148 (MQ=255)
                                                                                                                          ccTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATca                     >  1:297600/1‑148 (MQ=255)
                                                                                                                            tGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATcaca                   <  2:1261861/148‑1 (MQ=255)
                                                                                                                                atatTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAg               >  2:3207498/1‑148 (MQ=255)
                                                                                                                                   ttttGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCg            <  2:3024032/148‑1 (MQ=255)
                                                                                                                                   ttttGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCg            <  2:2604983/148‑1 (MQ=255)
                                                                                                                                    tttGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGc           >  2:158953/1‑148 (MQ=255)
                                                                                                                                        gATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTaaa       <  2:1758246/148‑1 (MQ=255)
                                                                                                                                        gATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTaaa       >  2:782378/1‑148 (MQ=255)
                                                                                                                                        gATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTaaa       >  1:98210/1‑148 (MQ=255)
                                                                                                                                         aTTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAg      <  2:3045/148‑1 (MQ=255)
                                                                                                                                         aTTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAg      <  1:890340/148‑1 (MQ=255)
                                                                                                                                          tttGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGa     >  2:1339072/1‑148 (MQ=255)
                                                                                                                                           ttGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGAc    >  2:1445721/1‑148 (MQ=255)
                                                                                                                                             gCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGt  <  2:365846/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCTGACGCCATGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGT  >  CP009273/2939168‑2939456

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: