Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,049,192 C→A 6.1% L137L (CTG→CTT)  PP_0907 ← RND transporter membrane fusion protein‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,049,1920CA 6.1% 304.7 / 4.8 98L137L (CTG→CTT) PP_0907RND transporter membrane fusion protein‑like protein
Reads supporting (aligned to +/- strand):  ref base C (41/51);  new base A (3/3);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

ATGCGTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCAGGCTCAGATCCACGCCCTGCA  >  NC_002947/1049043‑1049340
                                                                                                                                                     |                                                                                                                                                    
aTGCGTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCAACGCCCGGTCGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGccc                                                                                                                                                      >  1:1214738/1‑150 (MQ=255)
  gCGTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAg                                                                                                                                                    <  1:280468/150‑1 (MQ=255)
  gCGTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAg                                                                                                                                                    <  2:425018/150‑1 (MQ=255)
       tCTTGATCTTCTGCGCCTTAAGCTGCTCGGCCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCgg                                                                                                                                               >  1:426811/1‑150 (MQ=255)
       tCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTGGGCCTGGCCAAGGTCga                                                                                                                                               <  1:1932107/150‑2 (MQ=255)
         ttGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTTGGcc                                                                                                                                             >  1:1049377/1‑150 (MQ=255)
          tGATCTTCTGCGCCTTAAGCTGCTCGCCCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCa                                                                                                                                            <  2:1774354/150‑1 (MQ=255)
            aTCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTGGGCCTGGCCAAGGTCgacttcg                                                                                                                                          <  2:2204854/150‑7 (MQ=255)
             tcttctGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAggg                                                                                                                                         <  1:383199/150‑1 (MQ=255)
             tcttctGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAggg                                                                                                                                         >  2:815430/1‑150 (MQ=255)
              cttctGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGc                                                                                                                                        <  1:1723201/150‑1 (MQ=255)
                tctGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCgg                                                                                                                                      <  1:1861291/150‑1 (MQ=255)
                  tGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTa                                                                                                                                    >  1:2153042/1‑150 (MQ=255)
                   gcgcCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCg                                                                                                                                                >  2:1394319/1‑137 (MQ=255)
                   gcgcCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCg                                                                                                                                                <  1:1394319/137‑1 (MQ=255)
                   gcgcCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTAc                                                                                                                                   >  2:2116527/1‑150 (MQ=255)
                      ccTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCa                                                                                                                                                     >  2:1031782/1‑129 (MQ=255)
                      ccTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCa                                                                                                                                                     <  1:1031782/129‑1 (MQ=255)
                       cTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAg                                                                                                                               >  2:1986700/1‑150 (MQ=255)
                       cTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAg                                                                                                                               >  2:2236591/1‑150 (MQ=255)
                          aaGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAg                                                                                                                            <  1:124227/150‑1 (MQ=255)
                           aGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAgg                                                                                                                           <  2:2205932/150‑1 (MQ=255)
                           aGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAgg                                                                                                                           <  1:919994/150‑1 (MQ=255)
                            gctgctCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGCTGGCCGAATCGCCCACCAGCTGGCTGCCACGGCAGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAggg                                                                                                                          >  1:1223847/1‑150 (MQ=255)
                            gctgctCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAggg                                                                                                                          <  2:35257/150‑1 (MQ=255)
                             ctgctCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGc                                                                                                                         >  2:586854/1‑150 (MQ=255)
                              tgctCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCt                                                                                                                        >  2:1004997/1‑150 (MQ=255)
                                 tCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGtt                                                                                                                     >  2:2218059/1‑150 (MQ=255)
                                 tCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGtt                                                                                                                     >  2:466636/1‑150 (MQ=255)
                                 tCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGtt                                                                                                                     >  1:808000/1‑150 (MQ=255)
                                 tCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGtt                                                                                                                     <  2:1721122/150‑1 (MQ=255)
                                   gACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTc                                                                                                                   >  1:12321/1‑150 (MQ=255)
                                    accaccGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTcc                                                                                                                  <  1:2218059/150‑1 (MQ=255)
                                      caccGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTAcc                                                                                                                                  >  1:106171/1‑132 (MQ=255)
                                      caccGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTAcc                                                                                                                                  <  2:106171/132‑1 (MQ=255)
                                       accGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGACCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGa                                                                                                               >  2:431497/1‑150 (MQ=255)
                                                 tGCGCCCGTTCTGGGAGGAAAGGCGGGCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGa                                                                                                     <  2:1299180/150‑1 (MQ=255)
                                                  gcgcCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTTGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAg                                                                                                    <  2:1049377/150‑1 (MQ=255)
                                                     ccGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTg                                                                                                 <  2:968477/150‑1 (MQ=255)
                                                        gTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCATCTCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCaa                                                                                              <  1:840865/150‑1 (MQ=255)
                                                           cTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAagag                                                                                           <  1:815430/150‑1 (MQ=255)
                                                            tGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGc                                                                                          >  2:1945384/1‑150 (MQ=255)
                                                            tGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGc                                                                                          <  2:1223231/150‑1 (MQ=255)
                                                             gggAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCa                                                                                         <  1:1349511/150‑1 (MQ=255)
                                                             gggAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCa                                                                                         <  1:1004997/150‑1 (MQ=255)
                                                             gggAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCa                                                                                         <  2:426811/150‑1 (MQ=255)
                                                              ggacgaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCAAGGCCGAAGTCCACCTTGGCCAGGCACAGGTCGGCCAGGGCGGCACCGACCAGCGCTGCTTCCTGCTCGTCGTTGAGCTCCCAGAGCAg                                                                                        >  2:717467/1‑150 (MQ=255)
                                                               gaggaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGt                                                                                       >  1:1329163/1‑150 (MQ=255)
                                                               gaggaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGACCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGt                                                                                       <  1:431497/150‑1 (MQ=255)
                                                                aggaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGtt                                                                                      >  1:592537/1‑150 (MQ=255)
                                                                  gaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGt                                                                                                                                                  >  1:917884/1‑88 (MQ=255)
                                                                  gaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGt                                                                                                                                                  <  2:917884/88‑1 (MQ=255)
                                                                  gaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGc                                                                                    <  1:2236591/150‑1 (MQ=255)
                                                                    aaGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGc                                                                                                <  2:912836/136‑1 (MQ=255)
                                                                    aaGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGc                                                                                                >  1:912836/1‑136 (MQ=255)
                                                                     aGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCg                                                                                 <  1:586854/150‑1 (MQ=255)
                                                                       gCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGc                                                                               >  1:1867838/1‑150 (MQ=255)
                                                                         ggTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCt                                                                             >  1:1448702/1‑150 (MQ=255)
                                                                          gTCGAACTCGCCGCGTGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCtt                                                                            <  1:2305791/150‑1 (MQ=255)
                                                                            cGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTGGGCCTGGCCAAGGTc                                                                                                                                                 <  1:1523180/79‑1 (MQ=39)
                                                                            cGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTGGGCCTGGCCAAGGTc                                                                                                                                                 >  2:1523180/1‑79 (MQ=39)
                                                                             gAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGAc                                                                         >  1:630213/1‑150 (MQ=255)
                                                                                cTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCg                                                                                                                                       >  2:905405/1‑85 (MQ=255)
                                                                                cTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCg                                                                                                                                       <  1:905405/85‑1 (MQ=255)
                                                                                  cgccgcGGGAGATGGCCGAAGCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGtt                                                                    <  2:12321/150‑1 (MQ=255)
                                                                                   gccgcGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTg                                                                   >  1:1027467/1‑150 (MQ=255)
                                                                                     cgcgGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGc                                                                 <  1:463917/150‑1 (MQ=255)
                                                                                     cgcgGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTAGCCCGGc                                                                               <  1:878671/136‑1 (MQ=255)
                                                                                     cgcgGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTAGCCCGGc                                                                               >  2:878671/1‑136 (MQ=255)
                                                                                          gagaTGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAg                                                            <  1:213628/150‑1 (MQ=255)
                                                                                           agaTGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGt                                                           <  2:1214738/150‑1 (MQ=255)
                                                                                            gaTGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTc                                                          <  1:1945384/150‑1 (MQ=255)
                                                                                               ggCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGaa                                                       <  2:1448702/150‑1 (MQ=255)
                                                                                               ggCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGaa                                                       <  1:1986700/150‑1 (MQ=255)
                                                                                                gCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAg                                                      >  2:1492474/1‑150 (MQ=255)
                                                                                                   gAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGc                                                   >  1:1471426/1‑150 (MQ=255)
                                                                                                        gCCCACCAGCTGTCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGa                                              >  2:1843532/1‑150 (MQ=255)
                                                                                                            accaGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTg                                          <  1:2116527/150‑1 (MQ=255)
                                                                                                                     ctgCCACGGCCGAAGTCGACCTTGGCGAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCt                                 <  2:1471426/150‑1 (MQ=255)
                                                                                                                      tgCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTc                                                          <  1:1024618/124‑1 (MQ=255)
                                                                                                                      tgCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTc                                                          >  2:1024618/1‑124 (MQ=255)
                                                                                                                          aCGGCCGAGGTCGACCTAGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGAAAGAGGCAGAGTGCAGGGTGGGGACGGTGCCTTCCAc                            >  2:1191051/1‑150 (MQ=255)
                                                                                                                             gCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTc                         <  1:1400449/150‑1 (MQ=255)
                                                                                                                             gCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTc                         <  2:1873385/150‑1 (MQ=255)
                                                                                                                               cGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCa                       <  2:1947471/150‑1 (MQ=255)
                                                                                                                                gAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCAg                      <  2:865588/150‑1 (MQ=255)
                                                                                                                                    tCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCAGGCTc                  >  1:1619711/1‑150 (MQ=255)
                                                                                                                                     cGACCTGGGCCTGGCCAAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCAGGCTCa                 >  2:212040/1‑150 (MQ=255)
                                                                                                                                      gACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGAtc                                                                                                             <  1:1021549/57‑1 (MQ=255)
                                                                                                                                      gACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGAtc                                                                                                             >  2:1021549/1‑57 (MQ=255)
                                                                                                                                      gACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTg                                                                   <  2:650634/99‑1 (MQ=255)
                                                                                                                                      gACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTg                                                                   >  1:650634/1‑99 (MQ=255)
                                                                                                                                      gACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCAGGCTCAg                >  2:2123846/1‑150 (MQ=255)
                                                                                                                                      gACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCAGGCTCAg                >  1:2321411/1‑150 (MQ=255)
                                                                                                                                      gACCTGGGCCTGGCCAAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCAGGCTCAg                >  1:1380215/1‑150 (MQ=255)
                                                                                                                                        ccTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCAGGCTCAGAt              <  1:466636/150‑1 (MQ=255)
                                                                                                                                          ttGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTACCGGTGCCTTCCACCTCCAGGCTCAGATcc            >  1:519/1‑150 (MQ=255)
                                                                                                                                            ggCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTTAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCAGGCTCAGATCCAc          <  2:2321411/150‑1 (MQ=255)
                                                                                                                                                  gCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGCTTTGCCGATGCCTTCCACCCCCAGGCTCAGATCCACGCCCTg    <  2:1223847/150‑1 (MQ=255)
                                                                                                                                                    ccAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTTCCCGGCCTTGACTTGTTGGCCTGAGTCGAAGTGCAAGCATTTAACGATTCCTTCCAACTCCAGGCTCAGATGCACGCCCTGCa  >  1:22755/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGCGTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCAGGCTCAGATCCACGCCCTGCA  >  NC_002947/1049043‑1049340

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: ‑