Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 690,099 G→T 6.0% P105T (CCG→ACG)  miaB ← tRNA‑i(6)A37 methylthiotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273690,0990GT 6.0% 876.2 / 22.3 232P105T (CCG→ACG) miaBtRNA‑i(6)A37 methylthiotransferase
Reads supporting (aligned to +/- strand):  ref base G (104/114);  new base T (9/5);  total (113/119)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 233 total aligned reads displayed.

CATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTT  >  CP009273/689950‑690242
                                                                                                                                                     |                                                                                                                                               
cATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCgg                                                                                                                                                 >  2:1095126/1‑150 (MQ=255)
 atgatgGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGc                                                                                                                                                >  2:1291219/1‑150 (MQ=255)
 atgatgGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGc                                                                                                                                                >  2:3446054/1‑150 (MQ=255)
   gatgGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGTccc                                                                                                                                              >  2:568211/1‑150 (MQ=255)
        agaCAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGTCCCaaaaa                                                                                                                                         >  1:2367280/1‑150 (MQ=255)
        agaCAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCaaaaa                                                                                                                                         >  2:1123278/1‑150 (MQ=255)
         gaCAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAaat                                                                                                                                        <  1:2050126/150‑1 (MQ=255)
            aaaCGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAaataat                                                                                                                                     >  1:3122638/1‑150 (MQ=255)
             aaCGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGTCCCAAAAataata                                                                                                                                    >  2:2276793/1‑150 (MQ=255)
             aaCGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAataata                                                                                                                                    >  1:1008542/1‑150 (MQ=255)
              aCGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAatat                                                                                                                                   <  1:1445188/150‑1 (MQ=255)
                 cgGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGTCCCAAAAATAATATCGa                                                                                                                                <  1:3129554/150‑1 (MQ=255)
                 cgGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGa                                                                                                                                >  1:516049/1‑150 (MQ=255)
                 cgGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGa                                                                                                                                >  1:704469/1‑150 (MQ=255)
                 cgGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGa                                                                                                                                >  1:2807568/1‑150 (MQ=255)
                  ggTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGAc                                                                                                                               >  2:2155259/1‑150 (MQ=255)
                    tCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACAt                                                                                                                             >  2:1265642/1‑150 (MQ=255)
                    tCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACAt                                                                                                                             <  1:763791/150‑1 (MQ=255)
                     cGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATa                                                                                                                            >  1:1673489/1‑150 (MQ=255)
                     cGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATAACGACATa                                                                                                                            >  2:2918916/1‑150 (MQ=255)
                      ggCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGATGCTGATATCTACTACCGGGATGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCTGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAg                                                                                                                           <  1:1867873/150‑1 (MQ=255)
                       gCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGt                                                                                                                          >  1:1147115/1‑150 (MQ=255)
                       gCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGt                                                                                                                          <  1:3125074/150‑1 (MQ=255)
                       gCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGt                                                                                                                          >  1:1504877/1‑150 (MQ=255)
                              ggCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGgcgcg                                                                                                                   <  2:486179/150‑1 (MQ=255)
                                cgcgcgGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCt                                                                                                                 >  2:4237369/1‑150 (MQ=255)
                                     ggTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGa                                                                                                            >  2:88944/1‑150 (MQ=255)
                                     ggTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGa                                                                                                            >  2:89071/1‑150 (MQ=255)
                                     ggTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGa                                                                                                            <  1:3320832/150‑1 (MQ=255)
                                     ggTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGa                                                                                                            <  1:1123278/150‑1 (MQ=255)
                                       ttCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTTCAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAAt                                                                                                          <  1:4413311/150‑1 (MQ=255)
                                       ttCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGTCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAAt                                                                                                          >  1:261525/1‑150 (MQ=255)
                                          cGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAAtgtg                                                                                                       >  1:1528574/1‑150 (MQ=255)
                                            gTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCt                                                                                                     <  2:1759673/150‑1 (MQ=255)
                                             tAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTc                                                                                                    <  2:3382868/150‑1 (MQ=255)
                                                  ggTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGTCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCctt                                                                                               <  1:3853692/150‑1 (MQ=255)
                                                  ggTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCctt                                                                                               >  1:348617/1‑150 (MQ=255)
                                                  ggTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCctt                                                                                               <  1:182203/150‑1 (MQ=255)
                                                  ggTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCctt                                                                                               <  1:181740/150‑1 (MQ=255)
                                                     cAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCcttctt                                                                                            <  2:1937862/150‑1 (MQ=255)
                                                     cAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCcttctt                                                                                            <  2:1296613/150‑1 (MQ=255)
                                                            ctcGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGc                                                                                     <  1:1787142/150‑1 (MQ=255)
                                                              cGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATCTCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCa                                                                                   >  2:533517/1‑150 (MQ=255)
                                                              cGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCa                                                                                   >  2:1944258/1‑150 (MQ=255)
                                                                aTTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACg                                                                                 >  2:2628664/1‑150 (MQ=255)
                                                                   tCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAg                                                                              >  2:2682191/1‑150 (MQ=255)
                                                                      gggAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCa                                                                           >  2:4114255/1‑150 (MQ=255)
                                                                         aaGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCg                                                                        <  1:4082636/150‑1 (MQ=255)
                                                                         aaGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCg                                                                        >  2:762366/1‑150 (MQ=255)
                                                                         aaGCTGATATCAACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCACGCCTTCTTGCTATGCCACGCAGCCACCg                                                                        <  2:656126/150‑1 (MQ=255)
                                                                           gCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGAc                                                                      <  1:2638892/150‑1 (MQ=255)
                                                                            cTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACg                                                                     <  2:4544061/150‑1 (MQ=255)
                                                                            cTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATAGCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCGTGCGATGCCACGCAGCCACCGACg                                                                     <  2:3932129/150‑1 (MQ=255)
                                                                                 atCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGAt                                                                >  2:3333770/1‑150 (MQ=255)
                                                                                 atCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGAt                                                                <  1:2986134/150‑1 (MQ=255)
                                                                                 atCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACAGAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGAt                                                                <  2:3225830/150‑1 (MQ=255)
                                                                                  tCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATa                                                               <  2:3776708/150‑1 (MQ=255)
                                                                                        aCCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAgg                                                         <  1:2688050/150‑1 (MQ=255)
                                                                                        aCCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAgg                                                         <  1:394067/150‑1 (MQ=255)
                                                                                           gggCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCt                                                      <  1:2133407/150‑1 (MQ=255)
                                                                                            ggCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTg                                                     >  1:4122271/1‑150 (MQ=255)
                                                                                              cTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGa                                                   >  1:4306706/1‑150 (MQ=255)
                                                                                              cTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGa                                                   <  2:884656/150‑1 (MQ=255)
                                                                                                  ggTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGAttct                                               <  2:4456726/150‑1 (MQ=255)
                                                                                                  ggTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGAttct                                               >  1:1062204/1‑150 (MQ=255)
                                                                                                   gTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGAttctt                                              <  2:3099028/150‑1 (MQ=255)
                                                                                                   gTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGAttctt                                              <  2:2985586/150‑1 (MQ=255)
                                                                                                     cGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATtcttct                                            <  2:2491165/150‑1 (MQ=255)
                                                                                                      gCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTtctc                                           >  1:970789/1‑150 (MQ=255)
                                                                                                       ccACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTctct                                          <  1:4358532/150‑1 (MQ=255)
                                                                                                       ccACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTctct                                          <  1:1997779/150‑1 (MQ=255)
                                                                                                        cacgcacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCtt                                         >  1:218152/1‑150 (MQ=255)
                                                                                                        cacgcacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCtt                                         >  2:3918759/1‑150 (MQ=255)
                                                                                                         actcacgGAGTTGCTCATCTCCGGCAGTCGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGTTGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCttt                                        <  2:2552736/147‑1 (MQ=255)
                                                                                                         acgcacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCttt                                        <  1:4470178/150‑1 (MQ=255)
                                                                                                           gcacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTa                                      <  1:3783618/150‑1 (MQ=255)
                                                                                                           gcacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTa                                      >  2:2739370/1‑150 (MQ=255)
                                                                                                           gcacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTa                                      >  1:2582066/1‑150 (MQ=255)
                                                                                                            cacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTaa                                     >  2:1420456/1‑150 (MQ=255)
                                                                                                            cacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTaa                                     <  1:998778/150‑1 (MQ=255)
                                                                                                               ggAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAg                                  >  2:1554514/1‑150 (MQ=255)
                                                                                                                      atcatcTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAg                           >  1:3285190/1‑150 (MQ=255)
                                                                                                                        catcTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCg                         <  2:2823992/150‑1 (MQ=255)
                                                                                                                        catcTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCg                         >  1:3220940/1‑150 (MQ=255)
                                                                                                                         atcTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGa                        <  2:2631635/150‑1 (MQ=255)
                                                                                                                          tctcCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGAc                       >  2:2093427/1‑150 (MQ=255)
                                                                                                                          tctcCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACAGACGCCGATAATCAGGTCAGGATTCTTCTCTTTTAACATATTCCAGCGAc                       >  2:4299019/1‑150 (MQ=255)
                                                                                                                           ctcCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCGACCGACGCCGTTAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGAcc                      <  2:790574/150‑1 (MQ=255)
                                                                                                                           ctcCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGAcc                      >  2:4514055/1‑150 (MQ=255)
                                                                                                                              cGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCaa                   >  2:527370/1‑150 (MQ=255)
                                                                                                                                gCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACt                 >  2:1027711/1‑150 (MQ=255)
                                                                                                                                 cAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTg                <  2:1248441/150‑1 (MQ=255)
                                                                                                                                 cAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTg                >  2:384581/1‑150 (MQ=255)
                                                                                                                                 cAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTg                >  2:384527/1‑150 (MQ=255)
                                                                                                                                 cAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTg                <  2:1504877/150‑1 (MQ=255)
                                                                                                                                    aCGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATg             <  1:3745648/150‑1 (MQ=255)
                                                                                                                                     cGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATgg            <  2:1757302/150‑1 (MQ=255)
                                                                                                                                           aGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACt      >  2:2511902/1‑150 (MQ=255)
                                                                                                                                              gTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACtttt   <  1:2471848/150‑1 (MQ=255)
                                                                                                                                               tttGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACttttt  <  2:173258/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTT  >  CP009273/689950‑690242

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: