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,1580CT97.4% 250.1 / ‑1.3 79D354N (GAC→AAC) mrdApenicillin‑binding protein 2, transpeptidase involved in peptidoglycan synthesis
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (34/41);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

CGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATA  >  NC_000913/667010‑667302
                                                                                                                                                    |                                                                                                                                                
cgatcgatCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTaa                                                                                                                                                <  2:505092/151‑1 (MQ=255)
  atcgatcCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAAc                                                                                                                                              >  2:905339/1‑151 (MQ=255)
   tcgatcCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAAc                                                                                                                                              >  2:682264/1‑150 (MQ=255)
     gatcCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGc                                                                                                                                           <  2:330693/151‑1 (MQ=255)
        ccccATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGcgtc                                                                                                                                        <  2:1045013/151‑1 (MQ=255)
        ccccATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGcgt                                                                                                                                         <  1:586893/150‑1 (MQ=255)
         cccATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTTGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGcgtcg                                                                                                                                       <  2:157555/151‑1 (MQ=255)
            atatCGTATGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTAt                                                                                                                                    >  2:427449/1‑151 (MQ=255)
             tatCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTAtt                                                                                                                                   >  1:47232/1‑151 (MQ=255)
                cGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGTCACCAGCCTGGGTTAAACAGCGTCGTATTgcg                                                                                                                                <  1:241900/151‑1 (MQ=255)
                 gTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTgcgc                                                                                                                               >  1:127519/1‑151 (MQ=255)
                 gTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTgc                                                                                                                                 <  2:860645/149‑1 (MQ=255)
                   aGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGcgcg                                                                                                                              <  1:652987/150‑1 (MQ=255)
                    ggCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTTGTATTGCGCGTg                                                                                                                            <  1:103279/151‑1 (MQ=255)
                    ggCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGCTAAACAGCGTCGTATTGcgcg                                                                                                                              >  2:912796/1‑149 (MQ=255)
                      ccACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGAt                                                                                                                          <  2:931097/151‑1 (MQ=255)
                        aCCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATc                                                                                                                         >  1:98848/1‑150 (MQ=255)
                                  agaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCt                                                                                                              >  2:66688/1‑151 (MQ=255)
                                     agGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCaa                                                                                                           <  1:667677/151‑1 (MQ=255)
                                       gTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATg                                                                                                         >  1:90313/1‑151 (MQ=255)
                                        tATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGc                                                                                                        >  2:680166/1‑151 (MQ=255)
                                           ccGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCg                                                                                                      >  2:246622/1‑150 (MQ=255)
                                           ccGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCg                                                                                                      <  1:63403/150‑1 (MQ=255)
                                               aGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAAcc                                                                                                 <  1:983244/151‑1 (MQ=255)
                                                gATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCg                                                                                                >  2:537895/1‑151 (MQ=255)
                                                 aTTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGc                                                                                               >  1:746467/1‑151 (MQ=255)
                                                 aTTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAAcc                                                                                                 >  1:149666/1‑149 (MQ=255)
                                                    cttcCCGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCCGTCACGGTAACGTTTTTCTGAACCTGGCAGTTGCCACCAGCCGGGGCTAAACAGCGTCGTATTGCACGCGATCACCCCGGCGCACAGTGCCGaa                                                                                                    >  2:413637/1‑143 (MQ=255)
                                                       ccAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACAt                                                                                          >  1:77556/1‑150 (MQ=255)
                                                           cGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGt                                                                                     >  1:670341/1‑151 (MQ=255)
                                                            gctCTTATCACATTCAGAAGACCGTGGCCCCATTTTTTCCAGTGACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGCCGTAGTGCGCATGTTCACCCCGGAGCTCAAGGCCGAAACCGCCACATAGGGtt                                                                                    >  2:541641/3‑151 (MQ=255)
                                                              tCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGtt                                                                                    <  2:468420/149‑1 (MQ=255)
                                                                    gACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTaa                                                                                 <  2:292793/146‑1 (MQ=255)
                                                                        ttCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAc                                                                         <  1:162915/150‑1 (MQ=255)
                                                                          cAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAcgcg                                                                      >  2:449125/1‑151 (MQ=255)
                                                                           aGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCgg                                                                     >  1:454468/1‑151 (MQ=255)
                                                                              cccccGTGGCCCCATTTTTTCCAGTCACGATACCGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACACAGGCTTTAACTGTAGACGCggg                                                                    >  2:885896/3‑149 (MQ=255)
                                                                                cccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACATCTTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggagg                                                                 <  2:346902/150‑1 (MQ=255)
                                                                                cccGTGGCCCCATTTTTCCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCGCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggagg                                                                 <  2:380766/150‑1 (MQ=255)
                                                                                 ccGTGGCCTTTTTTTTTCCACCCACGATAATGTTTTTCGGACCGTGGCAGTTGCCCCCAGCCAGGCTTAACCAGGGTCGTATTGCGCGGGAGCACCCGGGCGCGAAAGGCTGAATCCGCCACATAGGGATTAACTGTAGACGCGggagg                                                                 <  2:385155/149‑1 (MQ=255)
                                                                                 ccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGCCGTATTCCGCGTGATCACCCCGGCGCTCAACGCCGAAACCGCCACACAGGGTTTAACTGTAAACGCGggaggc                                                                >  2:1046683/1‑149 (MQ=255)
                                                                                    tGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaa                                                             >  2:482144/1‑150 (MQ=255)
                                                                                     ggCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCCCCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAc                                                           <  1:839483/151‑1 (MQ=255)
                                                                                       ccccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAccc                                                         >  2:1061900/1‑151 (MQ=255)
                                                                                        cccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAACcccc                                                        >  2:736226/1‑151 (MQ=255)
                                                                                           aTTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACACCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtg                                                     >  2:69099/1‑151 (MQ=255)
                                                                                            ttttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgt                                                    >  2:651641/1‑151 (MQ=255)
                                                                                             tttttCACGCCACGATATCGTTTTTCCGAACCTGGCAGTGGCCCCCAGCCTGGGTTAACCAGCGTCGTATTGCGCGTGATCACCCCGGCACTGATTGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtg                                                     <  2:520554/149‑1 (MQ=255)
                                                                                              ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCGCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtg                                                   <  1:905631/150‑1 (MQ=255)
                                                                                              ttttCCAGCCACGATAACGTTTTTCCGAACCTGGCAGTTGCCCCCAGCCGGGGTTAAACAGCGTCGTATTGCGCGTGACCTCCCGGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgt                                                    <  1:211422/149‑1 (MQ=255)
                                                                                                  ccAGTCACGATAACTTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGgc                                               <  1:514551/150‑1 (MQ=255)
                                                                                                  ccAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGgcg                                              <  1:1039983/151‑1 (MQ=255)
                                                                                                   cAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGgcgc                                             >  2:697177/1‑151 (MQ=255)
                                                                                                    aGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGcgcg                                            <  1:910458/151‑1 (MQ=255)
                                                                                                         cGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTc                                        >  1:427116/1‑150 (MQ=255)
                                                                                                          gATAACGTTTTTCCGAACCTGGCCGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCAc                                      >  2:74987/1‑151 (MQ=255)
                                                                                                          gATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCAc                                      >  1:772853/1‑151 (MQ=255)
                                                                                                              aCGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGTTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGc                                  <  1:83063/151‑1 (MQ=255)
                                                                                                              aCGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGc                                  >  1:833031/1‑151 (MQ=255)
                                                                                                               cGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCg                                 <  1:331321/151‑1 (MQ=255)
                                                                                                                 tttttCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCgg                                >  2:885801/1‑150 (MQ=255)
                                                                                                                     tCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTa                            <  1:325685/150‑1 (MQ=255)
                                                                                                                     tCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTAt                           >  2:616370/1‑151 (MQ=255)
                                                                                                                      ccGAACCGGGCGGTTGCCCCCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCGGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTAtt                          <  1:402257/151‑1 (MQ=255)
                                                                                                                       cGACCCTGGCAGTTGCCCCCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc                         <  1:71654/151‑1 (MQ=255)
                                                                                                                        gAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTAAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgtgt                             <  2:76515/146‑1 (MQ=255)
                                                                                                                         aaCCGGGAAGTGGCCACCAGCGTGGGTTAAACAGCGCCTTATTGCGCGCGTTCACCCGGGCGCTCAATCCCCAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGTGGTAttg                         <  2:1004153/149‑2 (MQ=255)
                                                                                                                          aCCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGa                      <  1:190626/151‑1 (MQ=255)
                                                                                                                              ggCGGTTGCCCCCAGCCTGGTTTAAACGGCGTCGTATTGCGCGTGATCTCCCCGGCGCTAAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGt                  <  1:1010230/151‑1 (MQ=255)
                                                                                                                              ggCAGTTGCCCCCACCCTGGGTAAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCg                   <  1:825416/150‑1 (MQ=255)
                                                                                                                              ggCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATGGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGa                      <  2:122696/147‑1 (MQ=255)
                                                                                                                              ggCAGTCGCCACCTGCCTGGGTAAACCGGCGTCGGATTGCCCGGGATCCCCCCGCCGCTCAATTCCTAACCCGCCACATAGGGTTTAACTGTAGACGCGGGTGGATAAACCCCCTGTGTGGCGCGGTTCACTAGCGGTGTATTCGGAt                     <  2:602246/148‑1 (MQ=255)
                                                                                                                                cAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGGGGGATCCCCCCGGGGCTCAAGGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTc                <  2:492123/151‑1 (MQ=255)
                                                                                                                                   ttGCCACCAGCCTGGGTTAAACGGCGTCGTATCGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTca               <  2:205774/149‑1 (MQ=255)
                                                                                                                                     gCCACCAGCCTGGTTTAACCAGCGTCGTATGGCGCGTGATCCCCCCGGCGCTCAATGCCGAAGCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTcaacaa           <  2:469914/151‑1 (MQ=255)
                                                                                                                                     gCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGGATTCGGATCGTTcaaca            >  2:985937/1‑150 (MQ=255)
                                                                                                                                       caccaGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCATTGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAgg         <  1:112544/151‑1 (MQ=255)
                                                                                                                                               tGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaata  >  1:258583/1‑150 (MQ=255)
                                                                                                                                                gggTTAAACAGCGTCGTATTGCGCGTGACCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaata  <  1:1067905/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATA  >  NC_000913/667010‑667302

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