Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,525,437 T→G 5.2% intergenic (+14/+152) mrcA → / ← nudE penicillin‑binding protein 1a, murein transglycosylase and transpeptidase/adenosine nucleotide hydrolase; Ap3A/Ap2A/ADP‑ribose/NADH hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,525,4370TG 5.2% 472.3 / 11.1 154intergenic (+14/+152)mrcA/nudEpenicillin‑binding protein 1a, murein transglycosylase and transpeptidase/adenosine nucleotide hydrolase; Ap3A/Ap2A/ADP‑ribose/NADH hydrolase
Reads supporting (aligned to +/- strand):  ref base T (73/73);  new base G (5/3);  total (78/76)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

Only 100 of 154 total aligned reads displayed.

GCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCAC  >  NC_000913/3525290‑3525584
                                                                                                                                                   |                                                                                                                                                   
gcaCCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCg                                                                                                                                                   >  2:608626/1‑150 (MQ=255)
   ccGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATACTGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGgcg                                                                                                                                                <  1:1133674/150‑1 (MQ=255)
        cAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGACGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCttt                                                                                                                                           <  1:1807520/150‑1 (MQ=255)
        cAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCttt                                                                                                                                           <  1:1864575/150‑1 (MQ=255)
          gTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTc                                                                                                                                         >  1:2100071/1‑150 (MQ=255)
            tAGCTAATGGTGGCATCACCCGCGAAGAGTATTTCATCGTAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTTGGTACGACCTTTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAg                                                                                                                                       <  1:898950/150‑1 (MQ=255)
            tAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGGGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAg                                                                                                                                       <  1:213344/150‑1 (MQ=255)
              gCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGtt                                                                                                                                     >  1:1007550/1‑150 (MQ=255)
              gCTAATGGTGCCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACATAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAAATGTTCTGATTAAAAATGCGCTTCGGCGCCTTTTCAGtt                                                                                                                                     >  2:1593646/1‑150 (MQ=255)
               cTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGCGGTGGGAACGACCATTATCGATTATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGttt                                                                                                                                    >  2:142005/1‑150 (MQ=255)
                tAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTg                                                                                                                                   <  1:682765/150‑1 (MQ=255)
                tAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTg                                                                                                                                   <  1:1186901/150‑1 (MQ=255)
                  atgtgggCAACAGCCGCGAAGAGTATTTCATGGAAGGTACTCATGCTACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCt                                                                                                                                 <  2:440972/145‑1 (MQ=255)
                   tggtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTg                                                                                                                                <  1:1048864/150‑1 (MQ=255)
                    ggtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAATGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGa                                                                                                                               <  2:510411/150‑1 (MQ=255)
                     gtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTTCACCAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                                                                              >  2:41560/1‑150 (MQ=255)
                     gtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                                                                              <  2:1208998/150‑1 (MQ=255)
                            cAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                       <  1:515947/150‑1 (MQ=255)
                               ccGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACt                                                                                                                    >  2:704358/1‑150 (MQ=255)
                                cgcgAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACtt                                                                                                                   >  2:1029193/1‑150 (MQ=255)
                                    aaGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGttt                                                                                                               <  2:2157770/150‑1 (MQ=255)
                                     agagTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTa                                                                                                              <  1:421714/150‑1 (MQ=255)
                                       agTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGAACTTGTTTAtt                                                                                                            >  1:1938779/1‑149 (MQ=255)
                                        gTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGc                                                                                                           >  1:265228/1‑150 (MQ=255)
                                         tATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGAGCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGcc                                                                                                          <  2:2076068/150‑1 (MQ=255)
                                           tttCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCgg                                                                                                        <  2:886000/150‑1 (MQ=255)
                                           tttCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCgg                                                                                                        <  2:1018942/150‑1 (MQ=255)
                                              cATCGAAGGTACGCCGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATa                                                                                                     >  2:1131614/1‑150 (MQ=255)
                                              cATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATa                                                                                                     >  2:1183701/1‑150 (MQ=255)
                                                tCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACg                                                                                                   >  1:94629/1‑150 (MQ=255)
                                                  gAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGc                                                                                                 >  1:833230/1‑150 (MQ=255)
                                                  gAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGc                                                                                                 <  1:1099503/150‑1 (MQ=255)
                                                   aaGGCACGCCGCCGACACAACAGGCAGAGCACGAGGTGGGAACGACCAATATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAACAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCg                                                                                                <  1:332203/150‑1 (MQ=255)
                                                     gggaCGCCGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTg                                                                                              <  1:549591/147‑1 (MQ=255)
                                                      gTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGa                                                                                             <  1:1098909/150‑1 (MQ=255)
                                                       tACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGaa                                                                                            >  2:1866418/1‑150 (MQ=255)
                                                        aCGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAAc                                                                                           >  2:2185251/1‑150 (MQ=255)
                                                          gCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGAACTTGTTTATGCCGGATACGGCGTGAAcgc                                                                                         >  2:174614/1‑150 (MQ=255)
                                                           cAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAAcgcg                                                                                        >  2:2102968/1‑150 (MQ=255)
                                                           cAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAAcgcg                                                                                        <  1:1083167/150‑1 (MQ=255)
                                                           cAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAAcgcg                                                                                        >  1:1734601/1‑150 (MQ=255)
                                                            aGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGt                                                                                       <  2:488367/150‑1 (MQ=255)
                                                             gCCGACACCACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCTCAGGAATTGTTCTGATTAAAACGGCGCGTCGGCGCCTTTTGCGTTTGCGGACAAAGTGGACTTGTTTATGCCGGATACGGCGTGAACGCGtt                                                                                      >  1:1459180/1‑150 (MQ=255)
                                                              cagACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATAACGCGTTCAcgccgta                                                                                     <  2:1131162/148‑5 (MQ=255)
                                                               cGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGACTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTAt                                                                                    >  1:1614276/1‑150 (MQ=255)
                                                                gACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCTAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGTCGCCTTTTCAGTTTGCAGAAAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGtttac                                                                                   >  1:186146/1‑147 (MQ=255)
                                                                  cacaACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCg                                                                                 <  1:1567625/150‑1 (MQ=255)
                                                                    caacaGGCAGTGCACGAGGTGGGAACGACCATAATCGATAATGGCGAGGCACAGGAATTCTTCTGATTAAAAAGGCGCATCGGCGCCTTATCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGTATACGGCGTGAACGAGTAATCCGGc                                                                               >  1:1998049/1‑150 (MQ=255)
                                                                     aacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCTCTTCGGCGCGTTTTGAGTTTGCTGACATAGTGAACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCt                                                                              <  2:1059205/150‑1 (MQ=255)
                                                                       caGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTaa                                                                            >  1:1382961/1‑150 (MQ=255)
                                                                         ggCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACa                                                                          >  1:585/1‑150 (MQ=255)
                                                                         ggCAGTGCACGAGGAGGGAACGACCATTATCTATAATGGCGAGGCACATGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTAAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCTGGAACGCTTTATACGTCTaaaa                                                                          >  2:856485/1‑148 (MQ=255)
                                                                           cAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACaaa                                                                        >  2:797013/1‑150 (MQ=255)
                                                                               gCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACATGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCg                                                                    <  2:1978558/150‑1 (MQ=255)
                                                                               gCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGACGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCg                                                                    >  2:1274112/1‑150 (MQ=255)
                                                                               gCACGAGGTGGGAACGACCATTATCGATAATGGAGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTGCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCg                                                                    >  2:1057991/1‑150 (MQ=255)
                                                                                      gTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAAc                                                             <  2:2203723/150‑1 (MQ=255)
                                                                                        gggAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTc                                                           >  1:1740908/1‑150 (MQ=255)
                                                                                        gggAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTc                                                           >  1:646967/1‑150 (MQ=255)
                                                                                          gAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCaa                                                         >  2:1961348/1‑150 (MQ=255)
                                                                                          gAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCaa                                                         >  2:2214510/1‑150 (MQ=255)
                                                                                            aCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATa                                                       <  2:2068986/150‑1 (MQ=255)
                                                                                             cGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATaa                                                      >  2:1887238/1‑150 (MQ=255)
                                                                                               aCCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAAt                                                    >  2:942550/1‑150 (MQ=255)
                                                                                               aCCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAAt                                                    >  1:2025469/1‑150 (MQ=255)
                                                                                                ccATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAATGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAAtt                                                   >  1:379313/1‑150 (MQ=255)
                                                                                                ccATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAAtt                                                   >  1:943289/1‑150 (MQ=255)
                                                                                                    tATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAg                                               >  2:1811876/1‑150 (MQ=255)
                                                                                                       cGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTAAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGaaa                                            >  1:1253024/1‑150 (MQ=255)
                                                                                                         aTAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAAcc                                          <  1:532417/150‑1 (MQ=255)
                                                                                                            aTGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCa                                       <  1:1895744/150‑1 (MQ=255)
                                                                                                            aTGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCa                                       <  1:2129684/150‑1 (MQ=255)
                                                                                                            aTGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCa                                       <  1:1809977/150‑1 (MQ=255)
                                                                                                            aTGGCGAGGCACAGGAATTGTTCAGATTAAAAAGGCGCTTCGGCGCCTATTCAGTTTGCTGACAAAGTGCACTTGAATATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCa                                       <  1:1818578/150‑1 (MQ=255)
                                                                                                              ggCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATa                                     >  2:2123429/1‑150 (MQ=255)
                                                                                                                     cacaGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGa                              <  1:1396687/150‑1 (MQ=255)
                                                                                                                     cacaGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGa                              >  1:1150889/1‑150 (MQ=255)
                                                                                                                     cacaGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACAGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAACCCCCATAAGCCTGa                              >  1:1633120/1‑150 (MQ=255)
                                                                                                                      acaGGAATGGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGAt                             <  2:1459180/150‑1 (MQ=255)
                                                                                                                         ggAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGAGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCAGATaaa                          >  1:844105/1‑150 (MQ=255)
                                                                                                                              tGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtg                     <  1:797013/150‑1 (MQ=255)
                                                                                                                              tGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtg                     <  1:1183701/150‑1 (MQ=255)
                                                                                                                               gTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgt                    >  2:1553177/1‑150 (MQ=255)
                                                                                                                                ttCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgtg                   <  2:1734601/150‑1 (MQ=255)
                                                                                                                                ttCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgtg                   <  2:646967/150‑1 (MQ=255)
                                                                                                                                 tCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAGAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGATAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAATCATTGTGc                  <  1:109018/150‑1 (MQ=255)
                                                                                                                                 tCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGc                  <  2:2145595/150‑1 (MQ=255)
                                                                                                                                 tCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGc                  >  1:2102484/1‑150 (MQ=255)
                                                                                                                                   tGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCAt                >  1:554073/1‑150 (MQ=255)
                                                                                                                                    aactgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATc               >  1:910812/6‑150 (MQ=255)
                                                                                                                                     aTTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACGCAATAAATTGCAGAAACACCATAAGCCTGATAAACATTGTGAATCa              >  2:1116995/1‑150 (MQ=255)
                                                                                                                                      ttAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAg             <  2:1096613/150‑1 (MQ=255)
                                                                                                                                       tgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAgg            >  2:1138457/3‑150 (MQ=255)
                                                                                                                                       tgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAATCGTGAAAACTCAATAAATTGCAGAACCCCCATAAGCCTGATAAACATTGTGCATCAGGc           >  2:2179122/3‑150 (MQ=255)
                                                                                                                                        gaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGc           >  1:1921261/2‑150 (MQ=255)
                                                                                                                                            cggCGCTTCGGCGCCTTTTCAGTTAGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGCGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGGGCATCAGGCAAAc       <  2:336266/149‑1 (MQ=255)
                                                                                                                                               cgcTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAGGcttatggg    <  1:700578/150‑8 (MQ=255)
                                                                                                                                               cgcTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTc    >  1:1800458/1‑150 (MQ=255)
                                                                                                                                                gcTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCa   <  2:1267270/150‑1 (MQ=255)
                                                                                                                                                 cTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCAc  <  1:834898/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCAC  >  NC_000913/3525290‑3525584

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