Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,525,437 T→G 6.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 6.2% 384.5 / 9.0 138intergenic (+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 (69/57);  new base G (4/6);  total (75/63)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 140 total aligned reads displayed.

CAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACG  >  NC_000913/3525288‑3525585
                                                                                                                                                     |                                                                                                                                                    
       gggCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGcc                                                                                                                                               <  1:19281/150‑1 (MQ=255)
       gggCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCGTCGGCGcc                                                                                                                                               <  1:1437120/150‑1 (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)
                gCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGtt                                                                                                                                      >  1:134843/1‑150 (MQ=255)
                 cTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGttt                                                                                                                                     <  1:380493/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)
                       gtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGAGCATTATCGATAATGGCGAGGCACAGGAATTGGGCTGCTTAAAACGGCGCTTCGGCGCCGTTTCAGTTTGCTgcc                                                                                                                               >  1:270478/1‑148 (MQ=255)
                       gtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                                                                               <  1:1599444/150‑1 (MQ=255)
                         ggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACaa                                                                                                                             >  2:1122169/1‑150 (MQ=255)
                          gCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACAAAACAGGCAGTGCACGAGGTTGGAACGACCATTATCGAAAATTTCAATGCACAGTAATTGTTCTGATTAAAAAGGAGCTACGGCGACTTTTCAGTTTGCTGACaaa                                                                                                                            >  2:863931/1‑150 (MQ=255)
                           cAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGCCACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAg                                                                                                                           >  1:696334/1‑150 (MQ=255)
                              cAGCCGCGAAGCGTATTTCATCGACGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                        <  1:1057906/150‑1 (MQ=255)
                              cAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                        >  1:964688/1‑150 (MQ=255)
                              cAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                        <  1:1673956/150‑1 (MQ=255)
                              cAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                        <  2:557396/150‑1 (MQ=255)
                              cAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                                                                                        <  2:1683514/150‑1 (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)
                                        gagTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTAt                                                                                                              >  2:1805/1‑150 (MQ=255)
                                          gTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGc                                                                                                            >  1:872352/1‑150 (MQ=255)
                                           tATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGcc                                                                                                           <  1:355532/150‑1 (MQ=255)
                                            aTTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCATCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGccc                                                                                                          >  1:613034/1‑149 (MQ=255)
                                               tCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGAt                                                                                                       <  2:32591/150‑1 (MQ=255)
                                                cATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATa                                                                                                      >  2:1475768/1‑150 (MQ=255)
                                                  tCGAAGGTACGCAGCCGACACAACAGGCAGAGCACGAGGTGGTCACGCCCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGAATAAAAAGGCGCTTCGGCGCCTGTTCAGTTTGCTGACAAAGTGCACTTGTATATGCAGGATACg                                                                                                    >  1:1095136/1‑150 (MQ=255)
                                                       ggTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTg                                                                                               <  1:466270/150‑1 (MQ=255)
                                                           cGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAAcg                                                                                           <  2:1195118/150‑1 (MQ=255)
                                                             cacccGACCCAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCAGAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAAcgcg                                                                                         <  1:98033/147‑1 (MQ=255)
                                                              ctccGACACAAGATGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAACGGCGCTGCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGt                                                                                        <  1:1516355/148‑1 (MQ=255)
                                                               gCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGtt                                                                                       >  2:681571/1‑150 (MQ=255)
                                                               gCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGtt                                                                                       >  1:1027745/1‑150 (MQ=255)
                                                               gCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCGTCGGCGCGTTTTCAGTTTGCTGACAAAGTGCAATTGTTTATGCAGGATACGGCGAGAACGCGtt                                                                                       >  2:1395893/1‑150 (MQ=255)
                                                                ccGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTa                                                                                      <  1:519966/150‑1 (MQ=255)
                                                                 cGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCttt                                                                                                                                            <  1:174575/95‑1 (MQ=255)
                                                                 cGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTAt                                                                                     >  1:1593786/1‑150 (MQ=255)
                                                                    cacaACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCg                                                                                  >  2:1680255/1‑150 (MQ=255)
                                                                    cacaACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCg                                                                                  >  1:923505/1‑150 (MQ=255)
                                                                     acaacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCgg                                                                                 >  2:1402308/1‑150 (MQ=255)
                                                                      caacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGc                                                                                >  1:521397/1‑150 (MQ=255)
                                                                       atcaGGCAGTGCACGAGGTGGGAACGACCATTCTCGATAATGGCGATGCACAGGAATTGTTCTGATTAAAAAGGCGCCTAGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTTTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCt                                                                               <  2:902357/148‑1 (MQ=255)
                                                                         caGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTaa                                                                             >  2:187537/1‑150 (MQ=255)
                                                                          aGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAAc                                                                            >  1:1609996/1‑150 (MQ=255)
                                                                            gCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGAGTAAACAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACaa                                                                          <  1:355484/150‑1 (MQ=255)
                                                                            gCAGTCCACGAGGTGGGACCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTGCGGCGCCTTTTCAGATTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACaa                                                                          <  1:1093869/150‑1 (MQ=255)
                                                                             cAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCTCTTCGGCGCCTATTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACaaa                                                                         >  1:1150082/1‑150 (MQ=255)
                                                                                 gCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCg                                                                     >  1:991887/1‑150 (MQ=255)
                                                                                   aCGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTg                                                                   <  1:1206354/150‑1 (MQ=255)
                                                                                     gAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGaa                                                                 <  1:247473/150‑1 (MQ=255)
                                                                                         tGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACt                                                             <  2:33841/150‑1 (MQ=255)
                                                                                         ggggAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGGGTTATCCGGCTAACAAAATCGTGAAAACt                                                             <  1:1516316/149‑1 (MQ=255)
                                                                                             cctcgccATTATCGATAATGGCGAGGCCCAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGATAACTCAAt                                                         >  2:163324/6‑150 (MQ=255)
                                                                                             cctcgccATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAAt                                                         >  2:1333614/6‑150 (MQ=255)
                                                                                             aaCGACCATTATCGATCATTGCGAGGGACAGGAAGTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACTGGGTGAACGCGTTATCCGGCTAACTAAATCGTGAAAACTCatt                                                         <  2:1358054/150‑3 (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:701556/1‑150 (MQ=255)
                                                                                                      tATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAg                                                <  2:796553/150‑1 (MQ=255)
                                                                                                      tATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAg                                                <  1:958388/150‑1 (MQ=255)
                                                                                                       aTCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGa                                               >  2:471491/1‑150 (MQ=255)
                                                                                                            tAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGATGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATACGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAAccc                                          >  1:930200/1‑150 (MQ=255)
                                                                                                                ggCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATa                                      <  2:1505766/150‑1 (MQ=255)
                                                                                                                ggCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATa                                      >  2:896899/1‑150 (MQ=255)
                                                                                                                 gCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATaa                                     <  2:628429/150‑1 (MQ=255)
                                                                                                                    aGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGTCGACTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGcc                                  <  2:1036926/150‑1 (MQ=255)
                                                                                                                    aGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGcc                                  <  2:1030612/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)
                                                                                                                                  ttCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgtg                    <  2:1325045/150‑1 (MQ=255)
                                                                                                                                  ttCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgtg                    <  2:632672/150‑1 (MQ=255)
                                                                                                                                   tCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGc                   <  2:167329/150‑1 (MQ=255)
                                                                                                                                    cTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCa                  >  2:608013/1‑150 (MQ=255)
                                                                                                                                      gATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATc                <  2:1052260/150‑1 (MQ=255)
                                                                                                                                      gATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATc                >  2:81028/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)
                                                                                                                                          gaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTACAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGc            >  2:254646/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)
                                                                                                                                           aacaGGCGCTGCGGCGCCTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAGGctt          <  2:1274618/147‑3 (MQ=255)
                                                                                                                                           aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCa           <  2:1644524/150‑1 (MQ=255)
                                                                                                                                           aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCa           >  1:538880/1‑150 (MQ=255)
                                                                                                                                              aGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAAc        <  2:1268623/150‑1 (MQ=255)
                                                                                                                                               ggCGCGTCGTCGCCTTTTCAGTTTGCTGACAAAGTTCACTTGGTTATGCCGGATACGGCGTGAAAGCGTTATCAGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCACATCAGCCTGATGCACAATGTTTCTCAGGcttatg       <  2:1095136/150‑6 (MQ=255)
                                                                                                                                                   cTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCAc   <  2:105398/150‑1 (MQ=255)
                                                                                                                                                    ttCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACg  <  1:480334/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACG  >  NC_000913/3525288‑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: