Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,992,339 G→C 5.1% I259M (ATC→ATG PP_2616 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,992,3390GC 5.1% 242.0 / 4.0 79I259M (ATC→ATGPP_2616hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/31);  new base C (2/2);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

AGCGAGGCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTGGTGAATCAATGC  >  NC_002947/2992190‑2992488
                                                                                                                                                     |                                                                                                                                                      
aGCGAGGCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCatggcg                                                                                                                                                        >  1:140969/1‑150 (MQ=32)
 gCGAGGCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCatggcga                                                                                                                                                       >  1:487362/1‑150 (MQ=32)
     ggCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCatggcgatggc                                                                                                                                                   <  2:1236976/150‑1 (MQ=32)
      gCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCatggcgatggcg                                                                                                                                                  >  1:792309/1‑150 (MQ=32)
      gCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCatggcgatggcg                                                                                                                                                  >  2:517387/1‑150 (MQ=32)
        aaaaaCCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCGCCGCCATGGCGATGGCGGa                                                                                                                                                >  2:1272664/1‑150 (MQ=17)
               tCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGTCACCGCCATGGCGATGGCGGAAATCAcc                                                                                                                                         >  2:1389615/1‑150 (MQ=38)
               tCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGAGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCAcc                                                                                                                                         >  1:801370/1‑150 (MQ=255)
               tCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCGATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCAcc                                                                                                                                         >  2:876835/1‑150 (MQ=38)
                 gcgAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCaa                                                                                                                                       >  2:396596/1‑150 (MQ=39)
                 gcgAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCaa                                                                                                                                       >  1:1148310/1‑150 (MQ=38)
                        cGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTc                                                                                                                                <  1:396596/150‑1 (MQ=40)
                         gTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATCGCCATGGCggtg                                                                                                                                               >  2:58096/1‑132 (MQ=255)
                         gTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATCGCCATGGCggtg                                                                                                                                               <  1:58096/134‑3 (MQ=12)
                           tGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTccac                                                                                                                             <  2:140969/150‑1 (MQ=40)
                                tCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGTCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAggg                                                                                                                        <  1:1389615/150‑1 (MQ=255)
                                  aGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGcc                                                                                                                      >  1:695770/1‑150 (MQ=255)
                                  aGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGcc                                                                                                                      >  2:1287696/1‑150 (MQ=255)
                                    tggtgCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCAt                                                                                                                    <  2:1622204/150‑1 (MQ=255)
                                          tgtgACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCgg                                                                                                              >  2:228864/1‑150 (MQ=255)
                                                acGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATg                                                                                                        <  2:1606706/150‑1 (MQ=255)
                                                 cGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGc                                                                                                       <  1:1384540/150‑1 (MQ=255)
                                                    ttGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGa                                                                                                    >  1:589399/1‑150 (MQ=255)
                                                           ttGCTGCAAAATGCCGTTCAGCGACTTCGACACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTACACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCg                                                                                             >  1:1021734/1‑150 (MQ=255)
                                                           ttGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCg                                                                                             >  1:1919073/1‑150 (MQ=255)
                                                               tgcAAAATGCCTTTCAGCGACTCCAGCACCGCGTTTTTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACAGCCATGGCGATGGCGGAAATCACCAAGGCGCCCCACCAGGGCCATTTCCGGCTCATGCGGAAGTTGCGTggg                                                                                        >  2:1088661/1‑150 (MQ=17)
                                                                         cGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTgg                                                                              >  2:1344312/1‑150 (MQ=255)
                                                                          gTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTggg                                                                                        >  1:1321931/1‑139 (MQ=255)
                                                                          gTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTggg                                                                                        <  2:1321931/139‑1 (MQ=255)
                                                                          gTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGgt                                                                             <  1:876835/150‑1 (MQ=255)
                                                                          gTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCAATGCCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCGCATGCGGAAGTTGCGTGGGGCCACATTGgt                                                                             >  1:124668/1‑150 (MQ=255)
                                                                           ttCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAgg                                                                                                                         >  2:1034238/1‑106 (MQ=255)
                                                                           ttCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAgg                                                                                                                         <  1:1034238/106‑1 (MQ=255)
                                                                              aGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATg                                                                                                        <  2:121820/120‑1 (MQ=255)
                                                                              aGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATg                                                                                                        >  1:121820/1‑120 (MQ=255)
                                                                                     ccAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGc                                                                  >  1:272465/1‑150 (MQ=255)
                                                                                      cAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAg                                                                         >  1:1769091/1‑142 (MQ=255)
                                                                                      cAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAg                                                                         <  2:1769091/142‑1 (MQ=255)
                                                                                            cTCGGTTGTGATCAGGTGCAGGCGGTAGCTGGAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGActtc                                                           >  1:74828/1‑150 (MQ=255)
                                                                                             tCGGTTGTGATCAGGTGCAGGCGGTAGCTGTATATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGcc                                                                                                                      <  2:1233519/91‑1 (MQ=255)
                                                                                             tCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGcc                                                                                                                      >  1:1233519/1‑91 (MQ=255)
                                                                                                 ttGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCgg                                                      >  2:1221981/1‑150 (MQ=255)
                                                                                                 ttGTGATCAGGTGCAGGAGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCgg                                                      <  2:801370/150‑1 (MQ=255)
                                                                                                  tgtgATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCggg                                                     <  1:1344312/150‑1 (MQ=255)
                                                                                                   gtgATCAGGTGGAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCgggg                                                    <  2:95889/150‑1 (MQ=255)
                                                                                                       tCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAg                                                >  2:357204/1‑150 (MQ=255)
                                                                                                       tCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAg                                                >  2:1207788/1‑150 (MQ=255)
                                                                                                       tCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAg                                                <  2:695770/150‑1 (MQ=255)
                                                                                                       tCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCAAAGAGCGTCGTGTTGGGGCCACATTGGTGGAGGGGTCGCAGACTTCTTCGGGGGCAg                                                >  2:51791/1‑150 (MQ=17)
                                                                                                       tCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGAACGGAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAg                                                <  1:1088661/150‑1 (MQ=255)
                                                                                                          ggTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGGCTTCTTCGGGGGTAGGGc                                             >  2:755390/1‑150 (MQ=255)
                                                                                                           gTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGcc                                            <  2:589399/150‑1 (MQ=255)
                                                                                                             gCAGGCGGTAGTTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCcgc                                          <  1:731003/150‑1 (MQ=255)
                                                                                                             gCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCcgc                                          >  1:768074/1‑150 (MQ=255)
                                                                                                             gCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCAACCGTCGCCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCcgc                                          >  1:1320701/1‑150 (MQ=18)
                                                                                                               aGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCgcgc                                        >  2:108074/1‑150 (MQ=255)
                                                                                                                 gCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAg                                      >  2:190564/1‑150 (MQ=255)
                                                                                                                  cGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAg                                                                                                  <  1:1098459/90‑1 (MQ=255)
                                                                                                                  cGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAg                                                                                                  >  2:1098459/1‑90 (MQ=255)
                                                                                                                       gCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTcgcg                                <  2:768074/150‑1 (MQ=255)
                                                                                                                       gCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTcgcg                                <  1:1287696/150‑1 (MQ=255)
                                                                                                                         tGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGga                              <  1:1272664/150‑1 (MQ=255)
                                                                                                                          gTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTg                                                                               >  1:1695866/1‑100 (MQ=255)
                                                                                                                          gTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTg                                                                               <  2:1695866/100‑1 (MQ=255)
                                                                                                                          gTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCgggg                                                    <  1:303273/127‑1 (MQ=255)
                                                                                                                          gTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCgggg                                                    >  2:303273/1‑127 (MQ=255)
                                                                                                                          gTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCAATGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGGGTAGGGGTCGCAGACTTCTACGAAAGAAGGGCCGCGCAGTTCGCGgat                             >  2:449153/1‑150 (MQ=18)
                                                                                                                          gTAGATCCCGTAGGCCACCGCCATGCCGAGGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTACCGGCTCATGCGGAAGTTGCGTGGGTCCACATTGGTGTAGAGGGCGCAGACTCGTCCGGGGGTAGGGCCGCGCAGTGCGCGgat                             >  2:439183/1‑150 (MQ=18)
                                                                                                                           tAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAggg                                                                       <  2:1069528/107‑1 (MQ=255)
                                                                                                                           tAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGGGGGGCCACATTGGTGTAggg                                                                       >  1:1069528/1‑107 (MQ=255)
                                                                                                                             gATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCAGGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTGGCAGACTTCTTCGCGGGTAGGGACGCGCAGTTCGCGgatgat                          >  2:578661/1‑150 (MQ=255)
                                                                                                                               tCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTa                                                 >  2:275311/1‑125 (MQ=255)
                                                                                                                               tCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTa                                                 <  1:275311/125‑1 (MQ=255)
                                                                                                                                 ccGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCCTTGCTGGCTCATGCGAACGTTGCGGGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCg                      >  1:726257/1‑150 (MQ=17)
                                                                                                                                  cGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGt                     >  1:637977/1‑150 (MQ=255)
                                                                                                                                     aGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCCAGGACGACGGCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCa                  <  2:1021734/150‑1 (MQ=18)
                                                                                                                                       gCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCACAGTTCGCGGATGATCCCGTGCAGt                >  1:1811418/1‑150 (MQ=255)
                                                                                                                                         cACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTg              >  2:1163223/1‑150 (MQ=255)
                                                                                                                                          tccGCCATCGCCATGGCGGAAATCACCAAGGGGCTCCACCAGGGCAATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAggg                                                                       <  1:859605/91‑1 (MQ=37)
                                                                                                                                          tccGCCATCGCCATGGCGGAAATCACCAAGGGGCTCCACCAGGGCAATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAggg                                                                       >  2:859605/2‑92 (MQ=255)
                                                                                                                                          aCCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTgg             <  1:517387/150‑1 (MQ=255)
                                                                                                                                          aCCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTgg             <  1:108074/150‑1 (MQ=255)
                                                                                                                                              ccATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTGGTGaa         >  2:1354005/1‑150 (MQ=255)
                                                                                                                                              ccATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTGGTGaa         <  1:1221981/150‑1 (MQ=255)
                                                                                                                                                     gatggcgGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCAAATAGGCGTAGGGGTCGCAGACTTCTTCGGGGGTATGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTGGTGAATCAATGc  >  1:919565/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGCGAGGCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTGGTGAATCAATGC  >  NC_002947/2992190‑2992488

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: