Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,416,131:1 +G 8.1% intergenic (+29/+90) PP_3901 → / ← PP_3902 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,416,1311.G 8.1% 480.8 / 33.9 123intergenic (+29/+90)PP_3901/PP_3902hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (58/55);  new base G (4/6);  total (62/61)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 124 total aligned reads displayed.

GAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGAAAGTCG  >  NC_002947/4415983‑4416277
                                                                                                                                                      |                                                                                                                                                  
   tCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGCAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCgg                                                                                                                                                >  2:360375/1‑150 (MQ=255)
   tCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCgg                                                                                                                                                <  1:793465/150‑1 (MQ=255)
     ccAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGc                                                                                                                                              >  2:910304/1‑150 (MQ=255)
     ccAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGc                                                                                                                                              >  1:1535266/1‑150 (MQ=255)
       aGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaa                                                                                                                                            <  2:1506712/150‑3 (MQ=255)
       aGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaa                                                                                                                                            <  2:157067/150‑3 (MQ=255)
          aaGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaaggt                                                                                                                                         <  1:1650504/150‑6 (MQ=255)
             aCTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCg                                                                                                                                                 >  2:2179500/1‑139 (MQ=255)
             aCTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCg                                                                                                                                                 <  1:2179500/139‑1 (MQ=255)
                ggggCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGCTGGCAGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTgcgc                                                                                                                                   <  1:910304/150‑1 (MQ=255)
                 gggCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTGGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGcc                                                                                                                                  >  2:157542/1‑150 (MQ=255)
                       gcTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGACGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGaa                                                                                                                            >  1:2140500/1‑150 (MQ=255)
                         tcgtcgAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGaaaa                                                                                                                          <  1:1762597/150‑1 (MQ=255)
                         tcgtcgAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaa                                                                                                                                            <  2:876564/132‑3 (MQ=255)
                           gtcgAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGt                                                                                                                        <  1:1135718/150‑1 (MQ=255)
                            tcgAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTg                                                                                                                       >  1:2048821/1‑150 (MQ=255)
                             ggAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTCGCGCCATAGAAAAGTGa                                                                                                                      >  2:1657035/2‑150 (MQ=255)
                                   cATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATg                                                                                                                >  2:1688014/1‑150 (MQ=255)
                                   cATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATg                                                                                                                <  2:102905/150‑1 (MQ=255)
                                    aTGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATATAAAAGTGAATTATgg                                                                                                               <  2:1015228/150‑1 (MQ=255)
                                    aTGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATgg                                                                                                               >  2:542325/1‑150 (MQ=255)
                                    aTGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGACGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATgg                                                                                                               <  2:2140500/150‑1 (MQ=255)
                                     gggAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATggg                                                                                                              <  2:550789/149‑1 (MQ=255)
                                      ggAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGt                                                                                                             <  2:1150226/150‑1 (MQ=255)
                                         aGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaag                                                                                                                                           <  2:1841175/117‑4 (MQ=255)
                                         aGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaag                                                                                                                                           >  1:1841175/1‑114 (MQ=255)
                                           gCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTaccac                                                                                                        <  1:2254884/150‑1 (MQ=255)
                                           gCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTAATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTaccac                                                                                                        <  1:614938/150‑1 (MQ=255)
                                            cACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATa                                                                                                                               >  1:1472598/1‑126 (MQ=255)
                                            cACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTaccacc                                                                                                       >  2:1829268/1‑150 (MQ=255)
                                             aCCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTAccacca                                                                                                      >  1:1670463/1‑150 (MQ=255)
                                             aCCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTAccacca                                                                                                      >  2:2144341/1‑150 (MQ=255)
                                             aCCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCATCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTAacacca                                                                                                      >  2:2054823/1‑150 (MQ=255)
                                             aCCGAGGGCGTAAGCCTGGCCTATGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTAccacca                                                                                                      >  1:1508628/1‑150 (MQ=255)
                                                gAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATaa                                                                                                   >  2:831682/1‑150 (MQ=255)
                                                gAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATaa                                                                                                   <  1:1993164/150‑1 (MQ=255)
                                                gAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATaa                                                                                                   <  2:1928404/150‑1 (MQ=255)
                                                       tAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTAt                                                                                            <  1:2144341/150‑1 (MQ=255)
                                                       tAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTAt                                                                                            <  2:912219/150‑1 (MQ=255)
                                                           ccTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGc                                                                                        >  1:650303/1‑150 (MQ=255)
                                                                  tttGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCa                                                                                 >  1:952826/1‑150 (MQ=255)
                                                                   ttGAAGATGATGGCGGCGTTTGTGTAACCTGGGAAAAGCCGCTAGAAGCCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAt                                                                                <  2:1084770/150‑1 (MQ=255)
                                                                   ttGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAt                                                                                <  1:1829268/150‑1 (MQ=255)
                                                                    tGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCAAATCGCTCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAtt                                                                               <  1:151542/150‑1 (MQ=255)
                                                                     gAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAt                                                                                <  2:1331017/148‑1 (MQ=255)
                                                                     gAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAAGTATGGGTACCACCATAACATTTATGAGCAAACCCAt                                                                                >  1:1331017/1‑148 (MQ=255)
                                                                      aaGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGa                                                                             <  2:2048821/150‑1 (MQ=255)
                                                                           gatgGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAGGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACaa                                                                        >  2:1395993/1‑150 (MQ=255)
                                                                                  gcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGCCTGATTTCAACATCCGCATCTTGCCCGCCTCGAGCGGGCTTTTGAGGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                                 <  1:542325/150‑1 (MQ=255)
                                                                                   cgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAg                                                                <  2:48414/150‑1 (MQ=255)
                                                                                   cgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAg                                                                <  2:727869/150‑1 (MQ=255)
                                                                                    gTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGc                                                               >  1:1003957/1‑150 (MQ=255)
                                                                                        gtgtGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatac                                                           >  1:265147/1‑150 (MQ=255)
                                                                                         tgtgACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAg                                                                                                                              >  2:1056731/1‑82 (MQ=255)
                                                                                         tgtgACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAg                                                                                                                              <  1:1056731/82‑1 (MQ=255)
                                                                                         tgtgACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTCCTCCATAGCAACATTAATTTTGAGTACCACCCTAACATTCACGAGCTACCCCTCTCACAAAAAACATTCGcatact                                                          >  2:1496642/1‑150 (MQ=255)
                                                                                               cTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactcatact                                                    >  2:105365/1‑150 (MQ=255)
                                                                                               cTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactcatact                                                    >  2:2101361/1‑150 (MQ=255)
                                                                                               cTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCAACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATgg                                                                                                               >  2:1893237/1‑91 (MQ=255)
                                                                                                tGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGaaaa                                                                                                                          <  1:1843490/79‑1 (MQ=255)
                                                                                                tGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATTGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCatactcatactc                                                   <  1:360375/150‑1 (MQ=255)
                                                                                                tGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCatactcatactc                                                   >  2:2254488/1‑150 (MQ=255)
                                                                                                 gggAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCg                                                  >  1:2392407/1‑150 (MQ=255)
                                                                                                   gAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGCCATGTCAATTATGCGTATTTTTACAACATGTATTAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCt                                                <  1:2101361/150‑1 (MQ=255)
                                                                                                        gCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAatat                                                                   <  2:1374619/126‑1 (MQ=255)
                                                                                                           gCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactc                                                         >  1:642156/1‑133 (MQ=255)
                                                                                                           gCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactc                                                         <  2:642156/133‑1 (MQ=255)
                                                                                                             tAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAata                                                                    <  2:1250169/120‑1 (MQ=255)
                                                                                                             tAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACAATAACATTTATGAGCAAACCCATTGACACAAAata                                                                    >  1:1250169/1‑120 (MQ=255)
                                                                                                               gcagaCTGATTTCGACAAAGACACCTTGCCCGCCTGGAGCGGGGTTTTTTTGCGCAATAGAAACGTGAATTATGGGTACCACCATAACATTTATTAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGAt                                    <  2:330340/148‑1 (MQ=255)
                                                                                                               gaagaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGAt                                    >  1:2337589/1‑150 (MQ=255)
                                                                                                               gaagaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGAt                                    <  2:1508628/150‑1 (MQ=255)
                                                                                                                     tGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCTGATAGacac                              >  2:1965220/1‑150 (MQ=255)
                                                                                                                     tGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGacac                              <  1:1463536/150‑1 (MQ=255)
                                                                                                                     tGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGacac                              <  1:831682/150‑1 (MQ=255)
                                                                                                                       aTTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCg                            <  2:2088259/150‑1 (MQ=255)
                                                                                                                        tttCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactcatac                                                     <  1:460965/124‑1 (MQ=255)
                                                                                                                            gACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATggg                                                                                                              <  1:1980048/63‑1 (MQ=255)
                                                                                                                            gACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATggg                                                                                                              >  2:1980048/1‑63 (MQ=255)
                                                                                                                            gACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTAt                                                                                            <  1:2194337/81‑1 (MQ=255)
                                                                                                                            gACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTAt                                                                                            >  2:2194337/1‑81 (MQ=255)
                                                                                                                                  cccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATgg                                                                                                               >  2:1619673/1‑56 (MQ=255)
                                                                                                                                  cccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATgg                                                                                                               <  1:1619673/56‑1 (MQ=255)
                                                                                                                                  cccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGcaacaa                 >  2:952113/1‑150 (MQ=255)
                                                                                                                                   cctactcgCCCGTCTTAAGCTGGCTTTGGTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACaaa                <  2:1715017/143‑1 (MQ=255)
                                                                                                                                   ccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcata                                                            >  2:2043007/1‑106 (MQ=255)
                                                                                                                                   ccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcata                                                            <  1:2043007/106‑1 (MQ=255)
                                                                                                                                     aCCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAgg              >  1:535714/1‑150 (MQ=255)
                                                                                                                                      ccTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATTGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGc             <  1:872792/150‑1 (MQ=255)
                                                                                                                                       aaagCCCGCTCGAGGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGc                          <  2:693589/133‑1 (MQ=255)
                                                                                                                                        ttGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAg           <  1:105365/150‑1 (MQ=255)
                                                                                                                                             cccgcTCGAGGCGGGCTTTTTTTACGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACCCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGa        >  1:666804/5‑150 (MQ=255)
                                                                                                                                             ccGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactc                                                         >  2:1310609/1‑100 (MQ=255)
                                                                                                                                             ccGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactc                                                         <  1:1310609/100‑1 (MQ=255)
                                                                                                                                             ccGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGaa       >  2:2215131/1‑150 (MQ=255)
                                                                                                                                             ccGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGaa       >  1:1025836/1‑150 (MQ=255)
                                                                                                                                               gCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATa                                   >  2:1524732/1‑120 (MQ=255)
                                                                                                                                                 cTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGcaa                    <  1:103973/133‑1 (MQ=255)
                                                                                                                                                 cTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGAAAGTc   >  1:793833/1‑150 (MQ=255)
                                                                                                                                                  tCGAGCGGGCTTTTTTTGCGCCATAGAAAATTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGAAAGTCg  >  2:469534/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGAAAGTCG  >  NC_002947/4415983‑4416277

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: