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.97.3% 482.6 / 3.5 112coding (748/2313 nt)pbpCputative peptidoglycan enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base . (58/51);  total (58/54)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 143 total aligned reads displayed.

AGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTA  >  NC_000913/2644507‑2644686
                                                                                           |                                                                                        
aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACATCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         <  1:3428264/93‑1 (MQ=255)
aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         >  1:3720070/1‑93 (MQ=255)
aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         <  2:1197885/93‑1 (MQ=255)
aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         <  2:1419933/93‑1 (MQ=255)
aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         <  2:2031978/93‑1 (MQ=255)
aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         <  2:3101797/93‑1 (MQ=255)
 gACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggg                                                                                        <  1:1002064/93‑1 (MQ=255)
    cGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggc                                                                                     >  1:1590889/1‑92 (MQ=255)
    cGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggc                                                                                     <  2:2760509/93‑2 (MQ=255)
     ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                    >  1:1055407/1‑91 (MQ=255)
     ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                    <  1:3822277/93‑3 (MQ=255)
     ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                    <  2:1440965/93‑3 (MQ=255)
     ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                    >  1:1055408/1‑91 (MQ=255)
      gCATCCAACGGAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGgcc                                                                                   >  1:2979308/1‑93 (MQ=255)
       cATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggccag                                                                                  >  1:377199/1‑89 (MQ=255)
        aTCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGgccagc                                                                                >  2:2859907/1‑93 (MQ=255)
            aaCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGa                                                                            >  1:3680293/1‑93 (MQ=255)
            aaCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGa                                                                            <  2:4101626/93‑1 (MQ=255)
            aaCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGa                                                                            <  2:773767/93‑1 (MQ=255)
             aCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGAt                                                                           <  2:1701720/93‑1 (MQ=255)
              cGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATg                                                                          <  2:1342455/93‑1 (MQ=255)
              cGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATg                                                                          <  1:690181/93‑1 (MQ=255)
               gtagtaGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATgg                                                                         >  2:4269604/1‑93 (MQ=255)
               gtagtaGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATgg                                                                         >  2:2896713/1‑93 (MQ=255)
               gtagtaGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATgg                                                                         <  1:1342332/93‑1 (MQ=255)
                tagtagTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATggg                                                                        <  2:3227684/93‑1 (MQ=255)
                tagtagTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATggg                                                                        <  1:1280440/93‑1 (MQ=255)
                  gtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGtt                                                                      >  1:1680543/1‑93 (MQ=255)
                   tagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGCTGGGttc                                                                     <  2:2882893/93‑1 (MQ=255)
                   tagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttc                                                                     <  2:2777275/93‑1 (MQ=255)
                    agtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATTGGgtct                                                                    >  1:277879/1‑93 (MQ=255)
                    agtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttct                                                                    >  1:918077/1‑93 (MQ=255)
                    agtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttct                                                                    >  1:3323111/1‑93 (MQ=255)
                         ttttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCt                                                               <  2:4177999/93‑1 (MQ=255)
                            tGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGACAGCCAGATGGGTTCTTCCCTTGa                                                            >  1:4389936/1‑93 (MQ=255)
                             gTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGAc                                                           >  2:3312697/1‑93 (MQ=255)
                               cGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGAct                                                          <  1:817496/93‑1 (MQ=255)
                                gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGActct                                                        >  1:3087099/1‑93 (MQ=255)
                                gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGActct                                                        >  2:4730801/1‑93 (MQ=255)
                                gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTACCTTGActct                                                        <  2:840664/93‑1 (MQ=255)
                                 cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtt                                                       >  2:1035248/1‑93 (MQ=255)
                                 cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtt                                                       >  1:2250242/1‑93 (MQ=255)
                                      tcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAcc                                                  <  1:3564927/92‑1 (MQ=255)
                                      gcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAcc                                                  >  1:2390766/1‑93 (MQ=255)
                                       cttttGCCGCGCATCATGCGCGAAAACAGCAGTGACGGTTAAAGCTTTTGTCGGGAGCAAGCCAGATAAGTTCTACCCTTGCCTCTATAACCt                                                 <  2:3600520/93‑1 (MQ=255)
                                            tccGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCa                                            <  1:3945824/92‑1 (MQ=255)
                                            aCCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCa                                            >  2:3287253/1‑93 (MQ=255)
                                            aCCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCa                                            >  1:3849524/1‑93 (MQ=255)
                                            aCCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCa                                            <  1:341670/93‑1 (MQ=255)
                                             ccGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCAc                                           <  1:418527/93‑1 (MQ=255)
                                             ccGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCAc                                           <  2:2457215/93‑1 (MQ=255)
                                               gAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACgg                                         >  1:3528800/1‑93 (MQ=255)
                                               caGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTTTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACgg                                         <  2:1601015/92‑1 (MQ=255)
                                                 gCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGa                                       <  1:3992004/93‑1 (MQ=255)
                                                   atcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGAcc                                     <  1:3089513/93‑1 (MQ=255)
                                                   atcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGAcc                                     >  2:3469100/1‑93 (MQ=255)
                                                    tcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACca                                    <  2:4586236/93‑1 (MQ=255)
                                                    tcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACca                                    <  1:4118621/93‑1 (MQ=255)
                                                     catGCGCGAAAACAGCGGTGCCAGTTGCGGAATTTGTCGGGGGCCATCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcac                                   <  1:2739943/93‑1 (MQ=255)
                                                       tGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacac                                 <  1:2826900/93‑1 (MQ=255)
                                                       tGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacac                                 >  1:2413922/1‑93 (MQ=255)
                                                          gcgAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACAcc                              >  1:3139922/1‑93 (MQ=255)
                                                          gcgAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACAcc                              >  2:4065692/1‑93 (MQ=255)
                                                          gcgAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACAcc                              >  2:2585666/1‑93 (MQ=255)
                                                           cgAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCt                             <  1:2659393/93‑1 (MQ=255)
                                                           cgAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCt                             <  2:4039769/93‑1 (MQ=255)
                                                            gAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCtt                            <  1:2625181/93‑1 (MQ=255)
                                                            gAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCtt                            <  2:4533791/93‑1 (MQ=255)
                                                               aaCAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCa                         <  2:2611637/93‑1 (MQ=255)
                                                               aaCAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCa                         >  1:3413982/1‑93 (MQ=255)
                                                                 cAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACg                       >  2:1060164/1‑93 (MQ=255)
                                                                 cAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACg                       >  1:947383/1‑93 (MQ=255)
                                                                 aaGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACg                       <  1:445736/92‑1 (MQ=255)
                                                                  aGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACg                       <  2:4005545/93‑1 (MQ=255)
                                                                   gCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACTTGCCCACGGGACCACACACCTTGCAAGGc                     >  2:3764364/1‑93 (MQ=255)
                                                                   gCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACGCACCTTGCACGGc                     >  1:2914040/1‑93 (MQ=255)
                                                                   gCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGc                     >  1:4702766/1‑93 (MQ=255)
                                                                    cGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGcc                    <  2:2208757/93‑1 (MQ=255)
                                                                    cGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGcc                    <  2:1817418/93‑1 (MQ=255)
                                                                     ggTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCa                   >  2:77853/1‑93 (MQ=255)
                                                                     ggTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCa                   >  2:522454/1‑93 (MQ=255)
                                                                       tGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATc                 >  2:3266567/1‑93 (MQ=255)
                                                                       tGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATc                 <  1:1645968/93‑1 (MQ=255)
                                                                       tGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATc                 >  2:3249122/1‑93 (MQ=255)
                                                                       tGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATc                 <  2:3061541/93‑1 (MQ=255)
                                                                         ccAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGCACCACACACCTTGCACGGCCATCCg               >  1:2484428/1‑93 (MQ=255)
                                                                          cAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt              >  1:2917083/1‑93 (MQ=255)
                                                                          cAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt              <  2:1122840/93‑1 (MQ=255)
                                                                          cAGTTGCAGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt              >  1:3220427/1‑93 (MQ=255)
                                                                           aGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGtt             >  2:4617799/1‑93 (MQ=255)
                                                                           aGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGtt             <  2:1976334/93‑1 (MQ=255)
                                                                            gTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTc            >  1:1840711/1‑93 (MQ=255)
                                                                            gTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTc            >  1:4097824/1‑93 (MQ=255)
                                                                            gTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTc            >  1:4097825/1‑93 (MQ=255)
                                                                            gTTGCGGCAGTTGTCGGGGGCCAGCCAGATGGGTTCTGCCATTGACTCTTTTACCCGCTCACGGGACCACACACCTTGCCAGGCCATCCGTTc            <  2:3638960/93‑1 (MQ=255)
                                                                              tGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGa          <  2:691450/93‑1 (MQ=255)
                                                                                    aTTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACtt    >  1:1740866/1‑93 (MQ=255)
                                                                                    aTTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACtt    >  1:3364333/1‑93 (MQ=255)
                                                                                     tttGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACttt   >  1:417966/1‑93 (MQ=255)
                                                                                      ttGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTtta  <  2:4245085/93‑1 (MQ=255)
                                                                                           |                                                                                        
AGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTA  >  NC_000913/2644507‑2644686

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