Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,260,232 C→T 6.4% L393L (CTG→CTA)  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,2320CT 6.4% 192.2 / 2.9 62L393L (CTG→CTA) fruAfused fructose‑specific PTS enzymes: IIBcomponent/IIC components
Reads supporting (aligned to +/- strand):  ref base C (32/26);  new base T (2/2);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GCGGTGGCACCATACCTGCCGCCATAATCGCCGCCATCGGGCCATAGGTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCGGTTTCAGC  >  NC_000913/2260087‑2260376
                                                                                                                                                 |                                                                                                                                                
gCGGTGGCACCATACCTGCCGCCATAATCGCCGCCTTCGGGCCATAGGTTTGAGTACTCAGCAGACCCACCCCGAAGGCGTACGGTGCTTTGTTTACCGGACCGCCCATGACAGTACACAACTTGCCACCGAGGATCGCCCCccgcgg                                                                                                                                                >  2:341858/1‑146 (MQ=255)
        aCCATACCTGCCGCCATAATCGCCGCCATCGGGCCATAGGTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCAtt                                                                                                                                       >  1:269485/1‑149 (MQ=255)
          cATACCTGCCGCCATAATCGCCGCCATCGGGCCATAGGTTTGAGTACTCAGCAGCCCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCg                                                                                                                                     <  3:142439/149‑1 (MQ=255)
            tACCTGCCGCCATAATCGCCGCCATCGGGCCATAGGTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCa                                                                                                                                   <  3:263274/149‑1 (MQ=255)
             aCCTGCCGCCATAATCGCCGCCATCGGGCCATAGGTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCa                                                                                                                                         <  2:321678/142‑1 (MQ=255)
             aCCTGCCGCCATAATCGCCGCCATCGGGCCATAGGTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCa                                                                                                                                         >  1:321678/1‑142 (MQ=255)
               cTGCCGCCAAAATCGCCGCCATCGGGCCATAGGTATGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGATAACCGGACCGCCCATGTCAGTACACATCATGCCACCTAGGATCGCCCCCAGCAGAACCGCATTAGCAGta                                                                                                                                >  2:138189/1‑148 (MQ=255)
                   cgccATAATCGCCGCCATCGGGCCATAGGTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCa                                                                                                                            <  3:139368/149‑1 (MQ=255)
                    gccATAATCGCCGCCATCGGGCCATAGGTTTGAGTACTCATCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCAt                                                                                                                           <  1:153515/149‑1 (MQ=255)
                              ccgccATCGGGCCATAGGTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGc                                                                                                                 <  7:322952/149‑1 (MQ=255)
                              ccgccATCGGGCCATAGGTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGc                                                                                                                 >  1:49288/1‑149 (MQ=255)
                                              ggTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCtt                                                                                                 <  2:49288/149‑1 (MQ=255)
                                               gTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTc                                                                                                <  1:26027/149‑1 (MQ=255)
                                                 ttGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTc                                                                                                >  8:290108/1‑147 (MQ=255)
                                                  tGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCga                                                                                              <  1:341858/148‑1 (MQ=255)
                                                    aGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCgaga                                                                                            <  2:20165/148‑1 (MQ=255)
                                                     gTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGTAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAAt                                                                                          >  1:153223/1‑149 (MQ=255)
                                                     gTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGTAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAAt                                                                                          >  1:148118/1‑149 (MQ=255)
                                                     gTACTCAGCAGACCCACACCGAAGGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAAt                                                                                          <  6:193264/149‑1 (MQ=255)
                                                      tACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATg                                                                                         >  1:131369/1‑149 (MQ=255)
                                                       aCTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGc                                                                                        >  5:121812/1‑149 (MQ=255)
                                                         tCAGCAGACCCACACAGAATGCGTACGCGGCTTTGTTTACCGGACCGCCCATGTCAGTACCCATCATGCCCCCGAGGATAGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCAGCAGCCCGTGAGTCAGCCCTTCGAGAATGc                                                                                        >  3:147730/1‑147 (MQ=255)
                                                           agcagACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCa                                                                                      <  7:290108/147‑1 (MQ=255)
                                                           agcagACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAg                                                                                     >  4:315638/1‑148 (MQ=255)
                                                            gcagACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGc                                                                                    >  8:202121/1‑148 (MQ=255)
                                                                    acacCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGtt                                                                            >  1:89188/1‑148 (MQ=255)
                                                                      acCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTAc                                                                         <  6:121812/149‑1 (MQ=255)
                                                                        cGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTAcc                                                                        >  1:4186/1‑148 (MQ=255)
                                                                           aTGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGAt                                                                     >  1:329043/1‑148 (MQ=255)
                                                                            tGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCa                                                                   <  2:89188/149‑1 (MQ=255)
                                                                              cGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGTAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAgg                                                                 <  2:153223/149‑1 (MQ=255)
                                                                               gTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGTAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAgg                                                                 <  2:148118/148‑1 (MQ=255)
                                                                                tACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCa                                                                                                                                   >  2:106008/1‑81 (MQ=255)
                                                                                tACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCa                                                                                                                                   <  1:106008/81‑1 (MQ=255)
                                                                                tACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATAGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAgg                                                                 >  1:178418/1‑147 (MQ=255)
                                                                                   gctgctTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGat                                                            <  2:329043/149‑1 (MQ=255)
                                                                                     tgctTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGatca                                                          >  5:146010/1‑149 (MQ=255)
                                                                                            tttACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGc                                                      >  6:337651/1‑146 (MQ=255)
                                                                                            tttACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCa                                                    <  7:47050/148‑1 (MQ=255)
                                                                                              tACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGCATGCCAGCAACTGGGTTACCGATCAGGTAGATCATCGCCAgcc                                                 >  3:219423/1‑147 (MQ=255)
                                                                                                ccGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAgacc                                                <  7:202121/148‑1 (MQ=255)
                                                                                                ccGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCAAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAgaccg                                               >  6:166214/1‑149 (MQ=255)
                                                                                                 cGGACCGCCCATGTCAGTACACATCATGCCCCCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCACCTGGTTTACCGATCAGGTAGATCATCGCCAgaccga                                              >  5:117726/1‑149 (MQ=255)
                                                                                                 cGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAgaccga                                              <  3:315638/149‑1 (MQ=255)
                                                                                                      cGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGaccac                                          >  3:11885/1‑148 (MQ=255)
                                                                                                                aGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCAAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGaa                                <  5:166214/148‑1 (MQ=255)
                                                                                                                      catcatGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCACAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTTGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCg                         >  3:50161/1‑149 (MQ=255)
                                                                                                                               aCCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCa                <  4:11885/149‑1 (MQ=255)
                                                                                                                               aCCGAGGAGCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGGCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCa                <  2:269485/149‑1 (MQ=255)
                                                                                                                                ccGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGACTGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGCCTGGAAATTAGCGGAATGATCAg               >  6:1195/1‑149 (MQ=255)
                                                                                                                                 cGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAg               >  4:44040/1‑148 (MQ=255)
                                                                                                                                 cGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACAACCAGACTGGAAATTTGTGGAATGATCAg               >  1:199985/1‑148 (MQ=255)
                                                                                                                                   tggATCGCCCCCAGCAGATCGGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTGAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTATATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGAt            <  1:138189/148‑1 (MQ=255)
                                                                                                                                   aGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGGAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGAt            >  1:132777/1‑149 (MQ=255)
                                                                                                                                   aGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGAt            <  2:178418/149‑1 (MQ=255)
                                                                                                                                   aGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGAt            >  8:280689/1‑149 (MQ=255)
                                                                                                                                       tCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGCCCACCAGACTGGAAATTAGCGGAATGATCAGGATCGGt        >  1:207860/1‑149 (MQ=255)
                                                                                                                                         gCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCGGttt      <  2:132777/149‑1 (MQ=255)
                                                                                                                                           ccccAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTg                                                                               >  4:131022/1‑74 (MQ=255)
                                                                                                                                           ccccAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTg                                                                               <  3:131022/74‑1 (MQ=255)
                                                                                                                                           ccccAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGCCCGCCCACCCGACTGGAAATTAGCGGAATGCTCAGGATCGGTTTCa    >  1:297166/1‑149 (MQ=255)
                                                                                                                                            cccAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGAAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCGGTTTCAg   >  2:95658/1‑149 (MQ=255)
                                                                                                                                             ccAGCAGAACCGCATTCGCAGTCCCCATCGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCGGTTTCAgc  >  6:43406/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCGGTGGCACCATACCTGCCGCCATAATCGCCGCCATCGGGCCATAGGTTTGAGTACTCAGCAGACCCACACCGAATGCGTACGCTGCTTTGTTTACCGGACCGCCCATGTCAGTACACATCATGCCACCGAGGATCGCCCCCAGCAGAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCGGTTTCAGC  >  NC_000913/2260087‑2260376

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