Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,774,386 C→T 6.7% intergenic (+38/+28) PP_5743 → / ← PP_3335 TonB‑dependent receptor protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,774,3860CT 6.7% 306.2 / 9.6 105intergenic (+38/+28)PP_5743/PP_3335TonB‑dependent receptor protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (51/47);  new base T (4/3);  total (55/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

Only 100 of 109 total aligned reads displayed.

GAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGCGTGAACGGGGT  >  NC_002947/3774239‑3774533
                                                                                                                                                   |                                                                                                                                                     
gAAACGCTGGCCGGGGTGGGGTCCAGGTATTCGGCCGCACCAGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGAcggcgg                                                                                                                                                     <  2:79189/150‑1 (MQ=255)
   aCGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTGACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGcggccg                                                                                                                                                  >  2:1739829/1‑150 (MQ=255)
     gCTGGCCGTGGTGGTGTACAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCGTGTTGGAGCTTGGGCGAAGGGTGCAAGGAAGGagaccggg                                                                                                                                                >  1:408324/1‑149 (MQ=255)
         gCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCAcc                                                                                                                                            >  2:372601/1‑150 (MQ=255)
           cGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGt                                                                                                                                          >  2:1431486/1‑150 (MQ=255)
              ggtggtgTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCaa                                                                                                                                       <  2:1929314/150‑1 (MQ=255)
               gtggtgTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATACATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGTCCGGCCACCGTCAAc                                                                                                                                      >  2:458006/1‑150 (MQ=255)
                tggtgTCCAGGTATTCGGCCGCACCCGGCAGGCTGCCAGCCTGAACCACGCAGTTGGACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACt                                                                                                                                     >  2:865537/1‑150 (MQ=255)
                     tCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCa                                                                                                                                >  1:1924867/1‑150 (MQ=255)
                        aGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCTGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATgg                                                                                                                             <  1:926253/150‑1 (MQ=255)
                        aGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCGAGGACGGCGGCCGGCCACCGTCAACTGATCATgg                                                                                                                             <  2:450639/150‑1 (MQ=255)
                          gTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGtt                                                                                                                           >  2:1017658/1‑150 (MQ=255)
                                 gCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCgg                                                                                                                    <  1:1125880/150‑1 (MQ=255)
                                         gttCAGGCTGGCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCaa                                                                                                                                       >  2:536619/3‑123 (MQ=255)
                                         gttCAGGCTGGCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCaa                                                                                                                                       <  1:536619/121‑1 (MQ=255)
                                          gTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGca                                                                                                           >  1:2073124/1‑150 (MQ=255)
                                           tCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGcag                                                                                                          >  2:891342/1‑150 (MQ=255)
                                              ggCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGcagcca                                                                                                       <  1:396915/150‑1 (MQ=255)
                                                cTGCCAGCCTGAACCACGCAGTGGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGcagccagc                                                                                                     <  1:946158/150‑1 (MQ=255)
                                                   ccAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGcc                                                                                                  >  2:41905/1‑150 (MQ=255)
                                                    cAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGTGGCCGGCCGCCGTcctt                                                                                                                                      >  1:1724566/1‑110 (MQ=255)
                                                    cAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGTGGCCGGCCGCCGTcctt                                                                                                                                      <  2:1724566/113‑4 (MQ=255)
                                                        cTGAACCACGCAGTTGGACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAg                                                                                           <  1:865537/150‑1 (MQ=255)
                                                        cTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAg                                                                                           <  1:372601/150‑1 (MQ=255)
                                                        cTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGAACGGCGGCTGGCCACCGTCAACTGATTATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAg                                                                                           <  1:473742/150‑1 (MQ=255)
                                                          gAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCt                                                                                         >  1:2191779/1‑150 (MQ=255)
                                                              cACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGTTAACCTCATGTTGGCGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTAccc                                                                                     >  1:267892/1‑150 (MQ=255)
                                                              cACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGc                                                                                          >  1:1177225/1‑145 (MQ=255)
                                                              cACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGc                                                                                          <  2:1177225/145‑1 (MQ=255)
                                                              cACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTAccc                                                                                     >  1:341132/1‑150 (MQ=255)
                                                               aCGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCcg                                                                                    >  2:133719/1‑150 (MQ=255)
                                                                 gCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCcgcg                                                                                  >  1:2296161/1‑150 (MQ=255)
                                                                   aGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGcc                                                                                >  1:202077/1‑150 (MQ=255)
                                                                     ttGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCtt                                                                              <  1:1169320/150‑1 (MQ=255)
                                                                      tGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCttt                                                                             <  2:202077/150‑1 (MQ=255)
                                                                      tGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCttt                                                                             <  1:41905/150‑1 (MQ=255)
                                                                        aaCGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATgg                                                                                                                             <  2:1809691/102‑1 (MQ=255)
                                                                        aaCGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATgg                                                                                                                             >  1:1809691/1‑102 (MQ=255)
                                                                             gcgcATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAg                                                                                           >  2:2174294/1‑129 (MQ=255)
                                                                               gcATGGTTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAg                                                                    >  1:1328292/1‑150 (MQ=255)
                                                                               gcATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTc                                                                                                                                         >  2:389293/1‑83 (MQ=255)
                                                                               gcATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTc                                                                                                                                         <  1:389293/83‑1 (MQ=255)
                                                                                  tGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGTAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCgg                                                                 >  1:948490/1‑150 (MQ=255)
                                                                                   ggCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCggg                                                                <  2:2296161/150‑1 (MQ=255)
                                                                                      ttGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCgg                                                                                                                    <  1:849398/97‑1 (MQ=255)
                                                                                      ttGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCgg                                                                                                                    >  2:849398/1‑97 (MQ=255)
                                                                                        gAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTAc                                                                                       >  1:1090834/1‑122 (MQ=255)
                                                                                        gAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTAc                                                                                       <  2:1090834/122‑1 (MQ=255)
                                                                                         aaGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAg                                                          <  1:1431486/150‑1 (MQ=255)
                                                                                           gTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAgcgc                                                                                                   <  2:239618/109‑1 (MQ=255)
                                                                                           gTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAgcgc                                                                                                   >  1:239618/1‑109 (MQ=255)
                                                                                            tACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGgag                                                       <  1:2305745/150‑1 (MQ=255)
                                                                                             aCATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGgaga                                                      <  2:1924867/150‑1 (MQ=255)
                                                                                               atccatcGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAAc                                                                                                                                      <  2:2182556/70‑1 (MQ=255)
                                                                                               atccatcGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAAc                                                                                                                                      >  1:2182556/1‑70 (MQ=255)
                                                                                               atccatcGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAAc                                                    >  2:2343211/1‑150 (MQ=255)
                                                                                               atccatcGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAAc                                                    >  1:2106404/1‑150 (MQ=255)
                                                                                                   atcGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTAGCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTccc                                                >  1:22330/1‑150 (MQ=255)
                                                                                                   atcGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTccc                                                >  1:1545496/1‑150 (MQ=255)
                                                                                                   atcGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTccc                                                >  2:1344900/1‑150 (MQ=255)
                                                                                                    tcgtcgGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCa                                               >  2:1680942/1‑150 (MQ=255)
                                                                                                    tcgtcgGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCa                                               >  1:798761/1‑150 (MQ=255)
                                                                                                    tcgtcgGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCa                                               <  2:798761/150‑1 (MQ=255)
                                                                                                    tcgtcgGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCa                                               <  1:2068931/150‑1 (MQ=255)
                                                                                                     cgtcgGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAg                                              <  1:1017658/150‑1 (MQ=255)
                                                                                                            aaCCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACggtg                                       <  2:1814863/150‑1 (MQ=255)
                                                                                                              ccccATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATc                                                                                                                                 >  1:2401271/4‑60 (MQ=255)
                                                                                                              ccTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGGCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGGCTTTGCAATCGAGCGGCCTCCACGAGAACTCCCCGGACggttgt                                     <  2:642992/150‑1 (MQ=255)
                                                                                                              ccTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATc                                                                                                                                 <  2:2401271/60‑1 (MQ=255)
                                                                                                               cTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGTCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGTGACCGGCAGAACCTGATGG‑‑AAGCGACTCCCTCCGGAGCCATCGAGCGGGCTCCAGGAGAACTCACAGGCCGgtggtg                                    >  1:284931/1‑150 (MQ=255)
                                                                                                               cTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATg                                                                                                                              <  2:64955/62‑1 (MQ=255)
                                                                                                               cTCATGTTGGAGCTTGGGCAAAGGGTGAAAGGACGGCGGCCGGCCACCGTCAACTGATCATg                                                                                                                              >  1:64955/1‑62 (MQ=255)
                                                                                                                tCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGtggtgg                                   >  1:1815159/1‑150 (MQ=255)
                                                                                                                tCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGtggtgg                                   >  1:2294070/1‑150 (MQ=255)
                                                                                                                         aGCTTGGGCAAAGGGTGCAATCACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGccctc                          <  1:917751/150‑1 (MQ=255)
                                                                                                                          gCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGccctcc                         <  2:1815159/150‑1 (MQ=255)
                                                                                                                          gCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGccctcc                         <  2:1328292/150‑1 (MQ=255)
                                                                                                                            ttGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGccctcccc                       >  1:1639164/1‑150 (MQ=255)
                                                                                                                               ggCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGccctccccctc                    >  1:1737745/1‑150 (MQ=255)
                                                                                                                                gCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGa                                            >  1:1176579/1‑125 (MQ=255)
                                                                                                                                gCAAAGGGTGCAAGGACGGCGGCCGGCCACAGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGCGGCAAGccctccccctcc                   >  2:125960/1‑150 (MQ=255)
                                                                                                                                  aaaGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTg                 >  2:808753/1‑150 (MQ=255)
                                                                                                                                  aaaCGGGGCAAGGACGGCGACCGGCCACCGTCAACTGATCATGGATTAGCCGGGCACCGGCAGCCAGCGCCGG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCACGACGGTGGTGGCAAGCCCTCCCCCTCCTg                 >  2:601466/1‑150 (MQ=255)
                                                                                                                                    ggggTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTggg               <  2:1540262/149‑1 (MQ=255)
                                                                                                                                    aGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGCTCCTGGTATAGGGCTAAGCGCCTGATAGGCATGGCGAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTggg               <  1:891342/150‑1 (MQ=255)
                                                                                                                                    aGGGTGCAAGGACGGCGCCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGATAACTCCCAGGACGGTGGTGGCAAGCCCTCCCACTCCTggg               <  2:408324/150‑1 (MQ=255)
                                                                                                                                     gggTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTgggg              >  2:1689918/1‑150 (MQ=255)
                                                                                                                                      ggTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGc                                                                   >  2:680494/1‑96 (MQ=255)
                                                                                                                                      ggTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGc             >  1:1029907/1‑150 (MQ=255)
                                                                                                                                      ggTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGc             <  1:1334960/150‑1 (MQ=255)
                                                                                                                                        tGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGCGt           <  1:458006/150‑1 (MQ=255)
                                                                                                                                           aaGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGCGTGaa        >  1:2344448/1‑150 (MQ=255)
                                                                                                                                            ttgACGGTGGCCGGCCGCCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACggtggt                                     <  2:1843195/118‑1 (MQ=255)
                                                                                                                                            ttgACGGTGGCCGGCCGCCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACggtggt                                     >  1:1843195/3‑120 (MQ=255)
                                                                                                                                             ggACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGCGTGAACg      >  2:1247412/1‑150 (MQ=255)
                                                                                                                                              gACGGTGGCCGGCCGCCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTgg                <  2:1560380/139‑1 (MQ=255)
                                                                                                                                              gACGGTGGCCGGCCGCCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTACCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGCGTGAAcag     >  1:610719/1‑148 (MQ=255)
                                                                                                                                              gACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGCGTGAACgg     <  1:133719/150‑1 (MQ=255)
                                                                                                                                               aCGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGCGTGAACggg    <  2:2294070/150‑1 (MQ=255)
                                                                                                                                                 ggcggcCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGCGTGAACGGGGt  >  2:1116210/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGCCTG‑‑AAGCTACCCGCGCCTTTGCCATCGAGCGGGCTCCAGGAGAACTCCCAGGACGGTGGTGGCAAGCCCTCCCCCTCCTGGGGCGTGAACGGGGT  >  NC_002947/3774239‑3774533

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: