Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,711,818 T→A 6.3% L17Q (CTG→CAG) ‡ PP_4171 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,711,8180TA 6.3% 239.5 / 5.0 79L17Q (CTG→CAG) ‡PP_4171hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/32);  new base A (3/2);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GCCACGCGGTACTTTTCAAGACGTACGCCCGACACCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGAACTGATTCGCGCCAGCG  >  NC_002947/4711670‑4711964
                                                                                                                                                    |                                                                                                                                                   
gCCACGCGGTACTTTTCAAGACGTACGCCCGACACCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCt                                                                                                                                                     <  3:347427/149‑1 (MQ=255)
         tACTTTTCAAGACGTACGCCCGACACCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTcggc                                                                                                                                              >  6:202520/1‑147 (MQ=255)
         tACTTTTCAAGACGTACGCCCGACACCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTcggc                                                                                                                                              <  5:202520/147‑1 (MQ=255)
            ttttCAAGACGTACGCCCGACACCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCa                                                                                                                                         <  8:22998/149‑1 (MQ=255)
                  aGACGTACGCCCGACACCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCGAGCAGGGCGGTc                                                                                                                                          >  6:292882/1‑142 (MQ=255)
                  aGACGTACGCCCGACACCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCGAGCAGGGCGGTc                                                                                                                                          <  5:292882/142‑1 (MQ=255)
                   gACGTACGCCCGACACCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTg                                                                                                                                  >  1:346933/1‑149 (MQ=255)
                            ccGACACCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCgacgac                                                                                                                        >  4:196792/1‑149 (MQ=255)
                               acacCCACTAGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTg                                                                                                                     >  6:298386/1‑149 (MQ=255)
                                 acCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGTCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTcggcgg                                                                                                                                            >  6:139708/1‑125 (MQ=255)
                                 acCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGTCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTcggcgg                                                                                                                                            <  5:139708/125‑1 (MQ=255)
                                 acCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGt                                                                                                                    >  4:369245/1‑148 (MQ=255)
                                   ccACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCa                                                                                                                 <  6:59192/149‑1 (MQ=255)
                                    cACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAg                                                                                                                                        >  2:298371/1‑126 (MQ=255)
                                    cACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCAt                                                                                                                <  2:207954/149‑1 (MQ=255)
                                       tGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCGAGCAGGGCGGTCa                                                                                                                                         >  4:247331/1‑122 (MQ=255)
                                       tGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCGAGCAGGGCGGTCa                                                                                                                                         <  3:247331/122‑1 (MQ=255)
                                           aaaaCCCCCACGGTGAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGacaac                                                                                                         <  2:7969/149‑1 (MQ=255)
                                           aaaaCCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGacaa                                                                                                          >  5:203107/1‑148 (MQ=255)
                                           aaaaCCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGacaa                                                                                                          >  4:58028/1‑148 (MQ=255)
                                           aaaaCCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGacaa                                                                                                          >  8:99689/1‑148 (MQ=255)
                                              acacccACGGTCATTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCTGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACAcc                                                                                                      <  6:241190/146‑1 (MQ=255)
                                               cccccACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACAcc                                                                                                      >  1:123133/1‑148 (MQ=255)
                                                 cccACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGAc                                                                                                   <  5:196908/149‑1 (MQ=255)
                                                       gTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCca                                                                                             <  3:195609/149‑1 (MQ=255)
                                                            ttACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGt                                                                                        >  5:166227/1‑149 (MQ=255)
                                                              aCCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGtt                                                                                       <  3:196792/148‑1 (MQ=255)
                                                               ccATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTc                                                                                      <  2:289500/148‑1 (MQ=255)
                                                                 aTGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGa                                                                                   <  6:393854/149‑1 (MQ=255)
                                                                 aTGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGa                                                                                   >  8:69802/1‑149 (MQ=255)
                                                                 aTGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCGGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGa                                                                                   >  2:254599/1‑149 (MQ=255)
                                                                    ccTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCg                                                                                <  5:298386/149‑1 (MQ=255)
                                                                      tGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGc                                                                               >  6:436523/1‑148 (MQ=255)
                                                                         ttCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGCCCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCcaca                                                                           >  2:278563/1‑149 (MQ=255)
                                                                         ttCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCcaca                                                                           >  1:427570/1‑149 (MQ=255)
                                                                           cAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGc                                                                               >  7:361292/1‑143 (MQ=255)
                                                                           cAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGc                                                                               <  8:361292/143‑1 (MQ=255)
                                                                           cAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCAACGACAACACCGACCTGCCACAGGTTCTGACCGCCACACg                                                                         >  1:200013/1‑149 (MQ=255)
                                                                            aaTCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCacaccc                                                                        <  7:99689/149‑3 (MQ=255)
                                                                              tCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGGGGTCAGGCCGTGGCCGACGAAGAAGTTAATCCACAACACCGACAGGCCACAGGTTATGCCCCCCCCACGcc                                                                       >  1:142466/1‑149 (MQ=255)
                                                                              tCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCt                                                                      >  4:214831/1‑149 (MQ=255)
                                                                              tCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCt                                                                      >  8:75864/1‑149 (MQ=255)
                                                                                cTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTg                                                                     >  6:169055/1‑148 (MQ=255)
                                                                                cTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTg                                                                     <  5:169055/148‑1 (MQ=255)
                                                                                     aCAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGCCAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACaa                                                               <  1:315308/149‑1 (MQ=255)
                                                                                     aCAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACaa                                                               <  1:204820/149‑1 (MQ=255)
                                                                                       aaaGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATc                                                             <  1:118318/149‑1 (MQ=255)
                                                                                         aGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATccc                                                           >  8:49680/1‑149 (MQ=255)
                                                                                          ggCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATcccc                                                          >  2:370203/1‑149 (MQ=255)
                                                                                             cGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCccgc                                                       >  1:9444/1‑149 (MQ=255)
                                                                                                  cGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTg                                                  <  3:403702/149‑1 (MQ=255)
                                                                                                    tgtTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTc                                                <  2:346933/149‑1 (MQ=255)
                                                                                                     gtTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTcc                                               >  5:75938/1‑149 (MQ=255)
                                                                                                         ccTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGCCCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCggg                                           >  4:408926/1‑149 (MQ=255)
                                                                                                          cTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCggg                                           >  8:51779/1‑148 (MQ=255)
                                                                                                           tGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCggg                                           >  1:351824/1‑147 (MQ=255)
                                                                                                             gCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGgcag                                        <  5:436523/148‑1 (MQ=255)
                                                                                                             gCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGgca                                         <  6:203107/147‑1 (MQ=255)
                                                                                                                   ccccGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGAcc                                 >  6:386799/1‑149 (MQ=255)
                                                                                                                     ccGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGt                               <  7:51779/149‑1 (MQ=255)
                                                                                                                       gTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGCCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTAc                             <  7:319260/149‑1 (MQ=255)
                                                                                                                        tACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACt                            <  2:9444/149‑1 (MQ=255)
                                                                                                                         aCCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTc                           >  6:65547/1‑149 (MQ=255)
                                                                                                                            cgcgCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGa                         >  6:409390/1‑148 (MQ=255)
                                                                                                                             gcgcCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACa                       <  1:370203/149‑1 (MQ=255)
                                                                                                                                ccTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGc                     >  4:169903/1‑148 (MQ=255)
                                                                                                                                 ctgctgTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCg                    >  7:247494/1‑148 (MQ=255)
                                                                                                                                 ctgctgTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCg                    >  8:351107/1‑148 (MQ=255)
                                                                                                                                  tgctgTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGaa                  >  4:79714/1‑149 (MQ=255)
                                                                                                                                    cacggcctGACCGCCGAGCAGGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGAACt                >  4:217317/8‑149 (MQ=255)
                                                                                                                                      gTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGAACTGa              >  8:406577/1‑149 (MQ=255)
                                                                                                                                       ttATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGAACTg               <  3:214831/147‑1 (MQ=255)
                                                                                                                                         aTTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGAACTGAtt            <  3:369245/148‑1 (MQ=255)
                                                                                                                                          ttGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACACGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGAACTGATTc           <  8:247494/148‑1 (MQ=255)
                                                                                                                                          ttGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACAGACCTCCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCTCTGTCCAGCGCAGCATGACCGTACTCGACAGCGAACTGATTcg          >  2:54720/1‑149 (MQ=255)
                                                                                                                                            gACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGCCCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGACCTGATTcgcg        >  7:167554/1‑149 (MQ=255)
                                                                                                                                            gACCGCCCTGCTCGGAGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGAACTGATTcgcg        >  4:337984/1‑149 (MQ=255)
                                                                                                                                             aCCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCTGCATGACCGTACTCGACAGCGAACTGATTCgcgc       <  7:406577/149‑1 (MQ=255)
                                                                                                                                             aCCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGAACTGATTCgcgc       >  7:158164/1‑149 (MQ=255)
                                                                                                                                                  ccTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGAACTGATTCGCGCCAGCg  <  1:254599/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCACGCGGTACTTTTCAAGACGTACGCCCGACACCCACTGGCAAAACCCCCACGGTCAGTTACCATGCCTGTTTCAATCCTGAAACAAAGGCCGCTCCGTGTTCCCTGGCACGCCCCCGTACCCGCGCCTGCTGTTATTGACCGCCCTGCTCGGCGGTCAGGCCGTGGCCGACGACCTGTTCATCGACAACACCGACCTGCCACAGGTTCTGACCGCCACACGCCTGAAACAATCCCCCGCCGCTGTCCCGGGCAGCATGACCGTACTCGACAGCGAACTGATTCGCGCCAGCG  >  NC_002947/4711670‑4711964

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: