Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,361,762 C→G 6.2% G720R (GGC→CGC)  PP_2076 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,361,7620CG 6.2% 201.5 / 3.8 64G720R (GGC→CGC) PP_2076hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/25);  new base G (2/2);  total (37/27)
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.81e-01

CACTCTGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCCA  >  NC_002947/2361615‑2361911
                                                                                                                                                   |                                                                                                                                                     
cACTCTGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGt                                                                                                                                                     >  1:522612/1‑150 (MQ=255)
    ctGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcg                                                                                                                                                 >  1:200263/1‑150 (MQ=255)
       gctgctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                              >  2:939346/1‑150 (MQ=255)
        ctgctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgcg                                                                                                                                             >  1:1125803/1‑150 (MQ=255)
             tttGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGaccac                                                                                                                                        >  1:376447/1‑150 (MQ=255)
             tttGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGaccac                                                                                                                                        <  2:676546/150‑1 (MQ=255)
              ttGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCGGTACGGCGCGaacagt                                                                                                                                       <  2:724770/150‑3 (MQ=255)
                gTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTg                                                                                                                                     >  1:918377/1‑150 (MQ=255)
                 tAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTgg                                                                                                                                    >  2:843168/1‑150 (MQ=255)
                   aTGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTc                                                                                                                                  <  2:248189/150‑1 (MQ=255)
                      cTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCcac                                                                                                                               >  1:120067/1‑150 (MQ=255)
                       tCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCcaca                                                                                                                              <  2:918377/150‑1 (MQ=255)
                       tCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCcaca                                                                                                                              >  1:1026227/1‑150 (MQ=255)
                        cAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACAt                                                                                                                             >  2:575529/1‑150 (MQ=255)
                         aaTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATa                                                                                                                            >  2:534633/1‑150 (MQ=255)
                         aaTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATa                                                                                                                            >  2:1141929/1‑150 (MQ=255)
                          aTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGGTCGCGCGGTACGGCGCGaacagt                                                                                                                                       >  2:430157/1‑133 (MQ=21)
                          aTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGGTCGCGCGGTACGGCGCGaacagt                                                                                                                                       <  1:430157/138‑3 (MQ=21)
                            cGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACt                                                                                                                         <  1:575529/150‑1 (MQ=255)
                                  aGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTcc                                                                                                                   >  1:616564/1‑150 (MQ=255)
                                   gTGAAATCCCACCACGCGAGTACCAACCCCCCGGTCCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCCGGTCCACATAACTGGCTCCa                                                                                                                  >  1:876756/1‑150 (MQ=255)
                                     gAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAAc                                                                                                                <  1:534633/150‑1 (MQ=255)
                                            caccacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGGTACCACCTGGTCCACATAACTGGCTCCAACTTCGAAg                                                                                                         >  1:653931/1‑150 (MQ=255)
                                            caccacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAg                                                                                                         >  2:164949/1‑150 (MQ=255)
                                                    gAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCcgccgc                                                                                                 >  2:942308/1‑150 (MQ=255)
                                                           aaCCCCCCGGTGCCAGCTAAACTTGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTa                                                                                          <  1:1141929/150‑1 (MQ=255)
                                                           aaCCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTa                                                                                          <  2:427734/150‑1 (MQ=255)
                                                           aaCCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTa                                                                                          >  1:423932/1‑150 (MQ=255)
                                                                   ggTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATcgc                                                                                  <  2:522612/150‑1 (MQ=255)
                                                                     tGCCAGCTAAACTGGCGCCCCATAGTTTCGAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCgg                                                                                >  1:1188788/1‑150 (MQ=255)
                                                                     tGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGGACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCgg                                                                                <  2:593525/150‑1 (MQ=255)
                                                                          gtttagcTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTa                                                                                          <  1:1313835/129‑1 (MQ=255)
                                                                          gtttagcTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTa                                                                                          >  2:1313835/7‑135 (MQ=255)
                                                                                   gcgcCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCg                                                                  <  1:939346/150‑1 (MQ=255)
                                                                                   gcgcCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGCCATCGCGGCCGCCAGGAGTTCg                                                                  <  2:1125803/150‑1 (MQ=255)
                                                                                        ccATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTACTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTAGGTCtt                                                             >  2:627864/1‑150 (MQ=255)
                                                                                         aaTAGTATCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTc                                                            >  2:708037/2‑150 (MQ=255)
                                                                                                         cTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGt                                            >  2:1178840/1‑150 (MQ=40)
                                                                                                           agagCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCa                                          <  2:423932/150‑1 (MQ=40)
                                                                                                              gCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCaa                                                                                                                 <  2:828869/76‑1 (MQ=255)
                                                                                                              gCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCaa                                                                                                                 >  1:828869/1‑76 (MQ=255)
                                                                                                                 gCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTGCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGt                                    >  1:1244394/1‑150 (MQ=255)
                                                                                                                 gCCGCTGCTCCTCGTCCGGTCACACTGTTCGCACCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGt                                    >  1:205069/1‑150 (MQ=255)
                                                                                                                      tgctCCTCGTCCGGTAACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCTCGCATCACTTTCAGACAGTTcgcg                               >  2:458761/1‑150 (MQ=255)
                                                                                                                       gctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCACCCGTAGTCATCGCGGCTGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgcgc                              <  2:1026227/150‑1 (MQ=255)
                                                                                                                              gTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTc                                                               >  2:322488/1‑110 (MQ=255)
                                                                                                                              gTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTc                                                               <  1:322488/110‑1 (MQ=255)
                                                                                                                              gTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGa                       >  1:836162/1‑150 (MQ=255)
                                                                                                                               tCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGaa                      <  2:836162/150‑1 (MQ=255)
                                                                                                                               tCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGaa                      <  1:123003/150‑1 (MQ=255)
                                                                                                                                ccGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAg                     <  1:164949/150‑1 (MQ=255)
                                                                                                                                      acacTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATcgc                                                                                  >  1:1196059/1‑83 (MQ=40)
                                                                                                                                      acacTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATcgc                                                                                  <  2:1196059/83‑1 (MQ=40)
                                                                                                                                      acacTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATc               >  1:781438/1‑150 (MQ=40)
                                                                                                                                        acTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACt                                              >  2:1190592/1‑117 (MQ=40)
                                                                                                                                        acTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACt                                              <  1:1190592/117‑1 (MQ=40)
                                                                                                                                        acTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAg             <  1:1178840/150‑1 (MQ=255)
                                                                                                                                           ggtCGCGCGGTACGGCGCGAACACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCAt          >  2:509735/3‑150 (MQ=255)
                                                                                                                                            ttCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGcc                                                                                             >  2:163475/1‑66 (MQ=255)
                                                                                                                                            ttCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGcc                                                                                             <  1:163475/66‑1 (MQ=255)
                                                                                                                                            ttCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgcg                               >  2:366075/1‑128 (MQ=255)
                                                                                                                                            ttCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgcg                               <  1:366075/128‑1 (MQ=255)
                                                                                                                                               gcgcCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATcccc      <  2:120067/150‑1 (MQ=255)
                                                                                                                                                cgcCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGCCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCt     >  1:823122/1‑150 (MQ=255)
                                                                                                                                                   cGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCca  <  2:1319206/150‑1 (MQ=40)
                                                                                                                                                   cGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCca  <  2:205069/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CACTCTGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCCA  >  NC_002947/2361615‑2361911

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