Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,640,261 G→A 5.8% L147L (CTG→TTG)  PP_2310 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,640,2610GA 5.8% 203.9 / 2.5 69L147L (CTG→TTG) PP_2310methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (2/2);  total (35/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 = 6.33e-01

CCTCGTCGCCACTGGCTAGCGGTACGCGGGCCGACAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAGGATCCAGGCGCTCGGCAAT  >  NC_002947/2640114‑2640409
                                                                                                                                                   |                                                                                                                                                    
ccTCGTCGCCACTGGCTAGCGGTACGCGGGCCGACAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGca                                                                                                                                                    <  1:358277/150‑1 (MQ=255)
  tcgtcgCCACTGGCTAGCGGTACGCGGGCCGACAGGTCGCCGGTTTTTTCCAAGGCGATGGTCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGcaca                                                                                                                                                  >  2:1194629/1‑150 (MQ=255)
        ccACTGGCTAGCGGTACGCGGGCCGACAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACAcc                                                                                                                                            >  1:1023354/1‑150 (MQ=255)
        ccACTGGCTAGCGGTACGCGGGCCGACAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACAcc                                                                                                                                            >  1:1083081/1‑150 (MQ=255)
         cACTGGCTAGCGGTACGCGGGCCGACAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCAAGCATGGCCAGCATCAGGATCAGCACACACACcg                                                                                                                                           >  2:816635/1‑150 (MQ=255)
                           gggCCGACAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGc                                                                                                                                                     >  2:1361608/1‑122 (MQ=255)
                           gggCCGACAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGc                                                                                                                                                     <  1:1361608/122‑1 (MQ=255)
                           gggCCGACAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGAt                                                                                                                         <  1:271195/150‑1 (MQ=255)
                             gCCGACAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCg                                                                                                                       >  1:1375014/1‑150 (MQ=255)
                               cgccAGGTCGCCGGTTTTTTCCACGTCCAGCAGCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCgaa                                                                                                                     <  2:1375014/147‑1 (MQ=255)
                                  cAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCAGCGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCgaagaa                                                                                                                  <  2:86567/150‑1 (MQ=255)
                                        gCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGTGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGc                                                                                                                                  <  1:1068394/128‑1 (MQ=255)
                                        gCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGc                                                                                                                                  >  2:1068394/1‑128 (MQ=255)
                                        gCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGAATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAACACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTg                                                                                                            >  1:907670/1‑150 (MQ=255)
                                            gTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCACCGAGCATGGCCAGCCTGATCCTTAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGgcgc                                                                                                        >  2:749147/1‑150 (MQ=255)
                                             ttttttCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCa                                                                                                       <  1:889120/150‑1 (MQ=255)
                                                   ccACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCAAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTgg                                                                                                 <  1:816635/150‑1 (MQ=255)
                                                     aCGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCa                                                                                                       >  1:1351933/1‑142 (MQ=255)
                                                     aCGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCa                                                                                                       <  2:1351933/142‑1 (MQ=255)
                                                            gcaTCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTggg                                                                                        <  2:40451/150‑1 (MQ=255)
                                                            gcaTCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACCggg                                                                                        >  1:1208883/1‑150 (MQ=255)
                                                               tCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAagaaga                                                                                                                >  1:1000491/1‑123 (MQ=255)
                                                               tCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAagaaga                                                                                                                <  2:1000491/123‑1 (MQ=255)
                                                                  cGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGca                                                                                  <  1:516660/150‑1 (MQ=255)
                                                                   gTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACCCCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTCGGCGACTGCGcccccac                                                                                 >  2:1110216/1‑145 (MQ=255)
                                                                    tCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGcaca                                                                                <  1:1194629/150‑1 (MQ=255)
                                                                      cTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGcg                                                                                                                                                    >  1:799725/1‑79 (MQ=255)
                                                                      cTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGca                                                                                                                                                    <  2:799725/80‑1 (MQ=255)
                                                                      cTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGc                                                                              <  2:677184/150‑1 (MQ=255)
                                                                       ttGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGcc                                                                             >  2:322426/1‑150 (MQ=255)
                                                                         gAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACATGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCaa                                                                           >  1:632506/1‑150 (MQ=255)
                                                                           gggTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAAct                                                                         >  2:440909/1‑150 (MQ=255)
                                                                             gtgtTGAGTTGGCTGAGCAGGAAGCGGAGCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAActgc                                                                       >  1:749541/1‑150 (MQ=255)
                                                                               gtTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGc                                                                                       <  1:632831/132‑1 (MQ=255)
                                                                               gtTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGc                                                                                       >  2:632831/1‑132 (MQ=255)
                                                                                ttGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTgg                                                                    <  2:1208883/150‑1 (MQ=255)
                                                                                ttGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAACACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTgg                                                                                                 <  2:806920/121‑1 (MQ=255)
                                                                                ttGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAACACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTgg                                                                                                 >  1:806920/1‑121 (MQ=255)
                                                                                          tGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCgaa                                                                                                                     >  2:1176290/1‑91 (MQ=255)
                                                                                          tGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCgaa                                                                                                                     <  1:1176290/91‑1 (MQ=255)
                                                                                               aGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAAt                                                     >  2:1159027/1‑150 (MQ=255)
                                                                                                 gAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATAc                                                   >  1:930783/1‑150 (MQ=255)
                                                                                                            gcagTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGc                                                                                       <  2:1053393/103‑1 (MQ=255)
                                                                                                            gcagTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGc                                                                                       >  1:1053393/1‑103 (MQ=255)
                                                                                                                 tGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGTTGAAGCGATCGAAGCAGACCTGGCGCAGGCTGGGCGACTCGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACgcggc                                   <  2:1117577/150‑1 (MQ=255)
                                                                                                                   cGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGcggcgg                                 <  2:1083081/150‑1 (MQ=255)
                                                                                                                   cGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGcggcgg                                 <  1:1159027/150‑1 (MQ=255)
                                                                                                                      ggCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGtt                              <  1:440909/150‑1 (MQ=255)
                                                                                                                        cgcCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGtttt                            <  2:1185353/150‑1 (MQ=255)
                                                                                                                        cgcCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGtttt                            >  2:753870/1‑150 (MQ=255)
                                                                                                                         gcCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTc                           >  1:1319416/1‑150 (MQ=255)
                                                                                                                          ccGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCa                          >  1:1092544/1‑150 (MQ=255)
                                                                                                                            gatcATTGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCAcc                        <  2:930783/147‑1 (MQ=255)
                                                                                                                            gAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTTCGCCACGCGGCGGGTTTTCAcc                        <  2:602418/150‑1 (MQ=255)
                                                                                                                              gCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGACAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGa                      <  2:1023354/150‑1 (MQ=255)
                                                                                                                              gCATGGCCAGCATCAGGATCAACACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGa                      <  2:907670/150‑1 (MQ=255)
                                                                                                                               cATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAg                     >  2:227858/1‑150 (MQ=255)
                                                                                                                               cATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAg                     >  2:620080/1‑150 (MQ=255)
                                                                                                                               cATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAg                     <  1:322426/150‑1 (MQ=255)
                                                                                                                                    ccggCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAGGATcc                >  1:1422789/4‑150 (MQ=255)
                                                                                                                                          tCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAGGATCCAGgcgc          >  1:327617/1‑150 (MQ=255)
                                                                                                                                           cAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGTGCCAGCACAGCCAACTGCGGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAGGATCCAGGCGCt         >  1:461113/1‑150 (MQ=255)
                                                                                                                                            aGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCAtt                                                            <  2:1111223/98‑1 (MQ=255)
                                                                                                                                            aGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCAtt                                                            >  1:1111223/1‑98 (MQ=255)
                                                                                                                                             ggATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGAc                                                                                                               >  2:1064184/1‑46 (MQ=255)
                                                                                                                                             ggATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGAc                                                                                                               <  1:1064184/46‑1 (MQ=255)
                                                                                                                                             ggATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAGGATCCAGGCGCTCg       >  1:663357/1‑150 (MQ=255)
                                                                                                                                                 cAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGAGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAGGATCCAGGCGCTCGGCaa   <  2:749541/150‑1 (MQ=255)
                                                                                                                                                  aGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAGGATCCAGGCGCTCGGCAAt  <  2:1385140/150‑1 (MQ=255)
                                                                                                                                                   gCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCa                                                              >  1:1142676/1‑89 (MQ=255)
                                                                                                                                                   gCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCa                                                              <  2:1142676/89‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTCGTCGCCACTGGCTAGCGGTACGCGGGCCGACAGGTCGCCGGTTTTTTCCACGTGCAGCATCACGTCCTTGAGGGTGTTGAGTTGGCTGAGCAGGAAGCGGATCAGCAGTTGCGAGGCGCCGAGCATGGCCAGCATCAGGATCAGCACACACACCGCGTAATTGCTGAAGCGATCGAAGAAGACCTGGCGCAGGCTGGGCGACTGGGCCAGCACAGCCAACTGCTGGTCGCCATTGCGCAATACCTGTGCGCCACGCGGCGGGTTTTCACCGAGGATCCAGGCGCTCGGCAAT  >  NC_002947/2640114‑2640409

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