Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,074,135 G→A 6.7% V133M (GTG→ATG)  yihV → 6‑deoxy‑6‑sulphofructose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,074,1350GA 6.7% 187.5 / 2.9 61V133M (GTG→ATG) yihV6‑deoxy‑6‑sulphofructose kinase
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base A (2/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

ATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCA  >  NC_000913/4073998‑4074281
                                                                                                                                         |                                                                                                                                                  
ataACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGCGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGAt                                                                                                                                         >  2:87069/1‑149 (MQ=255)
 taACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGAt                                                                                                                                         <  1:145900/148‑1 (MQ=255)
   aCCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTa                                                                                                                                      >  2:98705/1‑149 (MQ=255)
    ccAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGCGCGTATAATCATGACCTACCCACGCCCGGATCTGCTGCATGACGGAGAGTGGGTGGAGGAAATTGAGTTCGAGCTGTGGGAGGGTGTGCTGGCAGagggag                                                                                                                                     >  7:173451/1‑144 (MQ=255)
     cAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTATGCTGGCAGATGTACg                                                                                                                                    >  5:10883/1‑149 (MQ=255)
     cAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTATGCTGGCAGATGTACg                                                                                                                                    >  5:14143/1‑149 (MQ=255)
      aGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGc                                                                                                                                   <  1:89348/149‑1 (MQ=255)
          gAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGc                                                                                                                               >  8:36013/1‑149 (MQ=255)
              tCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCacga                                                                                                                           >  7:209187/1‑149 (MQ=255)
                   gCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGgcg                                                                                                                      <  1:183934/149‑1 (MQ=255)
                   gCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGgcg                                                                                                                      >  3:64794/1‑149 (MQ=255)
                    cAATCCGCCATCATGGTGGATACCAAATGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGTGATGTTGTGCTGGCAGATGTACGCTGGCACGACGgcgc                                                                                                                     >  4:110040/1‑149 (MQ=255)
                     aaTCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCCGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAGTTCTCTCTGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGgcgc                                                                                                                     >  1:130714/1‑148 (MQ=255)
                      aTCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCt                                                                                                                    >  1:203374/1‑148 (MQ=255)
                           ccATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTaaaaaa                                                                                                              <  8:84761/149‑1 (MQ=255)
                            catcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAg                                                                                                             >  7:148644/1‑149 (MQ=255)
                                         aCCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGc                                                                                                >  4:229245/1‑149 (MQ=255)
                                           cAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGccc                                                                                              >  2:210119/1‑149 (MQ=255)
                                                 cGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAgg                                                                                        >  8:147925/1‑149 (MQ=255)
                                                  gAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGc                                                                                       <  5:186681/149‑1 (MQ=255)
                                                  gAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGc                                                                                       >  3:39590/1‑149 (MQ=255)
                                                    gCGGATTATCATTAACTACCCCATCCCGTATCTGCTTCCTGACTCAGAGTGGTTGGAGGAAATTCATTTCTCTCACTGGGAAGTTGTGCTGGCAGAGGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTCGGCCGTCAGGCgg                                                                                     <  3:197422/149‑1 (MQ=255)
                                                    gCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCgg                                                                                     <  1:215979/149‑1 (MQ=255)
                                                      cgATAATCATTACCTACCCCCGCCCGGAGCTGCTGCCTGACGCAGAGTGGCTGGACGAAATTGATACCTCTCAGTGGGATGGTGTGCTGGCAGATGGACGCTGGCACGACGGCGCTAAAACAGCCCTCACCCTGGCCCGTCAGGCgggg                                                                                   <  4:39590/148‑2 (MQ=255)
                                                       gATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGgtg                                                                                  <  1:98705/149‑1 (MQ=255)
                                                           aTCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTGGTGCTGGCAGATGGACGCTGGCACGACGGCGCTAAAACAGCCTTCACCCTGGCCCGTCAGGCGGGTGtgat                                                                              <  5:211271/149‑1 (MQ=255)
                                                             caatAACTTCCCCAGGCCGGATCTGCTGCCTGACGGAGTGTGGTTGGAGGTAATTGATTTCTCTGAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGTCGGCGCTAAAAAAGcc                                                                                                           <  4:192061/115‑1 (MQ=255)
                                                             cATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGcc                                                                                                           >  3:192061/1‑118 (MQ=255)
                                                              aTTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGAc                                                                           >  8:89289/1‑149 (MQ=255)
                                                              aTTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGAc                                                                           >  5:56269/1‑149 (MQ=255)
                                                                   cTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTc                                                                      <  8:148644/149‑1 (MQ=255)
                                                                       cccAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGa                                                                  >  7:12685/1‑149 (MQ=255)
                                                                           gCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACggg                                                              <  1:11940/149‑1 (MQ=255)
                                                                            cccGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTCTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGAACGTTCTGGACggg                                                              >  1:108990/1‑148 (MQ=255)
                                                                             ccGGATCTGCTGCGTGTCGCAGAGTGGTTGGAGGAAAGTGATTTATCTCTGTGTGATGGTGTGCTGGCAGATGTACGATGGCACGACGGGGCTATAAAAGCCTTCACCCTGGCCCTTCAGGCGGGTGTGATGAACGTTCTGGAAGGGGa                                                            <  1:107399/149‑1 (MQ=255)
                                                                                  tCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGATGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTa                                                       >  2:222372/1‑149 (MQ=255)
                                                                                      ctgcCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACAtt                                                        <  4:64794/144‑1 (MQ=255)
                                                                                      ctgcCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTATGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTAcgcc                                                   <  6:14143/149‑1 (MQ=255)
                                                                                      ctgcCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTATGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTAcgcc                                                   <  6:10883/149‑1 (MQ=255)
                                                                                          cTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAg                                               <  7:36013/149‑1 (MQ=255)
                                                                                           tGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGtt                                                                       <  8:24834/124‑1 (MQ=255)
                                                                                           tGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGtt                                                                       >  7:24834/1‑124 (MQ=255)
                                                                                                       gTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCa                                        >  2:207793/1‑143 (MQ=255)
                                                                                                       gTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCa                                        <  1:207793/143‑1 (MQ=255)
                                                                                                        ttGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCt                                 >  8:152843/1‑149 (MQ=255)
                                                                                                           gaggaAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCtggt                              >  1:128076/1‑149 (MQ=255)
                                                                                                           gaggaAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCtgg                               <  3:229245/148‑1 (MQ=255)
                                                                                                              gaAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGc                           <  1:210119/149‑1 (MQ=255)
                                                                                                                aaTTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCAt                         >  5:107847/1‑149 (MQ=255)
                                                                                                                     aTTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGAc                                                           <  3:153445/110‑1 (MQ=255)
                                                                                                                     aTTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGAc                                                           >  4:153445/1‑110 (MQ=255)
                                                                                                                     aTTTCTCTCAGTGGGAGGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCATTAAg                     <  8:44592/148‑1 (MQ=255)
                                                                                                                        tctctcAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGa                  <  7:217787/148‑1 (MQ=255)
                                                                                                                         ctctcAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGAt                 <  7:152843/148‑1 (MQ=255)
                                                                                                                          tctcAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGgt                                                                                   <  3:88228/81‑1 (MQ=255)
                                                                                                                          tctcAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGgt                                                                                   >  4:88228/1‑81 (MQ=255)
                                                                                                                            tcAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCAcg             >  7:21475/1‑149 (MQ=255)
                                                                                                                            tcAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCAcg             >  4:32589/1‑149 (MQ=255)
                                                                                                                               gTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCgg          >  4:111708/1‑149 (MQ=255)
                                                                                                                                tGGGATGTTGTGCTGGCAGATGTACGCTGGCACGGCGGCGCTAAAAAAGCCTTCACACTGGCCCGTCAGGAGGGTGTGATGACCGTTCTGGACGGGGAAATTACGCCGCAGGATATCAGTGAGCGGGTGGCATTAAGCGATCACGCTGc         >  3:196717/1‑149 (MQ=255)
                                                                                                                                tGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGc         >  6:122426/1‑149 (MQ=255)
                                                                                                                                       ttgtGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCa  <  1:87069/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
ATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCCTGGCCCGTCAGGCGGGTGTGATGACCGTTCTGGACGGGGACATTACGCCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCA  >  NC_000913/4073998‑4074281

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: ‑