Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,644,603 (G)6→5 coding (743/2313 nt) pbpC ← putative peptidoglycan enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,644,5980G.99.5% 917.3 / ‑3.5 200coding (748/2313 nt)pbpCputative peptidoglycan enzyme
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (89/110);  total (90/110)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 228 total aligned reads displayed.

AAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTT  >  NC_000913/2644506‑2644685
                                                                                            |                                                                                       
aaGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCg                                                                                         >  2:712832/1‑93 (MQ=255)
aaGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCg                                                                                         >  1:5095436/1‑93 (MQ=255)
  gACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggg                                                                                       <  1:6711672/93‑1 (MQ=255)
  gACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggg                                                                                       >  2:517417/1‑93 (MQ=255)
  gACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggg                                                                                       >  2:3363235/1‑93 (MQ=255)
  gACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggg                                                                                       >  1:202736/1‑93 (MQ=255)
   aCCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggg                                                                                      <  2:2991744/93‑1 (MQ=255)
   aCCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATCTGTCgggg                                                                                      >  1:478554/1‑93 (MQ=255)
      ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                   >  2:1020540/1‑91 (MQ=255)
      ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                   >  1:1032744/1‑91 (MQ=255)
      ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                   <  2:5468112/93‑3 (MQ=255)
      ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                   <  1:5649424/93‑3 (MQ=255)
       gCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGgcc                                                                                  >  2:1571775/1‑93 (MQ=255)
       gCATCCAACGTAGTAGTGATTTGGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcca                                                                                  <  1:3267953/93‑4 (MQ=255)
        cATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggccag                                                                                 >  1:6414817/1‑89 (MQ=255)
          tCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGgccagcc                                                                              <  2:6253040/93‑1 (MQ=255)
           ccAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGgccagcca                                                                             >  2:4270027/1‑93 (MQ=255)
           ccAACGTAGTAGTGATTTCGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGgccagcca                                                                             <  2:5027452/93‑1 (MQ=255)
            cAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGccagccag                                                                            <  1:3022291/93‑1 (MQ=255)
              aCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGAt                                                                          <  1:1273696/93‑1 (MQ=255)
               cGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATg                                                                         <  1:5564281/93‑1 (MQ=255)
                  agtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGt                                                                      <  2:2270170/93‑1 (MQ=255)
                   gtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGtt                                                                     >  2:2000906/1‑93 (MQ=255)
                   gtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGtt                                                                     >  1:2757346/1‑93 (MQ=255)
                   gtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGtt                                                                     <  1:3436016/93‑1 (MQ=255)
                     agtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttct                                                                   <  2:2060967/93‑1 (MQ=255)
                      gtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttctt                                                                  >  2:2182624/1‑93 (MQ=255)
                      gtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttctt                                                                  <  2:6413095/93‑1 (MQ=255)
                      gtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttctt                                                                  <  1:3844634/93‑1 (MQ=255)
                        gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTcc                                                                <  2:5718659/93‑1 (MQ=255)
                        gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTcc                                                                <  2:1312417/93‑1 (MQ=255)
                        gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTcc                                                                >  1:2360765/1‑93 (MQ=255)
                        gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTcc                                                                <  1:1857812/93‑1 (MQ=255)
                        gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTcc                                                                >  2:66400/1‑93 (MQ=255)
                        gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTcc                                                                >  2:371782/1‑93 (MQ=255)
                         aTTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTccc                                                               <  1:2405202/93‑1 (MQ=255)
                          ttttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGTATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCt                                                              <  2:1903607/93‑1 (MQ=255)
                          ttttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCt                                                              >  1:6208692/1‑93 (MQ=255)
                          ttttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCt                                                              <  2:4336207/93‑1 (MQ=255)
                           tttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCtt                                                             <  1:6036731/93‑1 (MQ=255)
                           tttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCtt                                                             >  2:2819995/1‑93 (MQ=255)
                              gTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGAc                                                          <  1:3787603/93‑1 (MQ=255)
                                 gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGActct                                                       >  1:5785750/1‑93 (MQ=255)
                                 gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGActct                                                       >  1:4312137/1‑93 (MQ=255)
                                 gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGActct                                                       <  1:1789548/93‑1 (MQ=255)
                                  cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtt                                                      >  2:5083261/1‑93 (MQ=255)
                                  cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtt                                                      <  2:454346/93‑1 (MQ=255)
                                  cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtt                                                      <  1:3205304/93‑1 (MQ=255)
                                   ttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCttt                                                     >  1:4289436/1‑93 (MQ=255)
                                   ttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCttt                                                     >  2:3503703/1‑93 (MQ=255)
                                   ttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCttt                                                     <  2:3870324/93‑1 (MQ=255)
                                   ttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCttt                                                     >  2:3856306/1‑93 (MQ=255)
                                     ttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTa                                                   <  1:199102/93‑1 (MQ=255)
                                       gcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAcc                                                 <  1:2471405/93‑1 (MQ=255)
                                       gcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAcc                                                 <  1:6648090/93‑1 (MQ=255)
                                       gcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAcc                                                 <  2:4893928/93‑1 (MQ=255)
                                       gcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAcc                                                 <  1:1006517/93‑1 (MQ=255)
                                          tttACCGAGCATCATGCGCGAAAACAGCGGTGCTAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGc                                              <  2:1182476/93‑1 (MQ=255)
                                            tACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTc                                            >  1:507447/1‑93 (MQ=255)
                                             aCCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCa                                           >  1:6003845/1‑93 (MQ=255)
                                               cGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACg                                         <  2:1445496/93‑1 (MQ=255)
                                                gAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACgg                                        >  1:2944176/1‑93 (MQ=255)
                                                gAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACgg                                        >  1:6476947/1‑93 (MQ=255)
                                                gAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACgg                                        >  2:4472477/1‑93 (MQ=255)
                                                  gCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACgggg                                      >  2:3244223/1‑92 (MQ=255)
                                                  gCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGa                                      <  2:6310656/93‑1 (MQ=255)
                                                    agcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGAcc                                    >  1:5309663/3‑93 (MQ=255)
                                                     tcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACca                                   >  2:6392458/1‑93 (MQ=255)
                                                      catGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcac                                  <  1:96591/93‑1 (MQ=255)
                                                       atGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcaca                                 >  2:2374042/1‑93 (MQ=255)
                                                       atGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcaca                                 >  2:6493583/1‑93 (MQ=255)
                                                        tGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacac                                <  1:842201/93‑1 (MQ=255)
                                                        tGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacac                                >  1:6255215/1‑93 (MQ=255)
                                                             gAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCtt                           <  2:5882833/93‑1 (MQ=255)
                                                             gAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCGCACGGGACCACACACCtt                           <  2:4068393/93‑1 (MQ=255)
                                                              aaaaCAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTg                          >  1:4739608/1‑93 (MQ=255)
                                                              aaaaCAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTg                          >  1:6407763/1‑93 (MQ=255)
                                                               aaaCAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCNTTTACCTGCTCACGGGACCAGACNCCTTGc                         >  2:505236/1‑93 (MQ=255)
                                                                aaCAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCa                        >  2:2577/1‑93 (MQ=255)
                                                                  cAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACg                      >  2:3636665/1‑93 (MQ=255)
                                                                   aGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACgg                     <  2:5597828/93‑1 (MQ=255)
                                                                     cGGTGCCAGTTGCGGCCTTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCCCACACTTTGCACGGcc                   >  1:2177338/1‑93 (MQ=255)
                                                                     cGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGcc                   <  2:6168869/93‑1 (MQ=255)
                                                                      ggTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCa                  <  1:5396471/93‑1 (MQ=255)
                                                                       gTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCAt                 <  2:4672663/93‑1 (MQ=255)
                                                                       gTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCAt                 >  1:4910073/1‑93 (MQ=255)
                                                                         gCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATcc               <  2:6289965/93‑1 (MQ=255)
                                                                           cAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt             <  1:244353/93‑1 (MQ=255)
                                                                           cAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt             <  1:4692366/93‑1 (MQ=255)
                                                                           cAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt             <  1:1629303/93‑1 (MQ=255)
                                                                           cAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt             >  1:2209650/1‑93 (MQ=255)
                                                                            cgtggcgcCATTTGTCGGGATCCAGCAAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGtt            <  1:4443201/89‑1 (MQ=255)
                                                                            aGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGtt            <  1:1244124/93‑1 (MQ=255)
                                                                             gTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTc           >  1:4481842/1‑93 (MQ=255)
                                                                             gTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTc           >  1:3306108/1‑93 (MQ=255)
                                                                               tGCGGCATTGGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCTATCCGTTCGa         <  1:5146964/93‑1 (MQ=255)
                                                                                gCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGac        >  2:4577172/1‑92 (MQ=255)
                                                                                   gCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCAc     <  1:6580286/93‑1 (MQ=255)
                                                                                    cATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACt    <  2:3131729/93‑1 (MQ=255)
                                                                                     aTTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACtt   >  1:408212/1‑93 (MQ=255)
                                                                                            |                                                                                       
AAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTT  >  NC_000913/2644506‑2644685

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