Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,525,438 C→A 5.1% intergenic (+15/+151) 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,4380CA 5.1% 420.3 / 6.9 138intergenic (+15/+151)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 C (70/60);  new base A (4/4);  total (74/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 139 total aligned reads displayed.

AGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACG  >  NC_000913/3525289‑3525585
                                                                                                                                                     |                                                                                                                                                   
      gggCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCGTCGGCGcc                                                                                                                                               <  1:1437120/150‑1 (MQ=255)
       ggCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCACTTCGGCGACt                                                                                                                                              >  2:116595/1‑150 (MQ=255)
         cAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCttt                                                                                                                                            <  2:1189361/150‑1 (MQ=255)
          aGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCCGAAGCGCCtttt                                                                                                                                           <  2:141368/150‑1 (MQ=255)
           gTTAGCTAATGGTGGCAACAGACGCTAAGAGTATTTCAACGAAGGTACGCATCAGACACCACACGCAGTGCACTAGGTTGGAACTAACATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCCGAAGCGCCttttt                                                                                                                                          <  2:335902/150‑2 (MQ=255)
              aGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGt                                                                                                                                       >  1:33841/1‑150 (MQ=255)
              aGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGt                                                                                                                                       <  1:1165289/150‑1 (MQ=255)
                cTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGttt                                                                                                                                     >  2:222510/1‑150 (MQ=255)
                    tggtggCAACAGCCGCGAAGAGTATTTCCTCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGcgg                                                                                                                                 <  2:523494/150‑3 (MQ=255)
                     ggtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGCCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGa                                                                                                                                >  1:188350/1‑150 (MQ=255)
                      gtggCAACAGCCGCGAAGAGTATTTCCTCGAAGGTACGCAGCCGACACAACAGGGAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                                                                               <  2:837251/150‑1 (MQ=255)
                      gtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                                                                               <  1:935175/150‑1 (MQ=255)
                        ggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACaa                                                                                                                             <  1:1316552/150‑1 (MQ=255)
                         gCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACAAAACAGGCAGTGCACGAGGTTGGAACGACCATTATCGAAAATTTCAATGCACAGTAATTGTTCTGATTAAAAAGGAGCTACGGCGACTTTTCAGTTTGCTGACaaa                                                                                                                            >  2:863931/1‑150 (MQ=255)
                          cAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGCCACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAg                                                                                                                           >  1:696334/1‑150 (MQ=255)
                             cAGCCGCGAAGCGTATTTCATCGACGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                        <  1:1057906/150‑1 (MQ=255)
                             cAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                        <  2:557396/150‑1 (MQ=255)
                             cAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                        >  1:964688/1‑150 (MQ=255)
                               gCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCAc                                                                                                                      >  1:105398/1‑150 (MQ=255)
                                ccGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCGCGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACt                                                                                                                     >  1:1255415/1‑150 (MQ=255)
                                ccGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACt                                                                                                                     >  1:379290/1‑150 (MQ=255)
                                 cgcgAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACtt                                                                                                                    >  2:38455/1‑150 (MQ=255)
                                  gcgAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTg                                                                                                                   >  2:1063714/1‑150 (MQ=255)
                                  gcgAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTg                                                                                                                   >  1:1612312/1‑150 (MQ=255)
                                     aaGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGATAAAGTGCACTTGttt                                                                                                                <  2:1318218/150‑1 (MQ=255)
                                       gagTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTAt                                                                                                              >  2:1805/1‑150 (MQ=255)
                                          tATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGcc                                                                                                           <  1:222510/150‑1 (MQ=255)
                                          tATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGcc                                                                                                           <  1:355532/150‑1 (MQ=255)
                                              tCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGAt                                                                                                       <  2:32591/150‑1 (MQ=255)
                                               cATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATa                                                                                                      >  2:1475768/1‑150 (MQ=255)
                                                 tCGAAGGTACGCAGCCGACACAACAGGCAGAGCACGAGGTGGTCACGCCCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGAATAAAAAGGCGCTTCGGCGCCTGTTCAGTTTGCTGACAAAGTGCACTTGTATATGCAGGATACg                                                                                                    >  1:1095136/1‑150 (MQ=255)
                                                  cGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGACCGCCCATTATCGAGAATGGCGAGGCACAGGAATTGGTCTGATTAAAACGGCGCTTCGGCGCCTTGTCCGGTAGCTGCCACAGTGCACTTGTTTATGCCGGATACgc                                                                                                   >  1:334966/1‑149 (MQ=255)
                                                   gAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGc                                                                                                  >  1:379010/1‑150 (MQ=255)
                                                        tACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGaa                                                                                             >  2:1387665/1‑150 (MQ=255)
                                                          cGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAAcg                                                                                           <  1:1558157/150‑1 (MQ=255)
                                                             ctccGACACAAGATGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAACGGCGCTGCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGt                                                                                        <  1:1516355/148‑1 (MQ=255)
                                                              gCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGtt                                                                                       >  1:1027745/1‑150 (MQ=255)
                                                              gCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGtt                                                                                       >  2:681571/1‑150 (MQ=255)
                                                              gCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCGTCGGCGCGTTTTCAGTTTGCTGACAAAGTGCAATTGTTTATGCAGGATACGGCGAGAACGCGtt                                                                                       >  2:1395893/1‑150 (MQ=255)
                                                               ccGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTa                                                                                      <  1:526550/150‑1 (MQ=255)
                                                               ccGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTa                                                                                      <  1:519966/150‑1 (MQ=255)
                                                                cGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCttt                                                                                                                                            >  2:174575/1‑95 (MQ=255)
                                                                cGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCttt                                                                                                                                            <  1:174575/95‑1 (MQ=255)
                                                                cGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTAt                                                                                     >  1:1593786/1‑150 (MQ=255)
                                                                   cacaACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCg                                                                                  >  1:923505/1‑150 (MQ=255)
                                                                   cacaACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCg                                                                                  >  2:1680255/1‑150 (MQ=255)
                                                                   cacaACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCg                                                                                  <  2:1229541/150‑1 (MQ=255)
                                                                    acaacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCgg                                                                                 >  2:1402308/1‑150 (MQ=255)
                                                                     caacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGc                                                                                <  2:160488/150‑1 (MQ=255)
                                                                      atcaGGCAGTGCACGAGGTGGGAACGACCATTCTCGATAATGGCGATGCACAGGAATTGTTCTGATTAAAAAGGCGCCTAGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTTTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCt                                                                               <  2:902357/148‑1 (MQ=255)
                                                                       acaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTa                                                                              >  2:26092/1‑150 (MQ=255)
                                                                        caGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTaa                                                                             >  2:187537/1‑150 (MQ=255)
                                                                         aGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAAc                                                                            >  1:1609996/1‑150 (MQ=255)
                                                                           gCAGTCCACGAGGTGGGACCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTGCGGCGCCTTTTCAGATTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACaa                                                                          <  1:1093869/150‑1 (MQ=255)
                                                                                gCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCg                                                                     >  1:991887/1‑150 (MQ=255)
                                                                                  aCGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTg                                                                   <  1:1206354/150‑1 (MQ=255)
                                                                                    gAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGaa                                                                 <  1:247473/150‑1 (MQ=255)
                                                                                       gTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAAc                                                              <  1:187537/150‑1 (MQ=255)
                                                                                        tGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACt                                                             <  2:33841/150‑1 (MQ=255)
                                                                                          ggAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTACGGCGCCTTTTCAGTTTGCTGAAAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCa                                                           >  1:1109910/1‑150 (MQ=255)
                                                                                            cctcgccATTATCGATAATGGCGAGGCCCAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGATAACTCAAt                                                         >  2:163324/6‑150 (MQ=255)
                                                                                            cctcgccATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAAt                                                         >  2:1333614/6‑150 (MQ=255)
                                                                                                 ccATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAAtt                                                    >  1:1268623/1‑150 (MQ=255)
                                                                                                  cATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTg                                                   <  2:1191280/150‑1 (MQ=255)
                                                                                                   aTTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCTTTATCCGGCTAACAAAATCGTGAAAACTCATTAAATTGc                                                  >  1:849416/1‑150 (MQ=255)
                                                                                                   aTTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGc                                                  >  2:269064/1‑150 (MQ=255)
                                                                                                     tATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAg                                                <  1:958388/150‑1 (MQ=255)
                                                                                                     tATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAg                                                <  2:796553/150‑1 (MQ=255)
                                                                                                      aTCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGa                                               >  2:471491/1‑150 (MQ=255)
                                                                                                               ggCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATa                                      >  2:896899/1‑150 (MQ=255)
                                                                                                                gCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATaa                                     <  2:628429/150‑1 (MQ=255)
                                                                                                                   aGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGTCGACTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGcc                                  <  2:1036926/150‑1 (MQ=255)
                                                                                                                        caGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATa                             >  1:437453/1‑150 (MQ=255)
                                                                                                                          ggAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATaaa                           >  2:800423/1‑150 (MQ=255)
                                                                                                                            aaTTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACa                         >  2:396765/1‑150 (MQ=255)
                                                                                                                            aaTTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACa                         <  1:384367/150‑1 (MQ=255)
                                                                                                                                 ttCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgtg                    <  2:632672/150‑1 (MQ=255)
                                                                                                                                 ttCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgtg                    <  2:1325045/150‑1 (MQ=255)
                                                                                                                                  tCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGc                   <  2:167329/150‑1 (MQ=255)
                                                                                                                                   cTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCa                  >  2:608013/1‑150 (MQ=255)
                                                                                                                                     gATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATc                >  2:81028/1‑150 (MQ=255)
                                                                                                                                     gATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATc                <  2:1052260/150‑1 (MQ=255)
                                                                                                                                     gATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATc                >  1:693069/1‑150 (MQ=255)
                                                                                                                                     gATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCATGATAAACATTGTGCATc                >  1:81006/1‑150 (MQ=255)
                                                                                                                                       ttAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGAGACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAg              <  2:270478/150‑1 (MQ=255)
                                                                                                                                        tAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAGTCGTGCAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAgg             >  1:830139/1‑150 (MQ=255)
                                                                                                                                         gaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCCGGc            >  1:1503919/2‑150 (MQ=255)
                                                                                                                                         gaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGc            >  1:196918/2‑150 (MQ=255)
                                                                                                                                         gaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCAGACAAAGTGCACATGTTTATGCCGGATACTGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAACCCCCAAAACCCTGATAAACATTGTGCATCAGGc            >  1:667244/2‑150 (MQ=255)
                                                                                                                                         aaaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGc            <  2:701556/150‑1 (MQ=255)
                                                                                                                                         aaaaaGGAGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGc            <  2:216025/150‑1 (MQ=255)
                                                                                                                                          aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCa           >  1:487333/1‑150 (MQ=255)
                                                                                                                                          aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCa           <  2:1644524/150‑1 (MQ=255)
                                                                                                                                          aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCa           >  1:538880/1‑150 (MQ=255)
                                                                                                                                             aGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAAc        <  2:1268623/150‑1 (MQ=255)
                                                                                                                                              ggCGCTTCGGCGCCTTTTCAGTTTGCTGACAAATTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACt       >  1:990784/1‑150 (MQ=255)
                                                                                                                                              ggCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACt       >  1:918773/1‑150 (MQ=255)
                                                                                                                                              ggCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACt       >  2:1033786/1‑150 (MQ=255)
                                                                                                                                              ggCGCGTCGTCGCCTTTTCAGTTTGCTGACAAAGTTCACTTGGTTATGCCGGATACGGCGTGAAAGCGTTATCAGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCACATCAGCCTGATGCACAATGTTTCTCAGGcttatg       <  2:1095136/150‑6 (MQ=255)
                                                                                                                                                   ttCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACg  <  1:480334/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACG  >  NC_000913/3525289‑3525585

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: