Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,317 G→A 10.1% G161D (GGC→GAC)  PP_3774 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,3170GA10.1% 188.7 / 8.3 69G161D (GGC→GAC) PP_3774hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base A (3/4);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

GGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAG  >  NC_002947/4301174‑4301462
                                                                                                                                               |                                                                                                                                                 
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                                                                                               <  2:67957/113‑1 (MQ=255)
                                                                                   tCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGAt                                                                                               >  1:67957/1‑113 (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)
                                                                                                                                            atggCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGa   <  7:71305/146‑1 (MQ=255)
                                                                                                                                            aaGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGa   <  8:6940/148‑1 (MQ=255)
                                                                                                                                             aGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAg  >  4:133315/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAG  >  NC_002947/4301174‑4301462

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: ‑