Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,686,077 A→G L171P (CTG→CCG)  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,0770AG98.6% 273.8 / NA 73L171P (CTG→CCG) fumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (33/39);  minor base C (1/0);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCG  >  NC_000913/1685932‑1686223
                                                                                                                                                 |                                                                                                                                                  
caggcTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCCTCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGgtgt                                                                                                                                                >  2:19108/1‑150 (MQ=255)
  ggcTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCtg                                                                                                                                             <  2:62488/151‑1 (MQ=255)
       tATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGGTTTCTCATTCGGTGTCTGTGTc                                                                                                                                         >  2:316833/1‑150 (MQ=255)
              tatGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGtttt                                                                                                                                 >  2:905984/1‑151 (MQ=255)
               atGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGGGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTTAGGGgttt                                                                                                                                 >  2:432784/1‑150 (MQ=255)
                   ttGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGc                                                                                                                             <  2:18036/150‑1 (MQ=255)
                   ttGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCt                                                                                                                            >  1:470891/1‑151 (MQ=255)
                    tGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTg                                                                                                                           >  1:781444/1‑151 (MQ=255)
                      agaTTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTg                                                                                                                           >  2:958020/1‑149 (MQ=255)
                           aTGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAAt                                                                                                                     <  2:421628/150‑1 (MQ=255)
                             gCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGACGTGGCAACCTGCAAGTGAGTACAACCAACTTTGACGATATCGGGAAAAAATCGGGGTATCTCATTGCACGTCCTGGTCTGGGTTTTAAGCTTAAGAATg                                                                                                                    >  2:1000181/1‑149 (MQ=255)
                                cGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGtt                                                                                                                <  1:797158/150‑1 (MQ=255)
                                cGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGt                                                                                                                 <  1:523807/149‑1 (MQ=255)
                                 gAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGttgc                                                                                                              <  2:445412/151‑1 (MQ=255)
                                     aTCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGttgcttgc                                                                                                          >  2:687669/1‑151 (MQ=255)
                                        gCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTgcgc                                                                                                        >  2:655369/1‑150 (MQ=255)
                                        gCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCa                                                                                                       <  1:874149/151‑1 (MQ=255)
                                        gCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCa                                                                                                       <  2:964951/151‑1 (MQ=255)
                                         cTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTgcgc                                                                                                        >  2:636222/1‑149 (MQ=255)
                                          tACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAgc                                                                                                     >  2:850864/1‑151 (MQ=255)
                                           aCCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAgc                                                                                                     <  2:795960/150‑1 (MQ=255)
                                             ccAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGcc                                                                                                  <  1:883415/151‑1 (MQ=255)
                                                  cGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCcagcag                                                                                             <  1:150760/151‑1 (MQ=255)
                                                   ggAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCagcagc                                                                                            <  2:930970/151‑1 (MQ=255)
                                                       aTCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGGTTTAAGCTGAGGGATGAGTTGCTTGCGCAGCGCCAGCAGcgcc                                                                                         >  1:600683/1‑150 (MQ=255)
                                                        tctcCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGATGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGcgccgc                                                                                       >  2:478625/1‑151 (MQ=255)
                                                           ccTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAAc                                                                                    <  1:131030/151‑1 (MQ=255)
                                                           ccTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTCGCTTGCGCAGCGCCAGCAGCGCCGCAAc                                                                                    <  1:430917/151‑1 (MQ=255)
                                                              gCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACg                                                                                   <  1:18881/149‑1 (MQ=255)
                                                              gCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTg                                                                                 <  2:38707/151‑1 (MQ=255)
                                                                cccAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCa                                                                               <  1:819463/151‑1 (MQ=255)
                                                                 ccAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCAt                                                                              <  2:28183/151‑1 (MQ=255)
                                                                   aGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCg                                                                            <  1:1066312/151‑1 (MQ=255)
                                                                     cGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATGTCCGCCAAACCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTGTTAAGCTGAGGAATGAGTTGCTTGCGCAAAGCCAGCTGCGCCGGCACGTGCATcccc                                                                          >  2:992956/1‑151 (MQ=255)
                                                                         aaCGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCg                                                                      <  2:726164/151‑1 (MQ=255)
                                                                               gTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGa                                                                <  2:418689/151‑1 (MQ=255)
                                                                               gTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGa                                                                <  1:1044557/151‑1 (MQ=255)
                                                                                 ggCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACa                                                              <  1:811756/151‑1 (MQ=255)
                                                                                   cATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATc                                                            >  2:883937/1‑151 (MQ=255)
                                                                                    aTCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCg                                                           <  1:848115/151‑1 (MQ=255)
                                                                                      ccTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGtt                                                         <  2:459445/151‑1 (MQ=255)
                                                                                        tGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTgg                                                       >  2:249418/1‑151 (MQ=255)
                                                                                         gCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTgg                                                       <  1:642337/150‑1 (MQ=255)
                                                                                          cAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGtttg                                                       <  2:130385/149‑3 (MQ=255)
                                                                                          cAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGa                                                      >  2:645646/1‑150 (MQ=255)
                                                                                            aGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACt                                                   >  2:986154/1‑151 (MQ=255)
                                                                                             gTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACtt                                                  <  1:684803/151‑1 (MQ=255)
                                                                                              tGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACttt                                                 >  1:513208/1‑151 (MQ=255)
                                                                                               gAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTg                                                <  1:1031298/151‑1 (MQ=255)
                                                                                                aGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTgg                                               >  1:22723/1‑151 (MQ=255)
                                                                                                 gTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGc                                              <  1:527440/151‑1 (MQ=255)
                                                                                                  tACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCt                                             >  2:545666/1‑151 (MQ=255)
                                                                                                   aCGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCt                                             <  1:659117/150‑1 (MQ=255)
                                                                                                           ttttGACGATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCAc                                    >  1:65019/1‑151 (MQ=255)
                                                                                                                aCGATATCGGCAAAAGCACGGGGTTTCTCATTCGGTGTCTGTGTCAGGGGTTTAAGCTGAGGAATGAGTTGCTTGCGAAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCAcgtcgt                               >  2:519905/1‑151 (MQ=255)
                                                                                                                  gATATCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCAcgtcgtc                              >  2:971339/1‑150 (MQ=255)
                                                                                                                   atatCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCAcgtcgtcgt                            <  1:38874/151‑1 (MQ=255)
                                                                                                                   atatCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCAcgtcgtcgt                            >  1:340003/1‑151 (MQ=255)
                                                                                                                     atCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTa                          <  2:474336/151‑1 (MQ=255)
                                                                                                                      tCGGCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGtt                           >  1:689860/1‑149 (MQ=255)
                                                                                                                         gCAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGt                      <  1:1010531/151‑1 (MQ=255)
                                                                                                                          cAAAAGCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTg                     >  2:855197/1‑151 (MQ=255)
                                                                                                                               gCACGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACtt                <  1:252558/151‑1 (MQ=255)
                                                                                                                                  cGGGATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTAc             <  1:791492/151‑1 (MQ=255)
                                                                                                                                  cGGGATTTCTCATTCGGTGTCTGTGTAAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGATAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTa              >  2:408506/1‑150 (MQ=255)
                                                                                                                                    ggATTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGt           >  2:737596/1‑151 (MQ=255)
                                                                                                                                      aTTTCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTc         <  1:520256/151‑1 (MQ=255)
                                                                                                                                        ttCTCATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTcc        <  2:190938/150‑1 (MQ=255)
                                                                                                                                         tctcATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCa       <  2:493739/150‑1 (MQ=255)
                                                                                                                                         tctcATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCAt      >  1:363688/1‑151 (MQ=255)
                                                                                                                                            cATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATccc   >  2:946492/1‑151 (MQ=255)
                                                                                                                                            cATTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATccc   >  1:414150/1‑151 (MQ=255)
                                                                                                                                             aTTCGGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCcg  >  1:379763/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCCCAGCGTTAACGGCGTGGCATCCTGCAAGTGAGTACGACCAATTTTGACGATATCGGCAAAAGCACGGGATTTCTCATTCAGTGTCTGTGTCAGGGTTTTAAGCTGAGGAATGAGTTGCTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCG  >  NC_000913/1685932‑1686223

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