Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,647,150 G→A A59V (GCT→GTT)  pbpC ← penicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,647,1500GA98.5% 253.6 / ‑4.9 66A59V (GCT→GTT) pbpCpenicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (34/31);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GAGATCTTGCCATGCTGCGCGCGCCACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGTTATCCAG  >  NC_000913/2647005‑2647296
                                                                                                                                                 |                                                                                                                                                  
gagaTCTTGCCATGCTGCGCGCGCCACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCg                                                                                                                                               >  1:666512/1‑151 (MQ=255)
       tGCCATGCTGCGCGCGCCACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCg                                                                                                                                            <  1:599849/147‑1 (MQ=255)
       tGCCATGCTGCGCGCGCCACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGaa                                                                                                                                          <  1:762463/149‑1 (MQ=255)
               tgcGCGCGCCACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAgc                                                                                                                                >  1:109332/1‑151 (MQ=255)
                gcgcgcgcCACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAgcg                                                                                                                               <  1:138076/151‑1 (MQ=255)
                gcgcgcgcCACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAgc                                                                                                                                <  1:1250028/150‑1 (MQ=255)
                 cgcgcgcCACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAgcgg                                                                                                                              <  2:269040/151‑1 (MQ=255)
                    gcgcCACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAATGGAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATAGGCGAAGCGCCAGAGCGGCGt                                                                                                                           >  2:859238/1‑151 (MQ=255)
                       ccACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTAct                                                                                                                        <  1:776782/151‑2 (MQ=255)
                            gagaATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCt                                                                                                                   >  1:971856/1‑151 (MQ=255)
                             agaATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCt                                                                                                                   >  1:1272517/1‑150 (MQ=255)
                             agaATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTg                                                                                                                  <  2:45006/151‑1 (MQ=255)
                                aTGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTgcgc                                                                                                               >  2:1046996/1‑151 (MQ=255)
                                   gATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCAc                                                                                                            <  2:1073588/151‑1 (MQ=255)
                                    aTTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACg                                                                                                           >  2:861481/1‑151 (MQ=255)
                                            ccGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTca                                                                                                   <  1:298308/151‑2 (MQ=255)
                                                aTGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCg                                                                                               >  2:610360/1‑151 (MQ=255)
                                                aTGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCg                                                                                               >  2:1110307/1‑151 (MQ=255)
                                                 tGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCgg                                                                                              >  2:138349/1‑151 (MQ=255)
                                                   cTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCggg                                                                                             >  1:1562/1‑150 (MQ=255)
                                                     tCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGAt                                                                                           <  1:1180446/150‑1 (MQ=255)
                                                       cAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTg                                                                                         >  2:261420/1‑150 (MQ=255)
                                                        aGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGAc                                                                                       >  2:216272/1‑151 (MQ=255)
                                                           aCCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTc                                                                                    <  2:423804/151‑1 (MQ=255)
                                                            ccAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTc                                                                                    >  2:349896/1‑150 (MQ=255)
                                                               gCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGc                                                                                >  1:905105/1‑151 (MQ=255)
                                                                  aTCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGc                                                                             >  2:868554/1‑151 (MQ=255)
                                                                      tCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGc                                                                             <  1:1059161/147‑1 (MQ=255)
                                                                      tCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGc                                                                             <  1:381593/147‑1 (MQ=255)
                                                                      tCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCgg                                                                           <  2:619620/149‑1 (MQ=255)
                                                                      tCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCgg                                                                           <  2:617755/149‑1 (MQ=255)
                                                                            ttGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGcc                                                                   >  1:930295/1‑151 (MQ=255)
                                                                             tGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGcc                                                                   >  2:888191/1‑150 (MQ=255)
                                                                                 cAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGt                                                              >  1:1236222/1‑151 (MQ=255)
                                                                                  aGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGtt                                                             <  2:840420/151‑1 (MQ=255)
                                                                                  aGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGtt                                                             >  2:216622/1‑151 (MQ=255)
                                                                                         cAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTg                                                       <  1:479696/150‑1 (MQ=255)
                                                                                          caGGTAACGGGGAGAAACATCTTCGAATGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTCCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGc                                                      >  1:1223208/2‑150 (MQ=255)
                                                                                            ggTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTcc                                                   <  2:443297/151‑1 (MQ=255)
                                                                                               aaCGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTcccc                                                 >  1:281870/1‑150 (MQ=255)
                                                                                               aaCGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCAGCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCa                                                <  1:491109/151‑1 (MQ=255)
                                                                                                aCGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCAc                                               <  1:120868/151‑1 (MQ=255)
                                                                                                 cGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACg                                              >  2:815227/1‑151 (MQ=255)
                                                                                                       gaAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCa                                         >  2:175873/1‑150 (MQ=255)
                                                                                                           cATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTACGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGaa                                     >  1:446724/1‑150 (MQ=255)
                                                                                                                tCGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGaaaa                                   <  2:963920/147‑1 (MQ=255)
                                                                                                                 cGATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAgg                                >  2:995365/1‑149 (MQ=255)
                                                                                                                  gATTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGgaga                             <  1:163663/151‑1 (MQ=255)
                                                                                                                   aTTGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGAGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGgaga                             >  1:948714/1‑150 (MQ=255)
                                                                                                                    ttGTTACCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGaa                            >  2:63798/1‑150 (MQ=255)
                                                                                                                        tACCAGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGGCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAggg                        <  2:730327/150‑1 (MQ=255)
                                                                                                                         aCCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCACGAAAAAGGAGAAAGGgc                       >  1:1056360/1‑150 (MQ=255)
                                                                                                                         aCCGGATAACGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGgcg                      <  2:699297/151‑1 (MQ=255)
                                                                                                                                aaCGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAATAAGGAGAAAGGGCGCGGCt                 >  1:399857/1‑149 (MQ=255)
                                                                                                                                aaCGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGc               <  2:678182/151‑1 (MQ=255)
                                                                                                                                 aCGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGcc              >  1:1069356/1‑151 (MQ=255)
                                                                                                                                  cGCCAGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCa             <  2:53390/151‑1 (MQ=255)
                                                                                                                                      aGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGt         >  2:1245708/1‑151 (MQ=255)
                                                                                                                                      aGATGCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGt         <  1:1110269/151‑1 (MQ=255)
                                                                                                                                          gCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGTTAt      <  1:137075/150‑1 (MQ=255)
                                                                                                                                          gCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGTTAt      <  2:90758/150‑1 (MQ=255)
                                                                                                                                          gCCGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGTTATc     <  1:868514/151‑1 (MQ=255)
                                                                                                                                           ccGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGTTATcc    >  2:1127425/1‑151 (MQ=255)
                                                                                                                                            cGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGTTATcc    >  2:276849/1‑150 (MQ=255)
                                                                                                                                            cGTCAACATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGTTATCca   <  1:351878/151‑1 (MQ=255)
                                                                                                                                              tCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGTTATCcag  <  2:466764/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GAGATCTTGCCATGCTGCGCGCGCCACCGAGAATGGATTCACCCCCGGATGCTTCCAGAACCAGCGATCTTCATAATTGATCAGCGCTTCAAGGTAACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGGTAGAGGCCAGAGTTTATCTGCTCCCCACGCTGCAAGAAAAAGGAGAAAGGGCGCGGCTGCCAACGTTATCCAG  >  NC_000913/2647005‑2647296

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