Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,525,439 G→A 5.9% intergenic (+16/+150) 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,4390GA 5.9% 472.7 / 12.3 153intergenic (+16/+150)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 G (72/72);  new base A (5/4);  total (77/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

Only 100 of 155 total aligned reads displayed.

GCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAT  >  NC_000913/3525290‑3525588
                                                                                                                                                     |                                                                                                                                                     
gcaCCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCg                                                                                                                                                       >  2:608626/1‑150 (MQ=255)
   ccGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATACTGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGgcg                                                                                                                                                    <  1:1133674/150‑1 (MQ=255)
   ccGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGgcg                                                                                                                                                    <  1:552025/150‑1 (MQ=255)
          gTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTc                                                                                                                                             >  1:2100071/1‑150 (MQ=255)
            tAGCTAATGGTGGCATCACCCGCGAAGAGTATTTCATCGTAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTTGGTACGACCTTTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAg                                                                                                                                           <  1:898950/150‑1 (MQ=255)
              gCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGtt                                                                                                                                         >  1:1007550/1‑150 (MQ=255)
              gCTAATGGTGCCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACATAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAAATGTTCTGATTAAAAATGCGCTTCGGCGCCTTTTCAGtt                                                                                                                                         >  2:1593646/1‑150 (MQ=255)
                tAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTg                                                                                                                                       <  1:1186901/150‑1 (MQ=255)
                  atgtgggCAACAGCCGCGAAGAGTATTTCATGGAAGGTACTCATGCTACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCt                                                                                                                                     <  2:440972/145‑1 (MQ=255)
                  aTGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGTGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCt                                                                                                                                     <  1:957529/150‑1 (MQ=255)
                   tggtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTg                                                                                                                                    <  1:1048864/150‑1 (MQ=255)
                    ggtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAATGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGa                                                                                                                                   <  2:510411/150‑1 (MQ=255)
                     gtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                                                                                  >  1:1389927/1‑150 (MQ=255)
                     gtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                                                                                  <  2:1208998/150‑1 (MQ=255)
                     ggggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                                                                                  <  1:689731/148‑1 (MQ=255)
                         cAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAg                                                                                                                              <  2:1040925/150‑1 (MQ=255)
                            cAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                           <  1:515947/150‑1 (MQ=255)
                            cAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                           <  1:2115387/150‑1 (MQ=255)
                               ccGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACt                                                                                                                        >  2:704358/1‑150 (MQ=255)
                                    aaGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGttt                                                                                                                   <  2:2157770/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)
                                         tATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGcc                                                                                                              <  1:422555/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:1945032/1‑150 (MQ=255)
                                                  gAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGc                                                                                                     >  1:833230/1‑150 (MQ=255)
                                                  gAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGc                                                                                                     <  1:1099503/150‑1 (MQ=255)
                                                     gttACGCAGCCGACACAACAGGCAGTGCACGAGGTGGTAACGACCATTATCGATAATGGCGAGGCACATGATTTGTTGTGATTAAAATGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTTTTTATGCCGGATACGGCGTg                                                                                                  <  1:18327/148‑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)
                                                               cGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTAt                                                                                        >  1:837457/1‑150 (MQ=255)
                                                               cGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGACTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTAt                                                                                        >  1:1614276/1‑150 (MQ=255)
                                                                gACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCTAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGTCGCCTTTTCAGTTTGCAGAAAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGtttac                                                                                       >  1:186146/1‑147 (MQ=255)
                                                                  cacaACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCg                                                                                     <  1:1567625/150‑1 (MQ=255)
                                                                  cacaACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCg                                                                                     <  2:1400436/150‑1 (MQ=255)
                                                                    caacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAATGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGc                                                                                   <  2:1081255/150‑1 (MQ=255)
                                                                    caacaGGCAGTGCACGAGGTGGGAACGACCATAATCGATAATGGCGAGGCACAGGAATTCTTCTGATTAAAAAGGCGCATCGGCGCCTTATCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGTATACGGCGTGAACGAGTAATCCGGc                                                                                   >  1:1998049/1‑150 (MQ=255)
                                                                     aacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCTCTTCGGCGCGTTTTGAGTTTGCTGACATAGTGAACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCt                                                                                  <  2:1059205/150‑1 (MQ=255)
                                                                     aacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCt                                                                                  <  2:1228604/150‑1 (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)
                                                                                  cGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGa                                                                     <  2:2040408/150‑1 (MQ=255)
                                                                                    aGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGaaa                                                                   <  1:2182349/150‑1 (MQ=255)
                                                                                      gTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAAc                                                                 <  2:2203723/150‑1 (MQ=255)
                                                                                       tGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACt                                                                >  2:360620/1‑150 (MQ=255)
                                                                                        gggAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTc                                                               >  1:646967/1‑150 (MQ=255)
                                                                                        gggAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTc                                                               >  1:1740908/1‑150 (MQ=255)
                                                                                          gAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCaa                                                             >  2:2214510/1‑150 (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)
                                                                                                  aTTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGc                                                     <  2:116488/150‑1 (MQ=255)
                                                                                                    tATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAg                                                   >  2:1811876/1‑150 (MQ=255)
                                                                                                       cGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTAAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGaaa                                                >  1:1253024/1‑150 (MQ=255)
                                                                                                        gATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCTCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAAc                                               >  1:622540/1‑150 (MQ=255)
                                                                                                         aTAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAAcc                                              <  1:532417/150‑1 (MQ=255)
                                                                                                            aTGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCa                                           <  1:2129684/150‑1 (MQ=255)
                                                                                                            aTGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCa                                           <  1:1895744/150‑1 (MQ=255)
                                                                                                            aTGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCa                                           <  1:1809977/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)
                                                                                                                       caGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATa                                >  2:1486347/1‑150 (MQ=255)
                                                                                                                       caGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATa                                <  1:935559/150‑1 (MQ=255)
                                                                                                                         ggAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGAGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCAGATaaa                              >  1:844105/1‑150 (MQ=255)
                                                                                                                              tGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtg                         <  1:1183701/150‑1 (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)
                                                                                                                                   tGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCAt                    >  1:554073/1‑150 (MQ=255)
                                                                                                                                     aTTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCa                  >  1:1809218/1‑150 (MQ=255)
                                                                                                                                     aTTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACGCAATAAATTGCAGAAACACCATAAGCCTGATAAACATTGTGAATCa                  >  2:1116995/1‑150 (MQ=255)
                                                                                                                                       tgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAgg                >  2:1138457/3‑150 (MQ=255)
                                                                                                                                       tgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAATCGTGAAAACTCAATAAATTGCAGAACCCCCATAAGCCTGATAAACATTGTGCATCAGGc               >  2:2179122/3‑150 (MQ=255)
                                                                                                                                        gaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGc               >  1:1921261/2‑150 (MQ=255)
                                                                                                                                        aaaaaGGTGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGc               >  1:815639/1‑150 (MQ=255)
                                                                                                                                            cggCGCTTCGGCGCCTTTTCAGTTAGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGCGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGGGCATCAGGCAAAc           <  2:336266/149‑1 (MQ=255)
                                                                                                                                               cgcTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAGGcttatggg        <  1:700578/150‑8 (MQ=255)
                                                                                                                                               cgcTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAGGcttatggg        <  1:1866418/150‑8 (MQ=255)
                                                                                                                                               cgcTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTc        >  1:1800458/1‑150 (MQ=255)
                                                                                                                                                 cTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCAc      <  1:360620/150‑1 (MQ=255)
                                                                                                                                                     ggCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAt  >  2:1038123/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAT  >  NC_000913/3525290‑3525588

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: