Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,929,914 T→G 7.6% intergenic (‑42/+111) PP_5195 ← / ← fbpA iron ABC transporter permease/iron ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,929,9140TG 7.6% 495.2 / 7.7 177intergenic (‑42/+111)PP_5195/fbpAiron ABC transporter permease/iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (64/95);  new base G (16/2);  total (80/97)
Fisher's exact test for biased strand distribution p-value = 8.21e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

CAGGTCGATGGATTGCCAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATAAAGCCTGGCGCTTAGTTGTAGCCAGCCCGATCCATCAGGCGAATGGCTT  >  NC_002947/5929765‑5930062
                                                                                                                                                     |                                                                                                                                                      
aggaagtggcggcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggcccttt                                                                                                                                                        >  1:584146‑M1/149‑150 (MQ=255)
 aGGTCGATGGATTGCCAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTc                                                                                                                                                       <  2:1124192/150‑1 (MQ=255)
   aagtggcggcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTcgc                                                                                                                                                    >  1:925034‑M1/146‑150 (MQ=255)
   aagtggcggcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTcgc                                                                                                                                                    >  2:1810127‑M1/146‑150 (MQ=255)
    tCGATGGAGTGCCAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTcgcg                                                                                                                                                   <  1:1767711/150‑1 (MQ=255)
    agtggcggcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTcgcg                                                                                                                                                   >  2:1982168‑M1/145‑150 (MQ=255)
     gtggcggcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCgg                                                                                                                                                  >  2:1860259‑M1/144‑150 (MQ=255)
     gtggcggcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCgg                                                                                                                                                  >  1:2111118‑M1/144‑150 (MQ=255)
     gtggcggcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCgg                                                                                                                                                  >  1:4077551‑M1/144‑150 (MQ=255)
     gtggcggcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCgg                                                                                                                                                  >  2:443596‑M1/144‑150 (MQ=255)
     cGATGGATTGCCAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCgg                                                                                                                                                  <  2:2494534/150‑1 (MQ=255)
     cGATGGATTGCCAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCgg                                                                                                                                                  <  2:3530455/150‑1 (MQ=255)
      gATGGATTGCCAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTGGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCggg                                                                                                                                                 >  2:3333430/1‑150 (MQ=255)
          ggcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTaaa                                                                                                                                             <  1:3629996‑M1/12‑1 (MQ=255)
           gcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAAc                                                                                                                                            >  1:3362082‑M1/138‑150 (MQ=255)
           gcgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaaacgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAAc                                                                                                                                            >  2:3257891‑M1/138‑150 (MQ=255)
           aTTGCCAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAAc                                                                                                                                            <  1:2680293/150‑1 (MQ=255)
            ttGCCAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAAcc                                                                                                                                           <  1:4439133/150‑1 (MQ=255)
            cgtggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAAcc                                                                                                                                           >  1:725337‑M1/137‑150 (MQ=255)
              tggggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAACCCg                                                                                                                                         >  1:3131178‑M1/135‑150 (MQ=255)
              gCCAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGGGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTGCGCGGGGGAAcccc                                                                                                                                         >  1:1556976/1‑149 (MQ=255)
               ccAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTGGCGGGGAAAcaccc                                                                                                                                        >  1:273214/1‑146 (MQ=255)
               ccAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGc                                                                                                                                        <  2:3284704/150‑1 (MQ=255)
               ccAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGc                                                                                                                                        <  2:3895536/150‑1 (MQ=255)
                gggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAACCCGCt                                                                                                                                       <  2:2638168‑M1/18‑1 (MQ=255)
                gggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAACCCGCt                                                                                                                                       <  2:522119‑M1/18‑1 (MQ=255)
                gggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAACCCGCt                                                                                                                                       <  1:2275603‑M1/18‑1 (MQ=255)
                 ggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAACCCGCTc                                                                                                                                      >  1:3626819‑M1/132‑150 (MQ=255)
                 ggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAACCCGCTc                                                                                                                                      >  1:1153348‑M1/132‑150 (MQ=255)
                 ggcagcttcaaggccgacagcattccggtggaaattgccggcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAACCCGCTc                                                                                                                                      >  1:3792140‑M1/132‑150 (MQ=255)
                 aGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTGCGCGGGTAAAcccccac                                                                                                                                      >  2:3863597/1‑146 (MQ=255)
                     caacaGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                  <  2:1462101/150‑1 (MQ=255)
                      aacaGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACg                                                                                                                                 <  2:1622279/150‑1 (MQ=255)
                       acaGAATCACTCTCATGGTCAGCATCACCAATGCGGCGTTCAGGGAGACCGGGAGGTACCATCGGCGTTGGGTGGTTTGGGGCAAGTCTGAAATCTCGAGTTCGATGGCGACGGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACgt                                                                                                                                <  2:686981/150‑2 (MQ=255)
                        cagcagCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACggg                                                                                                                               <  1:2541850/150‑1 (MQ=255)
                          gcagcaCGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGa                                                                                                                             <  2:3156971/150‑1 (MQ=255)
                              caCGCTCAGGGGCAGCAGCACCACGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGTGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTAGTTCGCGTGTAAACCCGCTCCTACGGGGACcgc                                                                                                                         <  1:2458739/150‑1 (MQ=255)
                              caCGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACcgc                                                                                                                         <  1:4226565/150‑1 (MQ=255)
                              caCGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACcgc                                                                                                                         <  1:138198/150‑1 (MQ=255)
                               aCGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACcgcg                                                                                                                        <  1:36602/150‑1 (MQ=255)
                                  cTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCaa                                                                                                                     <  1:1729224/150‑1 (MQ=255)
                                  cTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTGCGCGGGTAAAACCCCTCCAACGGGGACCCCGCac                                                                                                                     >  2:3822277/1‑149 (MQ=255)
                                   tCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAg                                                                                                                    <  2:4462917/150‑1 (MQ=255)
                                    gatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctttcgcgggtaaacccgctcctacaggatctgcactgccCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAgg                                                                                                                   >  1:2391191‑M1/102‑150 (MQ=255)
                                       gggCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGGAAACCCCCTCCGACGGGGCCCGCGCAAGGccg                                                                                                                >  2:2851061/1‑148 (MQ=255)
                                       gggCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGGTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTGGGGGGGGTACCCCCCCCCCACGGGGCCCccgccgggcgg                                                                                                                >  1:74337/1‑143 (MQ=255)
                                        ggCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGt                                                                                                               <  2:4200221/150‑1 (MQ=255)
                                           agcagcaCCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGGTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTGCGCGGGTGAACCCCCACCGACCGGGCCCCCCCAACGcccggca                                                                                                            >  2:534746/1‑144 (MQ=255)
                                             cagcaCCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAAAACCCGCCGACCGGGGCCCCCCCACGCAGTGCAGa                                                                                                          >  2:383737/1‑150 (MQ=255)
                                               gcaCCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTGGCGGGTAAAACCCCTCCTAACGGGACCCCCCCAGGCCGGGCAGaaa                                                                                                        >  2:1245944/1‑148 (MQ=255)
                                               gcaCCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATc                                                                                                        <  1:3333546/150‑1 (MQ=255)
                                               gcaCCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCCGGTAAACCCCCACCAACGGGGACCCCGCAAGGCAGTGCAGATc                                                                                                        >  2:1251076/1‑150 (MQ=255)
                                                  ccATGGCGTCGATCTGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGTTGTGGGGCACGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAATCCCGCTCCTACGGGGACCGCGCAAGGCTGTGAAGTTACTg                                                                                                     <  2:908491/150‑1 (MQ=255)
                                                    aaGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTa                                                                                                   <  2:674914/150‑1 (MQ=255)
                                                     aGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAg                                                                                                  <  1:2670173/150‑1 (MQ=255)
                                                         gcaagcgccaggctgaagccattcgcctgatggatcgggctggctacaactaagcgccaggctttatcctctttgggcccgttagggccctTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGc                                                                                              >  2:1204071‑M1/92‑150 (MQ=255)
                                                           cgATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCgg                                                                                            <  2:3464424/150‑1 (MQ=255)
                                                           cgATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCgg                                                                                            <  2:2569939/150‑1 (MQ=255)
                                                           cgATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAAACCGCACCAACGGGGCCCCCCCAAGGCACTGCACAACCAGTAGGAGCgg                                                                                            >  2:3812777/1‑150 (MQ=255)
                                                            gATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCggg                                                                                           <  1:2882217/150‑1 (MQ=255)
                                                            gATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCGCCCACCGGGACCGCCCAAGGCACGGCAGATCCGGGAGGATCggg                                                                                           >  1:2635519/1‑150 (MQ=255)
                                                                  ggAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTGGGCGGGTTAAAACCCACCAACGGGGGCCCCGCCAGGCCCGGCAGAAACATTATGAGCGGGTTTaac                                                                                     >  1:4294283/1‑148 (MQ=255)
                                                                  ggAGACCGGGAGGTACCAGCGGCGGTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGGGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCCCACCAACGGGGACCGCGCCAGGCCGTGCAGATCCGGGAGGAGCGGGTTTAcc                                                                                     >  2:1426358/1‑150 (MQ=255)
                                                                   gagaCCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTGGCGGGTAAAACCGCTACTACGGGGACAGCGCAAGGCAGTGCAGATAGTGTATGAGCGGGTTTAacc                                                                                   >  2:276175/1‑150 (MQ=255)
                                                                   gagaCCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTAccc                                                                                   <  2:1536470/150‑1 (MQ=255)
                                                                        cGGGAGGTACCAGCGGCGTTGGGGGTGTGGGGCAAGGCTGAAATCTCGTGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGaaa                                                                             <  1:2578561/150‑1 (MQ=255)
                                                                              gTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCACCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                        >  2:3348186/1‑150 (MQ=255)
                                                                              gTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCCCGGGTAAAACCCCACCAACCGGGACCGCCCAACGCAGTGCAGATCCGGTAGGAGCGGGTTTACCCGCGCAAGGGcc                                                                        >  2:289631/1‑150 (MQ=255)
                                                                              gTACCAGCGGCGGTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTAGCGGGTAAAACCGCTCCAAAGGGGTCCGCGCAATGCAGTGGAGATACTGTAGGTGGGGGTTTAACCGCGAAAggggc                                                                        >  2:961372/1‑148 (MQ=255)
                                                                                  ccgcggggTGGGGGGGTGTGGGGAACGGCTTAAATCTCGAGTGCGAGGGCGACCGCGGCGGGGTGGTTGGCGGGTAACCCCGCTCCTACGGGGACCGCGCACGGCAGTGCAGATCCTGTCGGAGGGGTTTTACCCGCGAAAGGGCCCTaa                                                                    <  1:1941163/148‑1 (MQ=255)
                                                                                  cAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCCGGTAAAACCCCACCAAAGGGGACCCCGCAACGCACGGCACATACGGTAGGAGCCGGGTTAACCGCGAAAGGGCCCGaa                                                                    >  2:1336060/1‑150 (MQ=255)
                                                                                   aGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAAACCGCTCCTACGGGGACCCCGCAAGGCCGTGGATATACTGTAGGAGCGGGTTTTCCCCCGAAAGGGCCCAAAc                                                                   >  2:2410940/1‑150 (MQ=255)
                                                                                    gcggcgTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTGGGCCGGGGAAACCCCCACCAAGGGGGCCCCCCCAGGCCGGGCACAAGCGGGAGGAGGGGGTTTTCCACCCAAAAGGGccagacg                                                                   >  1:1652163/1‑145 (MQ=255)
                                                                                     cggcgTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGTTGGCGACCGCTGCGCGGTCGTTCGCGGTTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACgg                                                                 <  1:3282787/150‑1 (MQ=255)
                                                                                           tGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCTGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCaa                                                           <  1:3757547/150‑1 (MQ=255)
                                                                                              ggggTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCAAaga                                                        <  2:1709411/150‑1 (MQ=255)
                                                                                              ggggTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCCCCCCAACGGGGGCCGCGCCAGGCCGGGCAGATCCGGGATGGGCGGGTTTAGCCGCGCAAGGGCCCGAACCGGCCCCAaga                                                        >  1:2940117/1‑150 (MQ=255)
                                                                                                ggTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAgg                                                      <  2:920039/150‑1 (MQ=255)
                                                                                                     ggggCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATaaa                                                 <  2:2467252/150‑1 (MQ=255)
                                                                                                       ggCATGGCTGAATTCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCTGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATAAAGc                                               <  1:3743907/150‑1 (MQ=255)
                                                                                                        gCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCCCGCCGAAGGGGACCCCCCAAGGCAGTGCCGATCCAGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGAACGGGCCCCAAGAGGATAAAGcc                                              >  2:733201/1‑150 (MQ=255)
                                                                                                          aaGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCCCGCCTACGGGGTCCGCGCAAGGCAGTGCAGATACGGTAGGAGCGGGTTTACCCGCCCAAGGGCCCTAACCGGCCCCAAGAGGATAAAGccag                                            >  2:2166036/1‑148 (MQ=255)
                                                                                                             gCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCCCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCAAACGGGCCCAAAGAGGATAAAGCCTGgcg                                         >  2:3300348/1‑150 (MQ=255)
                                                                                                              cTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCCCACCAACCGGGACCCCCCCAGGCCCTGCAGATCCGGTAGGAGCGGGTTTACCCCCCAAAGGGCCCCAACGGGCCCAAAGAGGATAAAACCAGgcgc                                        >  2:3629974/1‑150 (MQ=255)
                                                                                                               tGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAAACCCCTCCTAAGGGGCCCGCGCAAGGCAGGGCCGATACGGTAGGAGCGGGTTTTACCCCCCAATGGCCCAAAAGGGCCCAAAGAGGATAAAACCCGgcgcg                                       >  2:1582815/1‑149 (MQ=255)
                                                                                                               tGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTGCGCCGGTTAAACCCCACCCAACGGGGCCCCCCAACGCCCGGGACATACGGTTGGGGCCGGTGTGCCCGCCAAACGGCCCCAACCGGCCCCAAAAGGGTAAAACCGGgggcg                                       >  2:754670/1‑149 (MQ=255)
                                                                                                                 aaaTCTCGTGTGCGATGGCGACCGCTGCGCGGTCGTTCGGGGTTAAACCCGCTCCTAGGGGGACCGCGCAAGGCTGTGCAGATCGTGTCGGAGGGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATAAAGCCTGGCGCTTa                                     <  1:3968897/150‑1 (MQ=255)
                                                                                                                      tcgcgcgCGATGGCGACCGCTGCGCGGACGCTCGCGGGTCAACCCGCGCCCACGGGGACCGCGCAAGCCAGTGCTCATCCTGCAGGCGCGGGTTTACCCGCGAAAGGGCCCTCACGGGCCCAAAGCGGATAAAGCCTGGCGCTTAGTTGt                                <  1:830930/144‑1 (MQ=255)
                                                                                                                      tcGTGTGCGATGGCGACCGCTGCGCGGTCGTGCGCGGGTAACCCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACGCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATAAAGCCTGGCGCTTAGTTGt                                <  2:2278908/150‑1 (MQ=255)
                                                                                                                       cGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAGCGGGCCCAAAGAGGATAAAGCCTGGCGCTTAGTTGTa                               <  2:2236749/150‑1 (MQ=255)
                                                                                                                       cGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATAAAGCCTGGCGCTTAGTTGTa                               <  2:771059/150‑1 (MQ=255)
                                                                                                                              gATGGCGACCGCTGCGCGGTCGTTGGCGGGTAAAAACCCACCAAAGGGGAACGCCCAAATCAGTGCAGATACAGTAGGAGCGGGTTTAACCGCGAAATTGCCCAAAAGGGCCCAAAAATGATAAAAACTTGCGCTTATTTTTatacagcc                        >  2:1730569/1‑150 (MQ=255)
                                                                                                                              gATGGCGACCGCTGCGCGGTCGTGCGCGGGGGCACCCCCGCCCACCGGGCCCCCGCCCGGCAGGGCCGAGCCCGTCGGCGCGGGGTTGCCCGCCCAACGGCCCAAACCGGCCCCAAGAGGAGAAACCCAGGCGCGGAGGTGTagcccccc                        >  2:2008176/1‑146 (MQ=255)
                                                                                                                                    gACCGCTGCGCGGTCGTTCGCGGGTAAAAACGCAACAAAGGGGACCCCGCAAGGCAGTGCAGAAACTGTAGGAGCGGGTTTAACCGCGAAAGGGCCCCAACGGGCCCCAACAGGATAAAGCCGGGCGCGGGGTTGTAGCCACCCCGATcc                  >  2:3673970/1‑150 (MQ=255)
                                                                                                                                       cGCTGCGCGGTCGTTCGCGGGTACACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATAAAGCCTGGCGCTTAGTTGTAGCCAGCCCGATCCATc               <  1:1013417/150‑1 (MQ=255)
                                                                                                                                       cGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATAAAGCCTGGCGCTTAGTTGTAGCCAGCCCGATCCATc               <  1:1013376/150‑1 (MQ=255)
                                                                                                                                        gcggcgcGGTCGTTCGCGGGTAAACCCGCTCCTCCGGGGACCGCGCACGGCAGTGCAGATCCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATAAAGCCTGGCGCTTAGTTGTAGCCAGCCCGATCCATCa              <  1:3440214/147‑1 (MQ=255)
                                                                                                                                         cTGCGCGGTCGTTCGCGGGGAAACCCCCTCCGACCGGGACCCCGCAAGGCAGTGCCGATACAGTAGGAGCGGGTTTACCCCCCAAAGGGCCCGAACGGGCCCAAACATGATAAAGCCAGGCGCCTAGTTGTAGCCACCCCGCTCCATCac             >  2:4410487/1‑149 (MQ=255)
                                                                                                                                             gcgGTCGTTCGCGGGTAAAACCCCACCTACGGGGACCGCGCAAGGCAGTGCAGATCCGGTAGGGGCGGGTTTACCCCCGGAAGGGCCCTAACGGGCCCCAAGAGGCTAAAGCCAGGCGCTTAGTTGTAGCCAGCCCGATCCATCAGGCGa         >  2:4270971/1‑150 (MQ=255)
                                                                                                                                                   gTTCGCGGGTAACCCCGCTCCTACGGGGACCGCGAATGGCTGTGCAGATCCTGTAGGCGCGGTTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATAAAGCCTGGCGCTTAGTTGTAGCCAGCCCGATCCATCAGGCGAATGGCt   <  1:383737/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CAGGTCGATGGATTGCCAGGACAACAGCAGCACGCTCAGGGGCAGCAGCACCAAGGCGGCGATCAGGGAGACCGGGAGGTACCAGCGGCGTTGGGGGGTGTGGGGCAAGGCTGAAATCTCGAGTGCGATGGCGACCGCTGCGCGGTCGTTCGCGGGTAAACCCGCTCCTACGGGGACCGCGCAAGGCAGTGCAGATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCTAACGGGCCCAAAGAGGATAAAGCCTGGCGCTTAGTTGTAGCCAGCCCGATCCATCAGGCGAATGGCTT  >  NC_002947/5929765‑5930062

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