Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,802,077 T→G 5.1% D90A (GAC→GCC)  PP_4230 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,802,0770TG 5.1% 233.8 / 2.4 74D90A (GAC→GCC) PP_4230hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/30);  new base G (2/2);  total (42/32)
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.00e+00

CGAACAACACCGGATGGCCCTGTCGTCCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCTTCACGCAACACCACCCGCACATCATCGAA  >  NC_002947/4801931‑4802219
                                                                                                                                                  |                                                                                                                                              
cGAACAACACCGGATGGCCCTGTCGTCCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTca                                                                                                                                              <  8:190718/149‑1 (MQ=255)
   acaacaCCGGATGGCCCTGTCGTCCAGCATNACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTcagcc                                                                                                                                           >  1:35185/1‑149 (MQ=255)
   acaacaCCGGATGGCCCTGTCGTCCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTcagc                                                                                                                                            <  4:298312/148‑1 (MQ=255)
        aaccgATTGCCCTGTATTCCAGCAACACAAGGCAGAAATATCGTTTACCCAAAGGCCTCGCACGAAAACAAATACAAAGACGCGGACAACACCCAGGGCATATCACCCAACAGAATCTCTACTGCCTGCGCCTCGCTGTcacccagcc                                                                                                                                       >  8:163614/5‑147 (MQ=255)
          cGGATGGCCCTGTCGTCCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAg                                                                                                                                    <  3:239030/149‑1 (MQ=255)
                gCCCTGTCGTCCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCTTCGCTGGCTGACAGCGAGGCgcag                                                                                                                               <  6:364252/148‑4 (MQ=255)
                 cccTGTCGTCCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAg                                                                                                                                    >  4:349332/1‑142 (MQ=255)
                 cccTGTCGTCCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAg                                                                                                                                    <  3:349332/142‑1 (MQ=255)
                 cccTGTCGTCCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCAcc                                                                                                                             >  1:191541/1‑149 (MQ=255)
                  ccTGTCGTCCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCg                                                                                                                            >  4:259308/1‑149 (MQ=255)
                       cgtcCAGCATGACAAGGCAGGACTATCGTTGACCAACTGGCCTCGCGCGCAAGCAAATCCAGGGCCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGgccgc                                                                                                                       >  6:323654/1‑149 (MQ=255)
                         tcCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGgccgc                                                                                                                       >  3:163629/1‑147 (MQ=255)
                             gCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGc                                                                                                                  >  3:95417/1‑148 (MQ=255)
                              cATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGCCCCCCCCCAGGGCATATCACCCAACAGCATCGCTACGGCCTGCCCCTACCAGGACGCCAGCGAAGCCGCACCGGCCGCCAGGcg                                                                                                                 >  2:295333/1‑147 (MQ=255)
                                 gACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATgg                                                                                                             >  5:116250/1‑149 (MQ=255)
                                    aaGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGcc                                                                                                           >  5:384675/1‑148 (MQ=255)
                                       gCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATg                                                                                                       >  7:234156/1‑149 (MQ=255)
                                         aGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATg                                                                                                       >  3:185845/1‑147 (MQ=255)
                                         aGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATg                                                                                                       <  4:185845/147‑1 (MQ=255)
                                                tCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAAtgt                                                                                               >  5:327347/1‑148 (MQ=255)
                                                   ttGACCCACTGGCCTCGCGCGCAAGCAAATCCTGGGTCGCGGACGCCACCCATGGCATATCACCAATCAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGaaa                                                                                           <  7:163614/149‑1 (MQ=255)
                                                   ttGACACACTGGCCTCGCGCGCAAGCAAATCCAGGTTCGCGTACGCCACCCCGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGACAGCGAACCGGCACCGGCCGCCAGCCTATGTCCCATTCAAAATGTGaaa                                                                                           >  8:376802/1‑149 (MQ=255)
                                                     gACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAAtt                                                                                         >  7:12744/1‑149 (MQ=255)
                                                      aCCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGccc                                                                                                          <  6:218109/131‑1 (MQ=255)
                                                      aCCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGccc                                                                                                          >  5:218109/1‑131 (MQ=255)
                                                      aCCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAAttt                                                                                        <  4:163629/149‑1 (MQ=255)
                                                       cccACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTg                                                                                       <  2:191541/149‑1 (MQ=255)
                                                         cACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGa                                                                                     <  6:327347/149‑1 (MQ=255)
                                                         cACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGa                                                                                     <  1:295333/149‑1 (MQ=255)
                                                            tGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTg                                                                                  >  4:69530/1‑149 (MQ=255)
                                                            tGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTg                                                                                  >  4:75058/1‑149 (MQ=255)
                                                              gCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTgg                                                                                 <  3:259308/148‑1 (MQ=255)
                                                                     gcgcAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCaa                                                                                                  <  5:86006/124‑1 (MQ=255)
                                                                     gcgcAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCaa                                                                                                  >  6:86006/1‑124 (MQ=255)
                                                                           gcaaATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGAc                                                                                    >  6:208709/1‑132 (MQ=255)
                                                                           gcaaATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGAc                                                                                    <  5:208709/132‑1 (MQ=255)
                                                                           gcaaATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGG‑CCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCAcc                                                                  >  1:109058/1‑149 (MQ=255)
                                                                                tCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACcgcg                                                               >  1:147784/1‑148 (MQ=255)
                                                                                ccctgggTGGCGTCCGCGACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGga                                                                                >  5:386659/5‑131 (MQ=255)
                                                                                ccctgggTGGCGTCCGCGACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGga                                                                                <  6:386659/127‑1 (MQ=255)
                                                                                  cAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCa                                                             <  6:384675/148‑1 (MQ=255)
                                                                                  cAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAg                                                            <  2:147784/149‑1 (MQ=255)
                                                                                  cAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAg                                                            >  3:185886/1‑149 (MQ=255)
                                                                                   aGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCGCGGTGCAgg                                                           >  6:111803/1‑149 (MQ=255)
                                                                                       tCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCAccc                                                                        >  1:353857/1‑132 (MQ=255)
                                                                                       tCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCAccc                                                                        <  2:353857/132‑1 (MQ=255)
                                                                                       tCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAAc                                                       <  8:12744/149‑1 (MQ=255)
                                                                                         gcgGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCGCGGTGCAgggtggtc                                                     <  5:111803/149‑7 (MQ=255)
                                                                                               gCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCAAGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCCCCGCGCAGGCAACCCCCGGTc                                               >  7:289312/1‑149 (MQ=255)
                                                                                                 cACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGATACCACCCTGCACCGCCCAGGCAACCCCACGTCCa                                             >  6:303173/1‑149 (MQ=255)
                                                                                                 cACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCa                                             <  7:372497/149‑1 (MQ=255)
                                                                                                    ccAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGCCCACCCTGCACCGCGCAGGCAACCACAGGGCCCCGc                                          >  2:359887/1‑149 (MQ=255)
                                                                                                     cAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAg                                                            >  4:93717/1‑130 (MQ=255)
                                                                                                     cAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAg                                                            <  3:93717/130‑1 (MQ=255)
                                                                                                             atCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTccc                                 >  7:33037/1‑149 (MQ=255)
                                                                                                                aCCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCtt                              >  8:282651/1‑149 (MQ=255)
                                                                                                                aCCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCtt                              >  4:1874/1‑149 (MQ=255)
                                                                                                                  ccAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCATGCCTCGCTCCCCTTCa                            <  8:234156/149‑1 (MQ=255)
                                                                                                                  ccAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCTTc                             >  4:44660/1‑148 (MQ=255)
                                                                                                                   cAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACcgc                                                                >  8:137013/1‑112 (MQ=255)
                                                                                                                   cAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACcgc                                                                <  7:137013/112‑1 (MQ=255)
                                                                                                                    aaCAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCTTCACg                          >  7:205660/1‑149 (MQ=255)
                                                                                                                        gAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCTTCACGCAac                      >  5:18739/1‑149 (MQ=255)
                                                                                                                                 ctgcctgcGCCTCGCTGTCAGCAAGCTAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCAAAGGTCCACGCCTCGCTCCCCTTCTCGCAACAACAccc               <  7:376802/147‑1 (MQ=255)
                                                                                                                                    cctgcGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGaccacc                                                                         <  1:105148/86‑1 (MQ=255)
                                                                                                                                    cctgcGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGaccacc                                                                         >  2:105148/1‑86 (MQ=255)
                                                                                                                                      tgcGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCtt                              <  4:55447/127‑1 (MQ=255)
                                                                                                                                      tgcGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCtt                              >  3:55447/1‑127 (MQ=255)
                                                                                                                                        cgcCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCTTCACGCAACACCACCCGCAcatc        <  3:44660/147‑1 (MQ=255)
                                                                                                                                         gcttCGCTGGCTGACAGCGAGGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCt                                        <  4:295677/111‑1 (MQ=255)
                                                                                                                                         gcttCGCTGGCTGACAGCGAGGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCt                                        >  3:295677/4‑114 (MQ=255)
                                                                                                                                            tCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCTTCACGCAACACCACCCGCACATCATCGaa  >  6:334327/1‑149 (MQ=255)
                                                                                                                                              gCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCTTCACGCAACACCACCCGCACATCATCGaa  <  4:185886/147‑1 (MQ=255)
                                                                                                                                              gCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCTTCACGCAACACCACCCGCACATCATCGaa  <  6:18739/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CGAACAACACCGGATGGCCCTGTCGTCCAGCATGACAAGGCAGGACTATCGTTGACCCACTGGCCTCGCGCGCAAGCAAATCCAGGGTCGCGGACGCCACCCAGGGCATATCACCCAACAGAATCGCTACTGCCTGCGCCTCGCTGTCAGCCAGCGAAGCGGCACCGGCCGCCAGGCTATGGCCCATGCCCAATGTGAAATTTGGACTGGAGACCACCCTGCACCGCGCAGGCAACCCCAGGTCCACGCCTCGCTCCCCTTCACGCAACACCACCCGCACATCATCGAA  >  NC_002947/4801931‑4802219

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