Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,686,044 A→C 7.2% V182G (GTC→GGC)  fumC ← fumarate hydratase (fumarase C),aerobic Class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,686,0440AC 7.2% 177.8 / 5.1 66V182G (GTC→GGC) fumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base C (2/3);  total (31/35)
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.74e-01

GCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTT  >  NC_000913/1685897‑1686186
                                                                                                                                                   |                                                                                                                                              
gCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGAc                                                                                                                                               <  1:579954/149‑1 (MQ=255)
gCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGAc                                                                                                                                               >  2:1343878/1‑149 (MQ=255)
         ccAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGc                                                                                                                                      <  1:1534818/149‑1 (MQ=255)
                 cAGTTCCGCTACGTGATGCAGTCTGTATTCGATATGTTTGATATTATGCTCGAGCATCGCTAACCAGCCTGAAATCTCCTGCCCCAGCGTTAACGGCGGGGAATCCTGCAAGTGAGTACGACCAATTTTGCCGATATCGTCaaaattgg                                                                                                                              <  2:548221/149‑5 (MQ=255)
                  aGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACg                                                                                                                             <  1:654512/149‑1 (MQ=255)
                   gTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACgg                                                                                                                            >  1:1022656/1‑149 (MQ=255)
                   gTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACgg                                                                                                                            >  1:919688/1‑149 (MQ=255)
                    ttCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCCACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGCCGATATCGTCaaaattggtcg                                                                                                                           <  2:374897/149‑8 (MQ=255)
                        gCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGAttt                                                                                                                       <  1:418782/149‑1 (MQ=255)
                           aCGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTtctc                                                                                                                    <  2:1506833/149‑1 (MQ=255)
                            cGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCa                                                                                                                   <  2:945436/149‑1 (MQ=255)
                             gTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCAt                                                                                                                  >  1:1310791/1‑149 (MQ=255)
                               gAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTc                                                                                                                >  2:478442/1‑149 (MQ=255)
                               gAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAAGTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTc                                                                                                                >  2:1142480/1‑149 (MQ=255)
                                  gcaggcTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAgt                                                                                                             <  2:27239/149‑1 (MQ=255)
                                        tGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCtgt                                                                                                       >  1:690520/1‑149 (MQ=255)
                                        tGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCtgt                                                                                                       >  2:939835/1‑149 (MQ=255)
                                            ttCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCCTTCAGTGTCTGTGTCa                                                                                                   >  2:264537/1‑149 (MQ=255)
                                               gATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAggg                                                                                                <  2:428609/149‑1 (MQ=255)
                                               gATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAggg                                                                                                <  1:111139/149‑1 (MQ=255)
                                                  atGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGttt                                                                                             <  2:1310791/149‑1 (MQ=255)
                                                     tttGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTaa                                                                                          >  2:923976/1‑149 (MQ=255)
                                                      ttGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCCATTTTGGCGATATCGGCAAAAGCACGGGATTTCTCCTTCAGTGTCTGTGTCAGGGTTTTAAg                                                                                         >  1:1011258/1‑149 (MQ=255)
                                                          taTTATGCTCGAGCATCTCTACCCAGACTGAAATCTACTGCCACAGCTTTAACGGCGTGGCATCCTTCAAGTGATTACGACCAATTTTGAAGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTgt                                                                                     <  1:1205148/148‑2 (MQ=255)
                                                           tttATGGTCGAGCATCGCTACCCAGCCGGAAATATCCTGCCCAATCGTTAACGGCGTGGCATCCTGCAAGTTAGTAAGACCAATTTTGACGATATCGGCAAAAGCACTGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAg                                                                                    <  1:1166628/148‑1 (MQ=255)
                                                           aTTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCGCATTCAGTGTCTGTGTCAGGGTGTTAAGCTGAg                                                                                    >  1:785630/1‑149 (MQ=255)
                                                               tGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAAt                                                                                <  1:1466356/149‑1 (MQ=255)
                                                                gCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATg                                                                               >  2:1249949/1‑149 (MQ=255)
                                                                gCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATg                                                                               >  1:1018165/1‑149 (MQ=255)
                                                                 cTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGa                                                                              <  2:1541043/149‑1 (MQ=255)
                                                                  tCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAg                                                                             >  2:28846/1‑149 (MQ=255)
                                                                  tCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAg                                                                             <  1:28846/149‑1 (MQ=255)
                                                                    gAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGtt                                                                           >  2:40575/1‑149 (MQ=255)
                                                                    gAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGtt                                                                           >  1:271235/1‑149 (MQ=255)
                                                                     aGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCTTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGttg                                                                          <  2:1110530/149‑1 (MQ=255)
                                                                              aCCCAGCCGGAAATCTCCTGCCCCAGCGTTTACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAg                                                                 <  1:947977/149‑1 (MQ=255)
                                                                               cccAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAgc                                                                <  2:585870/149‑1 (MQ=255)
                                                                                 cAGCCGGAAATCTCCTGCCCCATCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAgcgc                                                              <  1:948503/149‑1 (MQ=255)
                                                                                    ccGGAAATCTCCTGCCCCAGCCTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGCGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCcag                                                           >  2:1210874/1‑149 (MQ=255)
                                                                                      ggAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCcagca                                                         >  2:1430003/1‑149 (MQ=255)
                                                                                             tcCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGCCGATATCGTCaaaat                                                                                                                                 <  1:625871/70‑2 (MQ=39)
                                                                                             tcCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGCCGATATCGTCaaaat                                                                                                                                 >  2:625871/1‑69 (MQ=39)
                                                                                             tcCTGCCCCAGCGATAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGCCGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTGTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGcgccgc                                                  >  1:269435/1‑149 (MQ=255)
                                                                                                  ccccAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGt                                             <  2:555217/149‑1 (MQ=255)
                                                                                                  ccccAGCGGTAACGGCGGGGCATCCTGCAAGTGAGGACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCCGGGTCAGGGTGTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGt                                             <  1:264537/149‑1 (MQ=255)
                                                                                                    ccAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGc                                           >  1:350300/1‑149 (MQ=255)
                                                                                                      cgCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCAt                                         <  1:1343878/148‑1 (MQ=255)
                                                                                                        cGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCg                                       <  1:177205/149‑1 (MQ=255)
                                                                                                           tAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGAGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCg                                    >  1:1445830/1‑149 (MQ=255)
                                                                                                              cGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCg                                 >  1:93953/1‑149 (MQ=255)
                                                                                                                 cGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGaa                              >  2:113301/1‑149 (MQ=255)
                                                                                                                  gTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCTGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGaaa                             <  1:1210874/149‑1 (MQ=255)
                                                                                                                  gTGGCATCCTGCAAGTGAGTACGACCAATATTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGGCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGACAGCAGCGCCGCAACGTGCATCGCCGTCGGaaa                             <  1:627421/149‑1 (MQ=255)
                                                                                                                    ggCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGa                           <  2:620888/149‑1 (MQ=255)
                                                                                                                               aGTGAGTACGACCAAATTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCCGTGGCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGaa                <  2:1015122/149‑1 (MQ=255)
                                                                                                                                 tGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACt              >  2:154263/1‑149 (MQ=255)
                                                                                                                                 tGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACt              <  2:470584/149‑1 (MQ=255)
                                                                                                                                         aCCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGc         >  2:1176520/1‑146 (MQ=255)
                                                                                                                                         aCCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGc         <  1:1176520/146‑1 (MQ=255)
                                                                                                                                         aCCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCttt      >  1:1517888/1‑149 (MQ=255)
                                                                                                                                         aCCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCttt      <  2:350300/149‑1 (MQ=255)
                                                                                                                                         aCCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCttt      >  2:1548751/1‑149 (MQ=255)
                                                                                                                                           cAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTg    <  1:1249949/149‑1 (MQ=255)
                                                                                                                                           cAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTg    >  2:1276994/1‑149 (MQ=255)
                                                                                                                                           cAATTTTGACGATATCGGCAAAAGAACGGGATTTCTCATTCAGAGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGAAGCGCCAGCAGCGCCGCAACTTGCATCGCCGTCGGAAAGACATCGTTGGCACTTTCGCTTTTg    >  1:591203/1‑149 (MQ=255)
                                                                                                                                             aTTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGtt  <  1:1216028/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTT  >  NC_000913/1685897‑1686186

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: