Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,310 T→C 10.8% S159P (TCT→CCT)  PP_3774 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,3100TC10.8% 162.5 / 10.0 65S159P (TCT→CCT) PP_3774hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base C (3/4);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

CGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCC  >  NC_002947/4301162‑4301454
                                                                                                                                                    |                                                                                                                                                
cgacgGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGt                                                                                                                                                  <  3:329783/149‑1 (MQ=255)
            ggTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTgg                                                                                                                                      <  5:271727/149‑1 (MQ=255)
              tAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGTGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGCCTTCAGACGTGGtg                                                                                                                                    <  7:396545/149‑2 (MQ=255)
                aaGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAAc                                                                                                                                  <  8:65031/149‑1 (MQ=255)
                aaGGGGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGCCTTCAGACGTGGtgtg                                                                                                                                  <  7:82014/149‑4 (MQ=255)
                 aGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGGGCGACGTTCAACATCATCATTCACACCACGTCTGAAGGCGTGGTAAcc                                                                                                                                 <  4:92071/149‑1 (MQ=255)
                   gtgtTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGt                                                                                                                               <  6:158134/149‑1 (MQ=255)
                     gtTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGa                                                                                                                             <  8:47838/149‑1 (MQ=255)
                            aaCAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGa                                                                                                                      <  2:352177/149‑1 (MQ=255)
                             aCAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGaa                                                                                                                     >  1:328027/1‑149 (MQ=255)
                              cAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAAt                                                                                                                    <  2:328027/149‑1 (MQ=255)
                              cAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAAt                                                                                                                    <  2:5411/149‑1 (MQ=255)
                              cAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGAGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACAGTGACGGCCGAAt                                                                                                                    >  1:182251/1‑149 (MQ=255)
                                 gggCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGc                                                                                                                  >  6:16578/1‑148 (MQ=255)
                                     aaCATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGCCTTCAGACGTGGtg                                                                                                                                    <  2:223672/126‑2 (MQ=255)
                                     aaCATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGCCTTCAGACGTGGtg                                                                                                                                    <  2:230604/126‑2 (MQ=255)
                                     aaCATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGCCTTCAGACGTGGtg                                                                                                                                    >  1:230604/1‑125 (MQ=255)
                                     aaCATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGCCTTCAGACGTGGtg                                                                                                                                    >  1:223672/1‑125 (MQ=255)
                                      ccATGACCCTGGGGTTGCGGCGCACCTGCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGGTCAACATCAGCATTCACCCCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACggg                                                                                                            <  8:295619/148‑1 (MQ=255)
                                           aCCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATc                                                                                                       <  2:96152/149‑1 (MQ=255)
                                             ccTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATcc                                                                                                      >  7:124028/1‑148 (MQ=255)
                                               tGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTg                                                                                                                                       <  6:269638/113‑1 (MQ=255)
                                               tGGGGATTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAACGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGAAACCGTGACGCCCGAAACCAACGGGGTATCCGc                                                                                                    >  3:6595/1‑148 (MQ=255)
                                                 gggTTTCGGCGCACCTTCTAGCAGAGCGACCCGGTGAATTCCACAGTCCTGAATGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCAGACCACGTCTGAAGGCGTTGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTtc                                                                                                 <  5:380054/149‑1 (MQ=255)
                                                  ggTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTtcg                                                                                                >  1:165720/1‑149 (MQ=255)
                                                        ggCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGc                                                                                          <  1:51501/149‑1 (MQ=255)
                                                           gcACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGAt                                                                                       <  5:139804/149‑1 (MQ=255)
                                                            cACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTTCCACTGGGTGCAACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGAtt                                                                                      <  8:9363/149‑1 (MQ=255)
                                                                ttCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTg                                                                                    >  1:324097/1‑147 (MQ=255)
                                                                ttCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTg                                                                                    >  3:323187/1‑147 (MQ=255)
                                                                        gaccgaccCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACa                                                                          <  3:357863/149‑1 (MQ=255)
                                                                          ccgaccAGTGGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCCCACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCTCTTCGTCGGGCGATTTGCAGCTCGACAGt                                                                        <  4:264699/149‑1 (MQ=255)
                                                                             accCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTtcgtcg                                                                                             >  4:2812/1‑125 (MQ=255)
                                                                             accCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTtcgtcg                                                                                             <  3:2812/125‑1 (MQ=255)
                                                                             accCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTAcc                                                                     >  7:297103/1‑149 (MQ=255)
                                                                                  gTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGc                                                                >  7:238254/1‑149 (MQ=255)
                                                                                  gTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGc                                                                >  6:350229/1‑149 (MQ=255)
                                                                                     aaTTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTc                                                             <  5:16578/149‑1 (MQ=255)
                                                                                      aTTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTTGGCTTc                                                             <  8:238254/148‑1 (MQ=255)
                                                                                       ttCCACAGTCCTGAAGGGGGTGCCTCTGGGTGCGACGTTCAACTTCAGGATTCACTCGATGTGTGTAGGCGTGGTAACCGTGACGGGCGAATGGAACGGGGTATCGGCTTCGTGGGGCGATTTGGAGCTGGACAGTACCTGGGCTTctc                                                           <  4:6595/149‑3 (MQ=255)
                                                                                       ttCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCa                                                            <  1:133403/148‑1 (MQ=255)
                                                                                       ttCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCa                                                            >  2:133403/1‑148 (MQ=255)
                                                                                       gtCCACAGTCCTGCAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCTTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCa                                                            <  8:401168/147‑1 (MQ=255)
                                                                                              gTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGtt                                                    <  8:124028/149‑1 (MQ=255)
                                                                                               tCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGAt                                                                                       >  1:67957/1‑113 (MQ=255)
                                                                                               tCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGAt                                                                                       <  2:67957/113‑1 (MQ=255)
                                                                                                 cTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGa                                                 >  7:111928/1‑149 (MQ=255)
                                                                                                  tGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAg                                                >  7:6940/1‑149 (MQ=255)
                                                                                                         ggTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAGCCGTGACGGCCGAATGCACCGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGGACCTGGGCTTCACGCCTGTTCGAGTTCCAcg                                         >  6:161973/1‑149 (MQ=255)
                                                                                                         ggTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCAcg                                         >  6:168853/1‑149 (MQ=255)
                                                                                                                cTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGt                                  >  2:92172/1‑149 (MQ=255)
                                                                                                                   ggTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTc                                 <  2:165720/147‑1 (MQ=255)
                                                                                                                    gTGCGACGTTCAACATCAGCCTTCACACCACGTCTGAAGGCGTGGTAACCGGGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCGGTTCGAGTTCCACGGCGGGGTCTAc                              <  8:268286/149‑1 (MQ=255)
                                                                                                                        gACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAAc                           >  6:388313/1‑148 (MQ=255)
                                                                                                                                aCATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGAc                  >  6:252473/1‑149 (MQ=255)
                                                                                                                                aCATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGAc                  >  8:286914/1‑149 (MQ=255)
                                                                                                                                   tCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTa                >  3:385947/1‑148 (MQ=255)
                                                                                                                                        aTTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGa          >  7:39681/1‑149 (MQ=255)
                                                                                                                                         ttCACACCACGTCTGAAGGCGTGGTAACCGTGACCGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCACGTCAGTACCTGGGCTTCACGCCTCTTCGAGTTCCTCGGGGGGGTCTGCAACCAGGTCGACTACAACGAc         >  8:71305/1‑149 (MQ=255)
                                                                                                                                          ggttaccacGCCTTCAGACGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACg        >  7:103027/5‑149 (MQ=255)
                                                                                                                                           cacaCCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACg        <  5:388313/148‑1 (MQ=255)
                                                                                                                                            acacCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCgggg                                   <  6:189744/120‑1 (MQ=255)
                                                                                                                                            acacCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCgggg                                   >  5:189744/1‑120 (MQ=255)
                                                                                                                                            acacCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGc       <  4:385947/148‑1 (MQ=255)
                                                                                                                                              accacGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAAc    >  5:238993/1‑149 (MQ=255)
                                                                                                                                                cacGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAcccc  >  4:223207/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCC  >  NC_002947/4301162‑4301454

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