Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,151,969 T→A 7.8% R168S (AGA→AGT pphC ← protein‑serine/threonine phosphatase PphC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,9690TA 7.8% 176.8 / 3.8 59R168S (AGA→AGTpphCprotein‑serine/threonine phosphatase PphC
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base A (2/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGA  >  NC_000913/2151823‑2152118
                                                                                                                                                  |                                                                                                                                                     
tGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTg                                                                                                                                                    >  2:308356/1‑150 (MQ=255)
       ctgcCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGATGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACgg                                                                                                                                             <  1:419195/150‑1 (MQ=255)
         gccgccAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttc                                                                                                                                           <  2:127350/150‑5 (MQ=255)
          ccgccAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttca                                                                                                                                          <  2:119626/150‑6 (MQ=255)
            gccAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCt                                                                                                                                        >  2:988391/1‑150 (MQ=255)
            gccAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttcacc                                                                                                                                        <  2:1220839/150‑8 (MQ=255)
               aGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcg                                                                                                                                     <  2:667556/150‑1 (MQ=255)
                 tCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcgtc                                                                                                                                   <  1:995586/150‑1 (MQ=255)
                    aTTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGt                                                                                                                                >  2:1026214/1‑150 (MQ=255)
                       gAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAAt                                                                                                                             >  2:1329953/1‑150 (MQ=255)
                       gAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGAATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAAt                                                                                                                             <  2:1391469/150‑1 (MQ=255)
                        aaaaaCGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAAc                                                                                                                                               <  2:998567/131‑1 (MQ=255)
                        aaaaaCGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAAc                                                                                                                                               >  1:998567/1‑131 (MQ=255)
                              ggggTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGc                                                                                                                      >  1:833361/1‑150 (MQ=255)
                                ggTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGt                                                                                                                    >  1:408633/1‑150 (MQ=255)
                                  tAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCa                                                                                                                  >  1:121487/1‑150 (MQ=255)
                                        aagGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCg                                                                                                                     <  1:439459/141‑1 (MQ=255)
                                        aagGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCg                                                                                                                     >  2:439459/1‑141 (MQ=255)
                                         agGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGc                                                                                                           <  2:121487/150‑1 (MQ=255)
                                          ggCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCa                                                                                                          <  2:875073/150‑1 (MQ=255)
                                          ggCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCa                                                                                                          <  2:1049736/150‑1 (MQ=255)
                                            cacaTGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCata                                                                                                        >  1:309188/1‑150 (MQ=255)
                                              caTGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTCCTGGTGAAACTTTCCAGACTGGAAACGGCAACTTCCACGCTAATGAAGTGCCTCATCTACGCatata                                                                                                      >  2:449044/1‑149 (MQ=255)
                                              caTGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATCCTGCAACATTGTGCGCACGCCCAGTGATGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATAtt                                                                                                      >  1:1132770/1‑150 (MQ=255)
                                                tGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCa                                                                                                    >  2:745546/1‑150 (MQ=255)
                                                           aTCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCgg                                                                                         >  1:541388/1‑150 (MQ=255)
                                                               aGCATATTTAACGCCTGTCGCTGGATAGCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCTTCGGTAATGAAGTGCGTCATTTTGGCATATTCACCCGCCATCGGCGTg                                                                                     <  2:1132770/150‑1 (MQ=255)
                                                               aGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCg                                                                                                                     <  1:188879/118‑1 (MQ=255)
                                                               aGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCg                                                                                                                     >  2:188879/1‑118 (MQ=255)
                                                                 cATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgag                                                                                   <  2:309188/150‑1 (MQ=255)
                                                                      ttAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCa                                                                              <  1:744027/150‑1 (MQ=255)
                                                                           gCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTCCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGc                                                                         <  1:439150/150‑1 (MQ=255)
                                                                               gTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGAc                                                                     <  1:745546/150‑1 (MQ=255)
                                                                                  gCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGt                                                                  >  2:884501/1‑150 (MQ=255)
                                                                                   cTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTg                                                                 <  2:207789/150‑1 (MQ=255)
                                                                                     ggATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCg                                                                   <  1:272245/146‑1 (MQ=255)
                                                                                     ggATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCg                                                                   >  2:272245/1‑146 (MQ=255)
                                                                                          ccATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCaaa                                                          >  1:572597/1‑150 (MQ=255)
                                                                                          ccATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACAGTGACCCaaa                                                          >  1:770010/1‑150 (MQ=255)
                                                                                           cATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAAt                                                         >  1:243839/1‑150 (MQ=255)
                                                                                            aTCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGTCAATTGCAGACCGTGACCCAAATc                                                        >  1:946690/1‑150 (MQ=255)
                                                                                            aTCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATc                                                        >  1:1283741/1‑150 (MQ=255)
                                                                                               cGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTac                                                     <  2:60128/150‑1 (MQ=255)
                                                                                               cGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGGCATGGTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTac                                                     <  2:596922/150‑1 (MQ=255)
                                                                                                gTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTacc                                                    <  1:175218/150‑1 (MQ=255)
                                                                                                      gctgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCAcg                                              >  1:418358/1‑150 (MQ=255)
                                                                                                           aaCTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACgccgcc                                         <  2:842541/150‑1 (MQ=255)
                                                                                                              tttGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATc                                      >  2:879973/1‑150 (MQ=255)
                                                                                                                 gtgCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGcc                                   <  1:1329953/150‑1 (MQ=255)
                                                                                                                 gtgCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGcc                                   <  1:54890/150‑1 (MQ=255)
                                                                                                                    cgcACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAAt                                >  1:555367/1‑150 (MQ=255)
                                                                                                                     gcACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATc                               >  1:764696/1‑150 (MQ=255)
                                                                                                                                       gccgTTTCCAGACTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCg             >  1:1329372/4‑150 (MQ=255)
                                                                                                                                          gTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCa          <  1:928763/150‑1 (MQ=255)
                                                                                                                                           tttCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAg         <  1:308356/150‑1 (MQ=255)
                                                                                                                                            ttCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAgg        <  2:541388/150‑1 (MQ=255)
                                                                                                                                               cAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGa     >  1:578996/1‑150 (MQ=255)
                                                                                                                                               cAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGa     >  2:251312/1‑150 (MQ=255)
                                                                                                                                                 gTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATg   <  2:555367/150‑1 (MQ=255)
                                                                                                                                                  tCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGa  <  2:833361/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGA  >  NC_000913/2151823‑2152118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: