Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 667,158 C→T D354N (GAC→AAC)  mrdA ← penicillin‑binding protein 2, transpeptidase involved in peptidoglycan synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913667,1580CT98.8% 294.2 / ‑3.9 81D354N (GAC→AAC) mrdApenicillin‑binding protein 2, transpeptidase involved in peptidoglycan synthesis
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (46/33);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

GATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAA  >  NC_000913/667011‑667303
                                                                                                                                                   |                                                                                                                                                 
gatcggtcCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTaa                                                                                                                                                 >  2:940444/1‑150 (MQ=255)
  tcgatcCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAAc                                                                                                                                               >  2:953617/1‑150 (MQ=255)
   cgatcCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAAc                                                                                                                                               >  2:501151/1‑149 (MQ=255)
    gatcCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGc                                                                                                                                            >  2:534552/1‑151 (MQ=255)
    gatcCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGc                                                                                                                                            <  2:28437/151‑1 (MQ=255)
     atcCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGcg                                                                                                                                           >  2:232403/1‑151 (MQ=255)
       ccccATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGcgtc                                                                                                                                         <  1:120646/151‑1 (MQ=255)
       ccccATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGcgtc                                                                                                                                         <  1:760084/151‑1 (MQ=255)
        cccATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGcgtcg                                                                                                                                        <  1:258525/151‑1 (MQ=255)
        cccATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGcgtcg                                                                                                                                        <  1:31340/151‑1 (MQ=255)
          cATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTa                                                                                                                                      <  1:473416/151‑1 (MQ=255)
           atatCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTAt                                                                                                                                     <  1:931064/151‑1 (MQ=255)
                 tAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGcgcg                                                                                                                               >  1:317240/1‑151 (MQ=255)
                   gaccAGCTAATAGAAGAAGGTAACCGGAGGTTCTCCCAGCGATCGTGTGCCAGTCAGACGCGCGTGGCCCCATTTTGTCCAGTCACGATAACGTTTTTCCGAACATGGCAGTTGCCACCAGCCTGGGTGAAGCAGCGTCGTATTGCGCGta                                                                                                                             <  1:471931/149‑2 (MQ=255)
                       aCCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATc                                                                                                                          >  2:602478/1‑150 (MQ=255)
                        ccTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCAc                                                                                                                        <  1:68416/151‑1 (MQ=255)
                          tGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCAccc                                                                                                                      >  2:78547/1‑151 (MQ=255)
                            aTAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCg                                                                                                                    >  2:1189224/1‑151 (MQ=255)
                             tAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCg                                                                                                                    >  1:254960/1‑150 (MQ=255)
                              agaagaAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGgc                                                                                                                  >  2:477262/1‑151 (MQ=255)
                              agaagaAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGgc                                                                                                                  <  2:1223744/151‑1 (MQ=255)
                               gaagaaGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGgcg                                                                                                                 >  1:859819/1‑151 (MQ=255)
                                 agaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGgcgc                                                                                                                >  1:563275/1‑150 (MQ=255)
                                   aagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTc                                                                                                              >  1:647071/1‑150 (MQ=255)
                                    agGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCaa                                                                                                            <  1:290101/151‑1 (MQ=255)
                                      gTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATg                                                                                                          >  2:551467/1‑151 (MQ=255)
                                       tATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGc                                                                                                         >  1:213090/1‑151 (MQ=255)
                                       tATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGc                                                                                                         >  1:250599/1‑151 (MQ=255)
                                            gCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaa                                                                                                     >  2:22096/1‑150 (MQ=255)
                                              aGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaaa                                                                                                    >  2:598477/1‑149 (MQ=255)
                                                     tcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGACCCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCca                                                                                              <  1:622419/148‑1 (MQ=255)
                                                           gATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAggg                                                                                       >  2:54393/1‑149 (MQ=255)
                                                           gATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGtt                                                                                     >  1:502476/1‑151 (MQ=255)
                                                           gATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGtt                                                                                     >  1:703488/1‑151 (MQ=255)
                                                           gATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGtt                                                                                     >  1:839854/1‑151 (MQ=255)
                                                                tgtgACATTCAGACGCCCGTGGCCCTTTTTTTTCCAGTCACGATAACGTTTTTCCGACCCTGGCAGTTCCCCCCAGCCTGGGTTAAACAGCGTCCTTTTGCGCGTAACCACCCCGGGGCTAAAGGCCGAAACCGCCACATAGGGtt                                                                                     <  1:1168625/146‑1 (MQ=255)
                                                                tgtgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACt                                                                                >  2:783372/1‑151 (MQ=255)
                                                                  tgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGt                                                                              >  2:115485/1‑151 (MQ=255)
                                                                    aCATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAg                                                                            <  2:311110/151‑1 (MQ=255)
                                                                     cATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGt                                                                              >  2:661197/1‑148 (MQ=255)
                                                                      aTTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAc                                                                          <  2:320717/151‑1 (MQ=255)
                                                                      aTTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAc                                                                          >  2:373973/1‑151 (MQ=255)
                                                                        tCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAcg                                                                         <  2:1057683/150‑1 (MQ=255)
                                                                                  gTGGCCCCATTTTTTCCAGTCACGTTAACGTTTTCCCGAACCTGGCAGTTGCCCCCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCAGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaa                                                              <  1:185710/151‑1 (MQ=255)
                                                                                    ggCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAc                                                            <  2:198888/151‑1 (MQ=255)
                                                                                       cccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAcccc                                                         <  2:995605/151‑1 (MQ=255)
                                                                                         cATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAccccc                                                        >  1:798190/1‑150 (MQ=255)
                                                                                         cATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCt                                                       >  2:752651/1‑151 (MQ=255)
                                                                                           ttttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgt                                                     >  2:403028/1‑151 (MQ=255)
                                                                                           tttttcCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAACCAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAAGGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtg                                                      <  2:548093/150‑1 (MQ=255)
                                                                                            tttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtg                                                    >  1:545726/1‑151 (MQ=255)
                                                                                             ttttCAATTCACGATAACTTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAACCAGCGTCGTATTGCCCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtg                                                    <  2:388216/150‑1 (MQ=255)
                                                                                              tttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACGGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtgt                                                   <  2:482548/150‑1 (MQ=255)
                                                                                              tttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtgtg                                                  >  2:51543/1‑151 (MQ=255)
                                                                                                tCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtgtg                                                  >  1:169576/1‑149 (MQ=255)
                                                                                                  cAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGgcgc                                              <  2:647241/151‑1 (MQ=255)
                                                                                                   aGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGcgcg                                             >  1:344448/1‑151 (MQ=255)
                                                                                                    gTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCgg                                            <  2:956012/151‑1 (MQ=255)
                                                                                                    gTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGCTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCgg                                            >  1:379893/1‑151 (MQ=255)
                                                                                                      cACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGt                                           <  2:1196433/150‑1 (MQ=255)
                                                                                                        cGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTc                                         >  2:359440/1‑150 (MQ=255)
                                                                                                        cGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTc                                         >  1:597750/1‑150 (MQ=255)
                                                                                                         gATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCa                                        <  2:933766/150‑1 (MQ=255)
                                                                                                               gTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCgg                                 >  1:545194/1‑151 (MQ=255)
                                                                                                               gTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCgg                                 >  1:1007546/1‑151 (MQ=255)
                                                                                                                     ccGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgcgt                              >  1:153658/1‑148 (MQ=255)
                                                                                                                      cGACCCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTAt                            <  1:78535/149‑1 (MQ=255)
                                                                                                                        aaCCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCgg                        <  2:1217610/151‑1 (MQ=255)
                                                                                                                        aaCCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCg                         <  1:1144305/150‑1 (MQ=255)
                                                                                                                          ccTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGAt                      >  2:736343/1‑151 (MQ=255)
                                                                                                                          ccTGGCAGTGGCCCCCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCGGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCAGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGAt                      <  1:677763/151‑1 (MQ=255)
                                                                                                                            tGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGGTCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTCTTCGGATCg                    >  1:388548/1‑151 (MQ=255)
                                                                                                                            tGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCg                    >  1:79272/1‑151 (MQ=255)
                                                                                                                              gCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGtt                  <  2:403370/151‑1 (MQ=255)
                                                                                                                                     ccaccaGCCTGGTTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTcaacaa            <  2:954727/150‑1 (MQ=255)
                                                                                                                                     ccaccaGCCGGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAg           <  1:386385/151‑1 (MQ=255)
                                                                                                                                       accaGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGc         >  2:14499/1‑151 (MQ=255)
                                                                                                                                             cTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGGGTATTCGGATCGTTCAACAAGGCGGaat    >  1:746215/1‑150 (MQ=255)
                                                                                                                                              tGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaata   >  1:848655/1‑150 (MQ=255)
                                                                                                                                               gggTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaataa  >  2:402410/1‑150 (MQ=255)
                                                                                                                                               gggTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaata   <  2:981779/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAA  >  NC_000913/667011‑667303

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