Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,939,313 G→T 74.7% 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,3130GT74.7% 444.6 / 257.5 325intergenic (‑112/+127)tcdA/mltAtRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA/membrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base G (48/34);  new base T (129/114);  total (177/148)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

Only 100 of 326 total aligned reads displayed.

ACGCTGACGCCATGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTG  >  CP009273/2939166‑2939462
                                                                                                                                                   |                                                                                                                                                     
  gCTGACGCCATGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTAttt                                                                                                                                                   >  1:1587238/1‑150 (MQ=255)
    tGACGCCATGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGt                                                                                                                                                 >  2:4965401/1‑150 (MQ=255)
     gACGCCATGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTa                                                                                                                                                >  2:120237/1‑150 (MQ=255)
         ccATGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGcc                                                                                                                                            >  2:1019496/1‑150 (MQ=255)
            tGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCACTATTTTGGATTTGCATGATTTGTAGGCCGGa                                                                                                                                         >  1:4352836/1‑150 (MQ=255)
            tGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGa                                                                                                                                         >  2:741268/1‑150 (MQ=255)
            tGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGa                                                                                                                                         >  2:3829507/1‑150 (MQ=255)
             gCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGAt                                                                                                                                        >  2:18850/1‑150 (MQ=255)
              cATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATa                                                                                                                                       >  2:133504/1‑150 (MQ=255)
              cATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATa                                                                                                                                       <  1:1505909/150‑1 (MQ=255)
                tCACTAATGACCCCAGACATAGCACCTCAGAATTAAACAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAg                                                                                                                                     <  2:1540009/150‑1 (MQ=255)
                 cACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAACTGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAgg                                                                                                                                    <  2:2994162/150‑1 (MQ=255)
                 cACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAgg                                                                                                                                    <  1:3231445/150‑1 (MQ=255)
                 cACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAgg                                                                                                                                    <  2:1841096/150‑1 (MQ=255)
                   cTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCg                                                                                                                                  <  1:3143613/150‑1 (MQ=255)
                    tAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGt                                                                                                                                 >  2:4434869/1‑150 (MQ=255)
                        tACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCAc                                                                                                                             <  2:3752276/150‑1 (MQ=255)
                        tACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATACCTCACCCTAAACGCATAAATGCCTGATGCGCTACTCTTATCAGGCCTACATAACTCCTGCAATATTTAGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCaa                                                                                                                             >  2:2397916/1‑149 (MQ=255)
                         aCCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                            >  1:1382196/1‑150 (MQ=255)
                             caGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                        <  2:4029033/150‑1 (MQ=255)
                                 cATAGCACCTCAGACTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                    <  2:1319525/150‑1 (MQ=255)
                                 cATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATACATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                    <  2:55256/150‑1 (MQ=255)
                                    aGCACCTCCGACTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATACATGCCTGATGCGCGACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTCTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                 <  2:617576/150‑1 (MQ=25)
                                     gCACCTCAGAATTAAAAAGGGCGATGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                <  1:1267013/150‑1 (MQ=255)
                                     gCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                >  2:3903345/1‑150 (MQ=255)
                                         cTCAGAAATAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa                                                                                                            <  2:4518052/150‑1 (MQ=255)
                                              aaTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg                                                                                                       <  2:4339042/150‑1 (MQ=255)
                                               aTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcgc                                                                                                      <  1:935070/150‑1 (MQ=255)
                                                   aaaaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCAc                                                                                                  <  2:2502398/150‑1 (MQ=255)
                                                         gCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAACAACGCGCACTTAGTc                                                                                            <  2:3565462/150‑1 (MQ=255)
                                                                   cccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAgg                                                                                  >  1:4108938/1‑150 (MQ=255)
                                                                     cTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAt                                                                                <  2:2277084/150‑1 (MQ=255)
                                                                     cTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAt                                                                                >  1:5049968/1‑150 (MQ=255)
                                                                      tCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAtg                                                                               >  1:699966/1‑150 (MQ=255)
                                                                       cACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAtgt                                                                              >  1:42885/1‑150 (MQ=255)
                                                                        aCCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAtgtg                                                                             >  1:26868/1‑150 (MQ=255)
                                                                          ccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAg                                                                           >  1:1802502/1‑150 (MQ=255)
                                                                              aaCGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTa                                                                       <  1:2321424/150‑1 (MQ=255)
                                                                              aaCGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTa                                                                       >  2:298305/1‑150 (MQ=255)
                                                                                cGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcc                                                                     >  1:2111828/1‑150 (MQ=255)
                                                                                cGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcc                                                                     >  1:29933/1‑150 (MQ=255)
                                                                                   aTAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCt                                                                  <  1:4507997/150‑1 (MQ=255)
                                                                                   aTAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCt                                                                  >  1:4057470/1‑150 (MQ=255)
                                                                                    tAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCCTTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTc                                                                 >  2:3080690/1‑150 (MQ=255)
                                                                                    tAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACACCTc                                                                 >  1:2614516/1‑150 (MQ=255)
                                                                                       aTGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAcc                                                              >  1:4845395/1‑150 (MQ=255)
                                                                                        tGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAccc                                                             <  1:3144013/150‑1 (MQ=255)
                                                                                         gCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCt                                                            >  1:4256904/1‑150 (MQ=255)
                                                                                         gCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCt                                                            >  1:1127332/1‑150 (MQ=255)
                                                                                         gCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCt                                                            <  2:144618/150‑1 (MQ=255)
                                                                                          ccTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTg                                                           >  2:1175526/1‑150 (MQ=255)
                                                                                            tGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGta                                                         <  2:2942333/150‑2 (MQ=255)
                                                                                            tGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTc                                                         >  1:406988/1‑150 (MQ=255)
                                                                                             gATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCa                                                        <  1:4900963/150‑1 (MQ=255)
                                                                                                gcgcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCatat                                                     >  1:528922/1‑150 (MQ=255)
                                                                                                gcgcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCatat                                                     <  2:3160208/150‑1 (MQ=255)
                                                                                                    tACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGt                                                 >  1:2534998/1‑150 (MQ=255)
                                                                                                    tACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGt                                                 >  1:5102353/1‑150 (MQ=255)
                                                                                                      cGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaa                                               >  1:2961954/1‑150 (MQ=255)
                                                                                                      cGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaa                                               >  1:4585089/1‑150 (MQ=255)
                                                                                                      cGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaa                                               <  1:903183/150‑1 (MQ=255)
                                                                                                       gCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaaa                                              >  1:1976783/1‑150 (MQ=255)
                                                                                                         ttATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaaaaa                                            >  2:2705084/1‑150 (MQ=255)
                                                                                                         ttATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaaaaa                                            >  2:4954395/1‑150 (MQ=255)
                                                                                                          tATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAAc                                           >  1:1955713/1‑150 (MQ=255)
                                                                                                           aTCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGGCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATCTCCGTAAAAACg                                          >  2:1682398/1‑150 (MQ=255)
                                                                                                             cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                        >  2:2150673/1‑150 (MQ=255)
                                                                                                             cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                        >  2:3964174/1‑150 (MQ=255)
                                                                                                              aGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCa                                       <  1:1035045/150‑1 (MQ=255)
                                                                                                              aGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCa                                       >  2:4802512/1‑150 (MQ=255)
                                                                                                              aGGCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCa                                       <  1:2764958/150‑1 (MQ=255)
                                                                                                               ggCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCAt                                      <  2:803371/150‑1 (MQ=255)
                                                                                                               ggCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCAt                                      >  1:2862191/1‑150 (MQ=255)
                                                                                                                gCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATa                                     >  2:2021149/1‑150 (MQ=255)
                                                                                                                gCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATa                                     <  2:3341522/150‑1 (MQ=255)
                                                                                                                gCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATa                                     <  2:1046694/150‑1 (MQ=255)
                                                                                                                gCCTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAACCCTCACCCTGTCATATCCGTAAAAACGGCATa                                     <  2:1938645/150‑1 (MQ=255)
                                                                                                                  cTACATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACa                                   <  1:3829507/150‑1 (MQ=255)
                                                                                                                     cATAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGaa                                >  1:3347378/1‑150 (MQ=255)
                                                                                                                     cATAACTACTGCAATATTTTGGATTTGCATTATTTGTAGTCCGGATAAGGCTTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGaa                                <  1:1022139/150‑1 (MQ=255)
                                                                                                                      aTAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAat                               <  2:1587238/150‑1 (MQ=255)
                                                                                                                       tAACTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAata                              <  1:1850428/150‑1 (MQ=255)
                                                                                                                       gaccTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAata                              <  2:4237246/147‑1 (MQ=255)
                                                                                                                        aaCTCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAatat                             >  2:2880694/1‑150 (MQ=255)
                                                                                                                           tCCTGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATcac                          >  1:3678910/1‑150 (MQ=255)
                                                                                                                              tGCAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATc                       >  2:1611318/1‑150 (MQ=255)
                                                                                                                                cAATATTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAg                     <  1:168993/150‑1 (MQ=255)
                                                                                                                                cAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAg                     <  1:4965388/150‑1 (MQ=255)
                                                                                                                                  atatTTTGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGcc                   <  1:2392851/150‑1 (MQ=255)
                                                                                                                                    atTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGc                 <  2:2590134/150‑1 (MQ=255)
                                                                                                                                     ttttGGATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGca                <  1:1178534/150‑2 (MQ=255)
                                                                                                                                     ttttGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCt                <  1:18651/150‑1 (MQ=255)
                                                                                                                                         ggATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAg            >  1:2867149/1‑150 (MQ=255)
                                                                                                                                          gATTTGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGTATACAGAATATCACATCAGCCGCTAAAGa           >  2:471951/1‑150 (MQ=255)
                                                                                                                                          gATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGa           >  2:5061409/1‑150 (MQ=255)
                                                                                                                                           aTTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGAc          <  1:4680726/150‑1 (MQ=255)
                                                                                                                                              tGCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGtt       <  1:4209052/150‑1 (MQ=255)
                                                                                                                                               gCATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTa      <  2:1263866/150‑1 (MQ=255)
                                                                                                                                                cATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTAc     >  1:3962450/1‑150 (MQ=255)
                                                                                                                                                  tGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCt   <  1:4000481/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACGCTGACGCCATGCATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTG  >  CP009273/2939166‑2939462

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: