Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,314 A→C 10.5% E160A (GAA→GCA) ‡ PP_3774 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,3140AC10.5% 187.7 / 13.8 67E160A (GAA→GCA) ‡PP_3774hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base C (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 = 9.14e-01

GGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGA  >  NC_002947/4301174‑4301461
                                                                                                                                            |                                                                                                                                                   
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)
                                                                                                                                            aaGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGa  <  8:6940/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGA  >  NC_002947/4301174‑4301461

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