Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,260,769 C→T 6.7% T214T (ACG→ACA)  fruA ← fused fructose‑specific PTS enzymes: IIBcomponent/IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,260,7690CT 6.7% 186.8 / 3.6 60T214T (ACG→ACA) fruAfused fructose‑specific PTS enzymes: IIBcomponent/IIC components
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base T (2/2);  total (28/32)
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.10e-01

TACCCGGCTCTTTAAACGCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCT  >  NC_000913/2260621‑2260913
                                                                                                                                                    |                                                                                                                                                
tACCCGGCTCTTTAAACGCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGcc                                                                                                                                                  <  2:471531/149‑1 (MQ=255)
    cGGCTCTTTAAACGCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGttt                                                                                                                                              <  3:409972/149‑1 (MQ=255)
    cGGCTCTTTAAACGCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGttt                                                                                                                                              >  2:399669/1‑149 (MQ=255)
     ggCTCTTTAAACGCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTg                                                                                                                                             >  2:192663/1‑149 (MQ=255)
      gCTCTTTAACCGCTTCGATACCAAAAGCAAAAGAAAGCCCGATCCAGAGACCGCCTGCAACGACCATCGGTAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGTGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGa                                                                                                                                            <  3:30261/149‑1 (MQ=255)
      gCTCTTTAAACGCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGa                                                                                                                                            <  1:159686/149‑1 (MQ=255)
      gCTCTTTAAACGCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGa                                                                                                                                            <  1:4878/149‑1 (MQ=255)
        tctTTAAACGCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGc                                                                                                                                          >  2:33295/1‑149 (MQ=255)
                 gCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACATACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCccgc                                                                                                                                 <  3:381096/149‑1 (MQ=255)
                 gCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGAGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCccgc                                                                                                                                 >  4:289835/1‑149 (MQ=255)
                  cTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCccgcc                                                                                                                                >  4:80044/1‑149 (MQ=255)
                          ccAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATa                                                                                                                        >  4:548631/1‑149 (MQ=255)
                           caaaagcaaaagAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATAc                                                                                                                       >  1:278054/1‑149 (MQ=255)
                            aaaagcaaaagAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACg                                                                                                                      <  4:274570/149‑1 (MQ=255)
                             aaagcaaaagAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACgg                                                                                                                     >  4:145728/1‑149 (MQ=255)
                                 caaaagAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGtt                                                                                                                 >  4:212947/1‑149 (MQ=255)
                                  aaaagAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTg                                                                                                                >  3:19642/1‑149 (MQ=255)
                                    aagaaagCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATTTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCt                                                                                                              <  2:26909/149‑1 (MQ=255)
                                     agaaagCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCtt                                                                                                             <  3:219346/149‑1 (MQ=255)
                                      aaaacCGCGATACACAGGCCGCCTACAAGGGCCATCGGCCGCATGTAAGAAACGCCCGACAGGAAGTGACTGTAGGCGCCTGCACCCTCTTTCTTAGTTTCAGTGGTCGCCGTTTGACTTTTGCCCGCCGGTTCATACGGCGTTGCTTc                                                                                                            <  4:27097/148‑1 (MQ=255)
                                            gcgATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAAc                                                                                                      >  3:326335/1‑149 (MQ=255)
                                                   acaGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTa                                                                                               >  3:373958/1‑149 (MQ=255)
                                                   acaGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTa                                                                                               >  1:396646/1‑149 (MQ=255)
                                                     aGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATc                                                                                             <  3:527025/149‑1 (MQ=255)
                                                            cTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCTGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTcc                                                                                      >  1:448483/1‑149 (MQ=255)
                                                            cTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCTGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTcc                                                                                      >  1:482720/1‑149 (MQ=255)
                                                            cTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCa                                                                                                           >  4:285160/1‑128 (MQ=255)
                                                            cTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCa                                                                                                           <  3:285160/128‑1 (MQ=255)
                                                             tGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCTGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCt                                                                                     <  2:448483/149‑1 (MQ=255)
                                                             tGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCTGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCt                                                                                     <  2:482720/149‑1 (MQ=255)
                                                                 aCGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTgcgc                                                                                 <  4:373958/149‑1 (MQ=255)
                                                                   gACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCgg                                                                               <  2:396646/149‑1 (MQ=255)
                                                                     ccATCGGCAGCATGTAAGAACCGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGtt                                                                             <  1:422465/149‑1 (MQ=255)
                                                                     ccAACGGCAGCATGTAAGAAAGGCCCGTCAGCAAGTGCCGGTATGCGCCTGCACCCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATGCGGCGTTGCTTGAGCAACCGCTTTATCCAGTTCCTGCGCGGtt                                                                             <  1:578255/149‑1 (MQ=255)
                                                                      cATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGttt                                                                            >  1:73415/1‑149 (MQ=255)
                                                                        tCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTttc                                                                          <  1:72259/149‑1 (MQ=255)
                                                                           gcagcaTGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTttcttc                                                                       >  2:147211/1‑149 (MQ=255)
                                                                            cagcaTGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCa                                                                      <  1:398302/149‑1 (MQ=255)
                                                                                    aaGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGa                                                              <  2:73415/149‑1 (MQ=255)
                                                                                     aGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGCTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGAc                                                             <  4:332611/149‑1 (MQ=255)
                                                                                      gAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAg                                                                                                          <  3:327037/103‑1 (MQ=255)
                                                                                      gAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAg                                                                                                          >  4:327037/1‑103 (MQ=255)
                                                                                            cccGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTaga                                                      <  1:204410/149‑1 (MQ=255)
                                                                                                   gCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGa                                               <  3:369012/149‑1 (MQ=255)
                                                                                                     aaGTGATGGTATTCGCCCGCACTCTCTTTCTTACTTTAAGTGGTCGCCGTTTGAGATTTGCCCGCCGGTTCATACCGCATTGCTTCAGCAACCGCTCTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATa                                             >  4:510751/1‑149 (MQ=255)
                                                                                                      aGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATAc                                            <  1:399669/149‑1 (MQ=255)
                                                                                                       gTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACa                                           <  2:278054/149‑1 (MQ=255)
                                                                                                              tATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGttt                                    >  3:427586/1‑149 (MQ=255)
                                                                                                                 gcgcCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTAcc                                 <  3:433417/149‑1 (MQ=255)
                                                                                                                     cTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCa                             >  4:188458/1‑149 (MQ=255)
                                                                                                                       gCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCaaa                           <  1:147211/149‑1 (MQ=255)
                                                                                                                                   ttACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTAt                                                                                              >  1:195211/1‑70 (MQ=255)
                                                                                                                                   ttACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTAt                                                                                              <  2:195211/70‑1 (MQ=255)
                                                                                                                                   ttACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGAt                <  1:390772/148‑1 (MQ=255)
                                                                                                                                   ttACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGAt                >  2:390772/1‑148 (MQ=255)
                                                                                                                                   ttACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATc               >  2:460552/1‑149 (MQ=255)
                                                                                                                                   ttACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATc               >  4:552439/1‑149 (MQ=255)
                                                                                                                                       tttCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACt           <  1:63709/149‑1 (MQ=255)
                                                                                                                                       tttCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACt           >  4:410656/1‑149 (MQ=255)
                                                                                                                                             tGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGata     >  4:409373/1‑149 (MQ=255)
                                                                                                                                                tCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCt  <  3:212947/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TACCCGGCTCTTTAAACGCTTCGATACCAAAAGCAAAAGAAAGCGCGATACACAGACCACCTGCAACGACCATCGGCAGCATGTAAGAAACGCCCGTCAGCAAGTGACGGTATGCGCCTGCACTCTCTTTCTTACTTTCAGTGGTCGCCGTTTGAGCTTTGCCCGCCGGTTCATACGGCGTTGCTTCAGCAACCGCTTTATCCAGTTCCTGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCT  >  NC_000913/2260621‑2260913

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