Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,489,781 T→A 8.5% M88K (ATG→AAG) ‡ PP_4825 → MarC family membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,489,7810TA 8.5% 228.9 / 8.7 82M88K (ATG→AAG) ‡PP_4825MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base A (3/4);  total (40/42)
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.52e-01

AACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAAAACTG  >  NC_002947/5489634‑5489929
                                                                                                                                                    |                                                                                                                                                    
aaCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATg                                                                                                                                                     >  5:236622/1‑149 (MQ=255)
 aCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATgg                                                                                                                                                    <  3:290305/149‑1 (MQ=255)
  cTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGc                                                                                                                                                   >  8:188875/1‑149 (MQ=255)
       cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                              >  4:234733/1‑149 (MQ=255)
       cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                              >  5:85013/1‑149 (MQ=255)
       cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGACCTGGGCATGGCCCAgg                                                                                                                                              >  8:2545/1‑149 (MQ=255)
        tGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                             >  2:320281/1‑149 (MQ=255)
          gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAggg                                                                                                                                             <  6:15922/147‑1 (MQ=255)
          gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAGGGc                                                                                                                                            <  7:102652/148‑1 (MQ=255)
          gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAGGGc                                                                                                                                            <  8:191658/148‑1 (MQ=255)
          gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGATCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                             <  1:26070/147‑1 (MQ=255)
              gTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGcc                                                                                                                                       >  3:407952/1‑149 (MQ=255)
                  cgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggcg                                                                                                                                   <  1:328346/149‑1 (MQ=255)
                  cgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggcg                                                                                                                                   <  1:329582/149‑1 (MQ=255)
                  cgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTGCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggcg                                                                                                                                   <  6:377828/149‑1 (MQ=255)
                            cTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAgg                                                                                                                            >  4:192031/1‑146 (MQ=255)
                            cTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAg                                                                                                                         <  8:35470/149‑1 (MQ=255)
                               ggTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGAc                                                                                                                       <  2:258715/148‑1 (MQ=255)
                                gTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACaa                                                                                                                     >  2:327579/1‑149 (MQ=255)
                                     cAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGt                                                                                                                  >  1:315189/1‑148 (MQ=255)
                                          tgtTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcagcag                                                                                                           <  6:224981/149‑1 (MQ=255)
                                            ttGCTGTACCTGTTCGGGCGGGAGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAgg                                                                                                          <  3:324740/148‑1 (MQ=255)
                                              gCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAt                                                                                                        <  2:315189/148‑1 (MQ=255)
                                               cTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACAGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAtgt                                                                                                      >  3:254704/1‑149 (MQ=255)
                                                 gTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGa                                                                                                    >  8:328423/1‑149 (MQ=255)
                                                  tACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGAc                                                                                                   <  6:445167/149‑1 (MQ=255)
                                                      tGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCAGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATc                                                                                               >  5:76109/1‑149 (MQ=255)
                                                       gTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCg                                                                                              >  3:410619/1‑149 (MQ=255)
                                                          cgggcgggCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTa                                                                                            <  3:192031/148‑1 (MQ=255)
                                                          cgggcgggCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTAc                                                                                           >  1:70293/1‑149 (MQ=255)
                                                           gggcgggcGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTAc                                                                                           <  2:54659/148‑1 (MQ=255)
                                                            ggcgggcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCg                                                                                         >  4:323935/1‑149 (MQ=255)
                                                             gcgggctTTGTTCGGTATTTTGGGCATCACGGCCTATGCGTTCGGCATCGGTGCTGGCAGGGTGCTGTTCATTTCGGCGCTGGGCTTGGCCCTGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGc                                                                                        <  1:161209/149‑1 (MQ=255)
                                                                ggcttttttCGTAATTTTCGGCATCACCGCCGTTGCCTTCCGCATCGGTGGTGGCAGCGTGGTGTTCATTTCGGCCCTGGGCATGCCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGa                                                                                     <  2:167789/142‑1 (MQ=255)
                                                                 gcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGAGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGAc                                                                                    <  2:8689/149‑1 (MQ=255)
                                                                   gATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCa                                                                                  >  5:387184/1‑149 (MQ=255)
                                                                         cGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGc                                                                            >  7:84029/1‑149 (MQ=255)
                                                                          ggAATTTTCGGCATCACCGCCGATGCCTTCCGCATGGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt                                                                           <  4:407952/149‑1 (MQ=255)
                                                                          ggAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGCCCATCGTACCGCTGACCATTcagat                                                                           >  8:309306/1‑145 (MQ=255)
                                                                                    ggCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGTCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTg                                                                   <  7:2545/148‑1 (MQ=255)
                                                                                    ggCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTgg                                                                  <  3:234733/149‑1 (MQ=255)
                                                                                     gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGc                                                                 >  3:75929/1‑149 (MQ=255)
                                                                                      cATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGcc                                                                >  6:354273/1‑149 (MQ=255)
                                                                                      cATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGcc                                                                <  6:399340/149‑1 (MQ=255)
                                                                                       aTCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGccc                                                               >  4:163797/1‑149 (MQ=255)
                                                                                       aTCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGccc                                                               <  7:188875/149‑1 (MQ=255)
                                                                                          aCCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCgg                                                            <  2:70293/149‑1 (MQ=255)
                                                                                          aCCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCgg                                                            >  8:408917/1‑149 (MQ=255)
                                                                                           ccgccgATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCgg                                                            >  5:364272/1‑148 (MQ=255)
                                                                                              ccgATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCAcc                                                        <  7:309306/149‑1 (MQ=255)
                                                                                               cgATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCAcc                                                        >  2:99575/1‑148 (MQ=255)
                                                                                                gATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCAt                                                      <  7:222604/149‑1 (MQ=255)
                                                                                                gATGCCTTCCGCATCGGTGCTGGCACCCTGCTGTTCATTTCGGCCCTGGCCATCGCCCAGGGCAAGCCGCCGGTAAAGGCAGAAAACGTTCAGCAGGATGTGACCATCGTACCGCGGCCCATTCCGCTGACCGTTGGCCCCGGCACCAt                                                      >  7:271419/1‑149 (MQ=255)
                                                                                                      ttCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGCCCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGc                                                >  7:367810/1‑149 (MQ=255)
                                                                                                       tCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGcc                                               >  4:416342/1‑149 (MQ=255)
                                                                                                        ccGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGcc                                               <  1:99575/148‑1 (MQ=255)
                                                                                                        cagCATCGGTGCTGGCAGCGTGCTGTTCATTGCGGCCCTGTGCATGGCCCAGGGACAGCCGGCTGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGcc                                               <  4:242142/146‑1 (MQ=255)
                                                                                                          gCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGAATCGACCAGGGCAACCCGGCGGTAAAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATAACCCTCACCGAAGGCCCCGGCACCAACCGTGCCctg                                            >  2:374652/1‑149 (MQ=255)
                                                                                                             tCGGTGCTGGAAGGGTGCTTTTCTTTTCGGCCCTGGGCATTGCCCAGGGCAAGCCGGCTTTAGAGGCAGACAACGTTCATCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTTGCCCCGGCACCATCGGTGCCctgct                                          <  1:374652/148‑1 (MQ=255)
                                                                                                                      gtaGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgggg                                <  5:397763/147‑1 (MQ=255)
                                                                                                                      gCAGCGTGCTGTTCATTTCGGCCGTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgggg                                <  6:85013/149‑1 (MQ=255)
                                                                                                                      gCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgggg                                <  3:151156/149‑1 (MQ=255)
                                                                                                                      gCAGCGTGCTGTTCATTTCGGCCCGGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgggg                                <  6:236622/149‑1 (MQ=255)
                                                                                                                        aGCGTGCTGTTCATTTCGGCCCTGGGCAGGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGt                              <  2:2781/149‑1 (MQ=255)
                                                                                                                             gCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATg                                   >  8:266782/1‑139 (MQ=255)
                                                                                                                             gCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATg                                   <  7:266782/139‑1 (MQ=255)
                                                                                                                               tGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGagccag                       >  6:203846/1‑149 (MQ=255)
                                                                                                                               tGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGagccag                       >  3:305494/1‑149 (MQ=255)
                                                                                                                               tGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGagccag                       >  5:159992/1‑149 (MQ=255)
                                                                                                                                       tCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACt                <  1:338515/148‑1 (MQ=255)
                                                                                                                                        ttgCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGa                            <  6:122347/133‑1 (MQ=255)
                                                                                                                                        ttgCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGa                            >  5:122347/3‑135 (MQ=255)
                                                                                                                                        ttgCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTgg              >  4:261347/3‑149 (MQ=255)
                                                                                                                                         ggCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCTCCATCGGTGCCCTGCTGGTGATGGGGGTGAGCGAGCCGCACTggg             <  8:84029/149‑1 (MQ=255)
                                                                                                                                          gCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGa            >  4:168016/1‑149 (MQ=255)
                                                                                                                                           cccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctgct                                          <  2:50302/118‑1 (MQ=255)
                                                                                                                                           cccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctgct                                          >  1:50302/1‑118 (MQ=255)
                                                                                                                                           cccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGAt           <  7:328423/149‑1 (MQ=255)
                                                                                                                                           cccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGAt           >  6:137545/1‑149 (MQ=255)
                                                                                                                                            ccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGAt           >  5:342641/1‑148 (MQ=255)
                                                                                                                                                  cATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGTTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAAAAc    <  1:320281/149‑1 (MQ=255)
                                                                                                                                                  cATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGaaaaa     >  6:55824/1‑148 (MQ=255)
                                                                                                                                                    tGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAAAActg  <  4:410619/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAAAACTG  >  NC_002947/5489634‑5489929

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: