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,1580CT96.2% 361.1 / ‑5.1 104D354N (GAC→AAC) mrdApenicillin‑binding protein 2, transpeptidase involved in peptidoglycan synthesis
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (52/48);  minor base A (0/4);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATC  >  NC_000913/667010‑667305
                                                                                                                                                    |                                                                                                                                                   
cgatcgatCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTaa                                                                                                                                                   <  2:1368031/151‑1 (MQ=255)
 gatcgatcCCCATATCGTAGGCCACCTGATAGAAGAAGGGATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTaa                                                                                                                                                   >  2:17193/1‑150 (MQ=255)
     gatcCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGc                                                                                                                                              >  2:1380978/1‑151 (MQ=255)
        ccccATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGcgtc                                                                                                                                           <  1:608443/151‑1 (MQ=255)
         cccATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGcgtcg                                                                                                                                          <  2:1286183/151‑1 (MQ=255)
              atCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTg                                                                                                                                     <  1:1217888/151‑1 (MQ=255)
               tCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTg                                                                                                                                     >  2:1236524/1‑150 (MQ=255)
                cGTAGGCCACCTGATAGAAGAAGGTTTCCGCAGATTCTTCCCGCGATCTTGTGACATTCATACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTg                                                                                                                                     <  2:1335982/149‑1 (MQ=255)
                  tAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGcgcg                                                                                                                                 >  2:374895/1‑151 (MQ=255)
                   aGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTTTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGt                                                                                                                                <  2:54007/151‑1 (MQ=255)
                   aGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGcgcg                                                                                                                                 >  2:311978/1‑150 (MQ=255)
                   aGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGcgcg                                                                                                                                 <  2:391958/150‑1 (MQ=255)
                        aCCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATc                                                                                                                            >  1:836880/1‑150 (MQ=255)
                        aCCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCa                                                                                                                           <  1:407360/151‑1 (MQ=255)
                          cTGATAGAAGAAGGTATCCGCCGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCAcc                                                                                                                         >  2:814541/1‑151 (MQ=255)
                              tAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCgg                                                                                                                     >  2:328188/1‑151 (MQ=255)
                               cgaagccggTATCCGGAGCTGTTTTCAGCTATCGTGTGACATTCAGACGCCCGTGGCCTCATTTTTTCCAGTAACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCATCGGCTAAACAGCGCCGTATTGCGTGTGATCcccc                                                                                                                        >  2:791373/8‑147 (MQ=255)
                                 aagaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGgcgc                                                                                                                  >  2:570536/1‑151 (MQ=255)
                                   gaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTc                                                                                                                >  1:1298064/1‑151 (MQ=255)
                                    aagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCa                                                                                                               <  1:1033642/151‑1 (MQ=255)
                                       gTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAAt                                                                                                             <  1:105201/150‑1 (MQ=255)
                                        tATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGc                                                                                                           >  1:583562/1‑151 (MQ=255)
                                         aTCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGcc                                                                                                          <  1:69565/151‑1 (MQ=255)
                                         aTCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGcc                                                                                                          >  1:1205503/1‑151 (MQ=255)
                                         aTCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGc                                                                                                           <  1:1381390/150‑1 (MQ=255)
                                         aTCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGc                                                                                                           <  2:1048848/150‑1 (MQ=255)
                                         aTCCGCAGATTATTCCAGCGATCTTGTGACGTTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGAGAACGTTTTTCCGAACCTGACAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATGGCTCGTGATCACCCCGGCGCTCAATg                                                                                                            >  2:736607/1‑149 (MQ=255)
                                            cGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGa                                                                                                        >  1:392090/1‑150 (MQ=255)
                                             gCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaaa                                                                                                      <  2:169229/151‑1 (MQ=255)
                                                   tcttcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGcc                                                                                                 >  2:303466/1‑150 (MQ=255)
                                                         aGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAgg                                                                                          <  1:692539/151‑1 (MQ=255)
                                                         aGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAgg                                                                                          >  2:349361/1‑151 (MQ=255)
                                                            gATTTTCTGACATTCAGACGCTCGAGGCCCGATTTTTTCCGAACAGAAAAACGTGTATCCACCCGCGGCATGTGCCCCCAGCCCGGGTTAAACAGCCTATTATTGCCGGTGATCGCCCGGGCCTTCAATGCCGAAACCGCCACATAgtgt                                                                                        <  2:855663/150‑1 (MQ=255)
                                                            gATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGtt                                                                                       >  1:250430/1‑151 (MQ=255)
                                                               cTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTa                                                                                     >  1:802289/1‑150 (MQ=255)
                                                                ttGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAAc                                                                                   >  2:692764/1‑151 (MQ=255)
                                                                 tgttACATTCAGACGCCCGTGGCCCCATTTTTCCCAGTCGCGATAACGTTTTTCCGAACCTGGCAGTTGCCCCCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTaa                                                                                    <  2:629653/149‑1 (MQ=255)
                                                                 tgtgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGt                                                                                        <  1:673824/145‑1 (MQ=255)
                                                                 tgtgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTa                                                                                     <  1:967746/148‑1 (MQ=255)
                                                                 tgtgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACt                                                                                  >  2:1228403/1‑151 (MQ=255)
                                                                  gttccATTCAGACGCCCGTGCCCCCTTTTTTCCCATTCACGATAACGTTTTCCCGACCCGGGCAGTGCCCCCAAGCCTGGGTAAAACAGCGTCGTATTGCGGTTTAGCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTaa                                                                                    <  2:644591/144‑1 (MQ=255)
                                                                    gACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGt                                                                                >  1:871494/1‑150 (MQ=255)
                                                                    gACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTa                                                                               <  1:1342102/151‑1 (MQ=255)
                                                                       aTTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAc                                                                            >  1:1067248/1‑151 (MQ=255)
                                                                           aGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCgg                                                                        >  2:708002/1‑151 (MQ=255)
                                                                            gACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCggg                                                                       >  1:333815/1‑151 (MQ=255)
                                                                               gCCCGTGTCCCCCTTTTTCCCTGTCTCTATAACGTTTTTCCGACCCGGTCAGTTGCCACCGCCCTGGGTAAAACGGCGTCGTATTGCGCGTAATCCCCCCTCAGCGCAATGCCGAAACCGCCACATAGGGTTTAACGGTAGACTCGggagg                                                                    <  2:1392375/151‑1 (MQ=255)
                                                                                cccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggagg                                                                    >  2:66133/1‑150 (MQ=255)
                                                                                 ccggagCCCCGTTTTTTCAAATAACGATACCATTTTTCCTACCGGCGGAGTGGCCACCATCCTGAGTTAAACAACTTCGTATTGGGCGTAATCGCCGCGGCCCTCAATACCGAAACCGCCACGTAGGGTTTAACTGTAGACGCGggagt                                                                    <  2:1202564/144‑2 (MQ=255)
                                                                                 ccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggagga                                                                   >  2:456938/1‑150 (MQ=255)
                                                                                 ccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggagga                                                                   >  2:706817/1‑150 (MQ=255)
                                                                                  cGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATa                                                                 <  1:690968/151‑1 (MQ=255)
                                                                                    tGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTAAATCCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaa                                                                <  1:1104282/150‑1 (MQ=255)
                                                                                     ggCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGACCCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaaa                                                               <  1:916351/150‑1 (MQ=255)
                                                                                     ggCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaaa                                                               >  1:373804/1‑150 (MQ=255)
                                                                                       ccccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGGAGTTGCCCCCAGCCTGGGTTAAACAGCGTCGTATGGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAcc                                                             <  1:721485/150‑1 (MQ=255)
                                                                                           aTTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtg                                                        >  2:443995/1‑151 (MQ=255)
                                                                                            ttttttCCAGTCACGATACCGTTTTTCCGAACCTGGCAGTTGCCCCCAGCCTGGTTTAAACAGCGTCGTATTGCGCTTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCt                                                         <  2:851207/149‑1 (MQ=255)
                                                                                             tttttCCAGTCGCGATAACGTTTTTCCGAACCTGGCAGTTGCCCCCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgt                                                       <  2:1014878/150‑1 (MQ=255)
                                                                                              ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtg                                                      >  2:200653/1‑150 (MQ=255)
                                                                                                    aGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGcgcg                                               >  2:888401/1‑151 (MQ=255)
                                                                                                      tCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGt                                             >  2:516951/1‑151 (MQ=255)
                                                                                                      tCACGATAACGCTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGACATCACCCCGACGCTCAATGCCGAAACCGCCACATAGGGCTTAACTGTAGACGCGGGAGGATACACCCCCTGTGTGGCGCGGt                                             >  2:751603/1‑151 (MQ=255)
                                                                                                         cGATACCGTTTTTCCGACCCGGGCAGTTGCCACCACCCTGGGTTAACCAGCGTCGTATTGCGTGTGATCCCCCCGGCGCTCATTGCCGAACCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTc                                           <  2:275780/150‑1 (MQ=255)
                                                                                                            tAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCAcc                                        >  1:1377407/1‑150 (MQ=255)
                                                                                                            tAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCAcc                                        >  1:858853/1‑150 (MQ=255)
                                                                                                             aaCGTTTTCCCGACCCTGGCATTGTCCACCAGCGTGGGTTAAACAGCGTCGTATTGCGGGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAg                                      <  2:238611/151‑1 (MQ=255)
                                                                                                               cGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCg                                    >  1:859696/1‑151 (MQ=255)
                                                                                                               cGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGGGTGGCGCGGTTCACCAGCg                                    >  2:103563/1‑151 (MQ=255)
                                                                                                                gTTTTTCGGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGCGATAACCCCGGCGATCAATGCCGCAACCACCACATAGGGTTTAACTGTAGACGCAGGAGGATAAACCCCCTGTGTGCCGTGGTTCACCAGCgg                                   >  2:762831/1‑151 (MQ=255)
                                                                                                                gTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCgg                                   >  1:1336855/1‑151 (MQ=255)
                                                                                                                  ttttCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgt                                  <  1:457098/150‑1 (MQ=255)
                                                                                                                  ttttCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgtg                                 >  2:80330/1‑151 (MQ=255)
                                                                                                                  ttttCCGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGACGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACCCGGGAGGATAAACCCCCTGTGTCGCGCGGTTCACCAGCGgtg                                 >  2:1144744/1‑151 (MQ=255)
                                                                                                                  ttttCCGAACCTGGCAGTTCCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAACCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgt                                  <  2:1154555/150‑1 (MQ=255)
                                                                                                                   tttCCGACCCTGGCAGTTGCCCCCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTATAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgtg                                 <  1:547871/150‑1 (MQ=255)
                                                                                                                      ccGAACCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAACCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTa                               <  2:864644/149‑1 (MQ=255)
                                                                                                                         aaCCTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCgg                          <  1:116906/151‑1 (MQ=255)
                                                                                                                          aCCTGGCAGTTGCCACCAGCCTGGGTTAAACAACGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACACAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGGGGCGCCGCTCACCATCGGTGTATTc                            >  1:900839/1‑148 (MQ=255)
                                                                                                                           ccgggCAGTTGCCACCAGCCTGGGTAAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGAt                        <  2:863764/148‑1 (MQ=255)
                                                                                                                           ccTGGCAGTTGCCACCAGCCTGGGTTAAACAGGGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGa                         <  1:621082/150‑1 (MQ=255)
                                                                                                                           ccTGGCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGa                         <  1:1278214/150‑1 (MQ=255)
                                                                                                                           ccTGGCAGTGGCCCCCAGCCGGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGAt                        <  2:424678/151‑1 (MQ=255)
                                                                                                                               gCAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGCCGCTCAATGCCAAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                       <  2:1166147/148‑1 (MQ=255)
                                                                                                                                cAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTc                   >  1:492219/1‑151 (MQ=255)
                                                                                                                                cAGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTc                   >  1:969188/1‑151 (MQ=255)
                                                                                                                                 aGTTGCCACCAGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTc                   >  1:1349799/1‑150 (MQ=255)
                                                                                                                                    tGCCACCCGCCTGGGTTAAACAGCGTCGTATGGCGGGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTcaac                <  2:1041955/150‑1 (MQ=255)
                                                                                                                                     gCCCCCAGCCTGGTTTAAACAGCGTCGGATTGCGCGTGATCCCCCGGGGGCTAAATGCCAAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTcaaca               <  1:1220717/150‑1 (MQ=255)
                                                                                                                                      cccccaGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCACAATGCCGAATCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAATCCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTcaac                <  1:17195/145‑1 (MQ=255)
                                                                                                                                       ccccaGCCTGGGTTAAACGGCGTGGTATGGCGGGTGTGCGCCCCGTGGTTGAGTGCCGAAACCGCCACATAGGGTGTAGCTGGAGACGCGGGTGGATAAACTCCCTGTGTGGCGGCGTTCACCAGCGATGTATTCGGATCGtt                    <  2:782676/141‑1 (MQ=255)
                                                                                                                                       caccaGCCTGGTTTAAACAGCGTCGTATTGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAgg            <  1:419386/151‑1 (MQ=255)
                                                                                                                                       caccaGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGGGTATTCGGATCGTTCAACAAgg            >  2:311573/1‑151 (MQ=255)
                                                                                                                                        accaGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAgg            >  2:1425949/1‑150 (MQ=255)
                                                                                                                                         ccaGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCg          >  1:1197947/1‑151 (MQ=255)
                                                                                                                                         ccaGCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCg          >  2:736052/1‑151 (MQ=255)
                                                                                                                                            gCCTGGGTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGa        <  1:595816/150‑1 (MQ=255)
                                                                                                                                                gggTTAAACAGGGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaataat   <  2:275523/151‑1 (MQ=255)
                                                                                                                                                gggTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaataat   <  1:1319788/151‑1 (MQ=255)
                                                                                                                                                 ggTTAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATc  <  1:618675/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATC  >  NC_000913/667010‑667305

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: ‑