Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,313 G→T 10.5% E160* (GAA→TAA) ‡ PP_3774 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,3130GT10.5% 186.8 / 12.1 67E160* (GAA→TAA) ‡PP_3774hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base T (3/4);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

GGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCC  >  NC_002947/4301174‑4301454
                                                                                                                                           |                                                                                                                                             
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)
                                                                                                                                          tGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTc                     <  5:173351/124‑1 (MQ=255)
                                                                                                                                          tGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTc                     >  6:173351/1‑124 (MQ=255)
                                                                                                                                           |                                                                                                                                             
GGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCC  >  NC_002947/4301174‑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: ‑