Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,992,336 G→C 5.1% A260A (GCC→GCG PP_2616 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,992,3360GC 5.1% 243.6 / 4.0 79A260A (GCC→GCGPP_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

GTCAGCGAGGCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTGGTGAA  >  NC_002947/2992187‑2992481
                                                                                                                                                     |                                                                                                                                                  
gTCAGCGAGGCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAACATGCGGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCACGTAGGCCCCCGCCatg                                                                                                                                                    >  1:1613271/1‑150 (MQ=255)
gTCAGCGAGGCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCatg                                                                                                                                                    <  2:652276/150‑1 (MQ=32)
   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                                                                                                                                        <  1:58096/134‑3 (MQ=12)
                            gTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATCGCCATGGCggtg                                                                                                                                        >  2:58096/1‑132 (MQ=255)
                              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)
                                                              ttGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCatg                                                                                                                                                    >  1:1450094/1‑88 (MQ=35)
                                                              ttGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCatg                                                                                                                                                    <  2:1450094/88‑1 (MQ=35)
                                                              ttGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCg                                                                                      >  1:1919073/1‑150 (MQ=255)
                                                                  tgcAAAATGCCTTTCAGCGACTCCAGCACCGCGTTTTTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACAGCCATGGCGATGGCGGAAATCACCAAGGCGCCCCACCAGGGCCATTTCCGGCTCATGCGGAAGTTGCGTggg                                                                                 >  2:1088661/1‑150 (MQ=17)
                                                                            cGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTgg                                                                       >  2:1344312/1‑150 (MQ=255)
                                                                             gTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTggg                                                                                 <  2:1321931/139‑1 (MQ=255)
                                                                             gTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTggg                                                                                 >  1:1321931/1‑139 (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:695770/150‑1 (MQ=255)
                                                                                                          tCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAg                                         >  2:357204/1‑150 (MQ=255)
                                                                                                          tCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAg                                         >  2:1207788/1‑150 (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                         <  1:1287696/150‑1 (MQ=255)
                                                                                                                          gCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTcgcg                         <  2:768074/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                                             >  2:303273/1‑127 (MQ=255)
                                                                                                                             gTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCgggg                                             <  1:303273/127‑1 (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  <  1:1221981/150‑1 (MQ=255)
                                                                                                                                                 ccATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTGGTGaa  >  2:1354005/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTCAGCGAGGCAAAAACCTCGCGAATTCGTTGAACTCAGTGGTGCTGTGACACGGTTGGCTATTGCTGCAAAATGCCGTTCAGCGACTCCAGCACCTCGGTTGTGATCAGGTGCAGGCGGTAGCTGTAGATCCCGTAGGCCACCGCCATGGCGATGGCGGAAATCACCAAGGGGCTCCACCAGGGCCATTGCCGGCTCATGCGGAAGTTGCGTGGGGCCACATTGGTGTAGGGGTCGCAGACTTCTTCGGGGGTAGGGCCGCGCAGTTCGCGGATGATCCCGTGCAGTTGGTGAA  >  NC_002947/2992187‑2992481

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: