Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,223,035 A→G intergenic (‑136/+29) pbpG ← / ← yohC penicillin‑binding protein 7/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,223,0350AG99.7% 1118.3 / ‑5.1 325intergenic (‑136/+29)pbpG/yohCpenicillin‑binding protein 7/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (146/177);  total (147/178)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

Only 100 of 330 total aligned reads displayed.

GATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGT  >  NC_000913/2222945‑2223124
                                                                                          |                                                                                         
gATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCg                                                                                         >  2:362576/1‑93 (MQ=255)
  tttACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGcc                                                                                       <  1:856611/93‑1 (MQ=255)
  tttACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGcc                                                                                       <  2:6685554/93‑1 (MQ=255)
   ttACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCa                                                                                      >  2:8166781/1‑93 (MQ=255)
    tACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAg                                                                                     <  2:2854528/93‑1 (MQ=255)
    tACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAg                                                                                     >  2:5003462/1‑93 (MQ=255)
        gAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAAt                                                                                 <  2:3270782/93‑1 (MQ=255)
          ggCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGc                                                                               >  2:3425491/1‑93 (MQ=255)
             ttgcgcCGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCATCAGCAATGCAAc                                                                            >  2:5969349/3‑93 (MQ=255)
              aGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACg                                                                           >  1:4827523/1‑93 (MQ=255)
                cGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTc                                                                         >  2:5851663/1‑93 (MQ=255)
                 gACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCa                                                                        <  1:730795/93‑1 (MQ=255)
                  aCGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGGATTGCTGACATTTGCGCGAGCAATGCAACGTcca                                                                       >  1:7859608/1‑91 (MQ=255)
                  aCGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCaa                                                                       <  1:4008883/93‑1 (MQ=255)
                    gCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAg                                                                     >  1:7997598/1‑93 (MQ=255)
                     ccATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGa                                                                    >  1:8333831/1‑93 (MQ=255)
                      catcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGaa                                                                   <  2:2467642/93‑1 (MQ=255)
                       atcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAAc                                                                  >  2:6327326/1‑93 (MQ=255)
                       atcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAAc                                                                  <  2:766753/93‑1 (MQ=255)
                        tcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACt                                                                 >  1:6296452/1‑93 (MQ=255)
                        tcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACt                                                                 <  2:5758676/93‑1 (MQ=255)
                        tcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACt                                                                 >  1:4181022/1‑93 (MQ=255)
                        tcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACt                                                                 <  2:4722062/93‑1 (MQ=255)
                          atAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACttt                                                               <  1:1020320/93‑1 (MQ=255)
                          atAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACttt                                                               <  2:6447528/93‑1 (MQ=255)
                          atAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACttt                                                               >  2:3899303/1‑93 (MQ=255)
                           tAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTc                                                              <  2:4518503/93‑1 (MQ=255)
                           tAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTc                                                              >  1:6409753/1‑93 (MQ=255)
                            aaCGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCt                                                             <  2:1857763/93‑1 (MQ=255)
                            aaCGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCt                                                             <  2:3031821/93‑1 (MQ=255)
                             aCGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTa                                                            <  1:6896929/93‑1 (MQ=255)
                             aCGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTa                                                            >  1:7165657/1‑93 (MQ=255)
                               gAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGa                                                          >  1:5713932/1‑93 (MQ=255)
                                aGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGaa                                                         <  1:1361915/93‑1 (MQ=255)
                                  cAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACa                                                       <  1:3310880/93‑1 (MQ=255)
                                   aaaaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACaa                                                      >  2:4284261/1‑93 (MQ=255)
                                    aaaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACaaa                                                     >  2:664753/1‑93 (MQ=255)
                                     aaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAAc                                                    <  1:6801437/93‑1 (MQ=255)
                                     aaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAAc                                                    <  2:859176/93‑1 (MQ=255)
                                     aaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAAc                                                    >  2:514671/1‑93 (MQ=255)
                                     aaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAAc                                                    >  1:6969219/1‑93 (MQ=255)
                                       aGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGa                                                  >  1:149572/1‑93 (MQ=255)
                                       aGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGa                                                  >  2:2638238/1‑93 (MQ=255)
                                        gTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGAt                                                 <  2:4530461/93‑1 (MQ=255)
                                           cGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAAc                                              <  1:3378400/93‑1 (MQ=255)
                                           cGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAAc                                              <  1:3705880/93‑1 (MQ=255)
                                           cGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAAc                                              <  1:4931300/93‑1 (MQ=255)
                                           cGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAAc                                              <  1:4676802/93‑1 (MQ=255)
                                            gATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAAcc                                             >  2:708533/1‑93 (MQ=255)
                                               cGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTa                                          <  2:5630688/93‑1 (MQ=255)
                                               cGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTa                                          <  2:8516702/93‑1 (MQ=255)
                                               cGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTa                                          >  1:8685399/1‑93 (MQ=255)
                                                gTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAg                                         >  2:1483422/1‑93 (MQ=255)
                                                gTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAg                                         >  1:5536547/1‑93 (MQ=255)
                                                 tCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGc                                        >  1:1420916/1‑93 (MQ=255)
                                                  cAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGcc                                       >  1:1515422/1‑93 (MQ=255)
                                                   aaaGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGccc                                      >  2:5350724/1‑93 (MQ=255)
                                                   aaaGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGccc                                      >  2:7368396/1‑93 (MQ=255)
                                                    aaGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGcccc                                     >  2:721972/1‑93 (MQ=255)
                                                     aGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCa                                    >  2:6075647/1‑93 (MQ=255)
                                                      ggAGAATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCAt                                   <  2:4595356/93‑1 (MQ=255)
                                                        agaATCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAg                                 <  1:3004084/93‑1 (MQ=255)
                                                          aaTCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGta                               <  1:8751278/93‑1 (MQ=255)
                                                           aTCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGtat                              <  2:8647504/93‑1 (MQ=255)
                                                           aTCGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGtat                              <  1:3876108/93‑1 (MQ=255)
                                                             cGTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATaa                            >  2:3135753/1‑93 (MQ=255)
                                                             cGTGAGGAAATGCTGCATTGCTGACATTTCCGCCAGCAATGCAACGTCAAAGACCTGTCTAGAACAAACGATAACCGTAGCCCCATAGTATaa                            <  2:1203181/93‑1 (MQ=255)
                                                              ttGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAAc                           >  2:4818861/2‑93 (MQ=255)
                                                              gTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAAc                           <  2:6564043/93‑1 (MQ=255)
                                                              gTGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAAc                           >  1:8602592/1‑93 (MQ=255)
                                                               tGAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACg                          >  2:5451071/1‑93 (MQ=255)
                                                                gAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACgg                         >  1:3192292/1‑93 (MQ=255)
                                                                gAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACgg                         >  1:6651955/1‑93 (MQ=255)
                                                                gAGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACgg                         >  2:2029019/1‑93 (MQ=255)
                                                                 aGGAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGt                        <  2:3995033/93‑1 (MQ=255)
                                                                  ggAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGtt                       <  2:5323420/93‑1 (MQ=255)
                                                                  ggAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGtt                       <  1:5259620/93‑1 (MQ=255)
                                                                  ggAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGtt                       <  1:968007/93‑1 (MQ=255)
                                                                   gcaaTGCTGCATTGCTGACCTTTGCGCAAGCAATGCAACGTCAACGAACTTTCTAGAACAAACGATAACCGGAGCCCCATAGTATAGAGGttg                      >  1:7750349/3‑92 (MQ=255)
                                                                   gAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTa                      <  2:801177/93‑1 (MQ=255)
                                                                   gAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTa                      <  1:7015932/93‑1 (MQ=255)
                                                                   gAAATGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTa                      >  2:7549990/1‑93 (MQ=255)
                                                                    aaaTGCTTCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAg                     <  1:8316587/93‑1 (MQ=255)
                                                                    aaaTGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAg                     <  2:7746412/93‑1 (MQ=255)
                                                                    aaaTGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAg                     <  1:4114904/93‑1 (MQ=255)
                                                                    aaaTGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAg                     <  1:2565766/93‑1 (MQ=255)
                                                                     aaTGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAgg                    <  1:1688643/93‑1 (MQ=255)
                                                                      aTGCTGCATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAggg                   >  1:3786808/1‑93 (MQ=255)
                                                                       tgctgcATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGc                  <  1:5116678/93‑1 (MQ=255)
                                                                         ctgcATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGa                >  1:7699397/1‑93 (MQ=255)
                                                                           gcATTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagc              <  1:6096103/93‑1 (MQ=255)
                                                                             aTTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagcag            >  1:4589002/1‑93 (MQ=255)
                                                                             aTTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagcag            <  1:278572/93‑1 (MQ=255)
                                                                             aTTGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagcag            >  1:7508457/1‑93 (MQ=255)
                                                                              ttGCTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTGAGGGCGagcagc           >  2:8053042/1‑93 (MQ=255)
                                                                                gCTGACATTTGCGCCAGCAATGCAACATCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCAc         <  1:2234391/93‑1 (MQ=255)
                                                                                 cTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACt        <  2:675539/93‑1 (MQ=255)
                                                                                 cTGACATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACt        <  1:8200163/93‑1 (MQ=255)
                                                                                    aCATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTcc     >  1:10610/1‑93 (MQ=255)
                                                                                     cATTTGCGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCa    >  1:4502537/1‑93 (MQ=255)
                                                                                          |                                                                                         
GATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGT  >  NC_000913/2222945‑2223124

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑