Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,838,680 A→G 6.0% I195I (ATT→ATC)  ascG ← asc operon transcriptional repressor; prpBC operon repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,838,6800AG 6.0% 209.9 / 3.1 67I195I (ATT→ATC) ascGasc operon transcriptional repressor; prpBC operon repressor
Reads supporting (aligned to +/- strand):  ref base A (28/35);  new base G (2/2);  total (30/37)
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.86e-01

CTAACGCGCTAAATTTAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGGCTGCTGTTTTTGCGCAGGC  >  NC_000913/2838532‑2838816
                                                                                                                                                    |                                                                                                                                        
cTAACGCGCTAAATTTAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCa                                                                                                                                          <  3:255831/149‑1 (MQ=255)
   aCGCGCTAAATTTAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCaa                                                                                                                                         <  6:332357/147‑1 (MQ=255)
    cgcgCTAAATTTAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGc                                                                                                                                      <  4:177192/149‑1 (MQ=255)
      cgcTAAATTTCGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCtg                                                                                                                                    <  1:271447/149‑1 (MQ=255)
         tAAATTTAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCtggtg                                                                                                                                 >  2:57156/1‑149 (MQ=255)
           aaTTTAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTggg                                                                                                                               <  4:102513/149‑1 (MQ=255)
           aaTTTAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTggg                                                                                                                               <  4:90888/149‑1 (MQ=255)
           aaTTTAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTACATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCtggtgg                                                                                                                                <  2:30161/148‑1 (MQ=255)
            aTTTAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTgggg                                                                                                                              <  4:115644/149‑1 (MQ=255)
              ttAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTggggg                                                                                                                             >  1:364654/1‑148 (MQ=255)
               tAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGa                                                                                                                            <  3:96687/148‑1 (MQ=255)
                    ccACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCAt                                                                                                                      >  1:267175/1‑149 (MQ=255)
                    ccACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCAt                                                                                                                      <  5:177442/149‑1 (MQ=255)
                     cACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCTGCAAGACGTTCAATGCTGGTGGGGGAATCCATc                                                                                                                     <  2:236163/149‑1 (MQ=255)
                     cACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATc                                                                                                                     >  8:179457/1‑149 (MQ=255)
                     cACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATc                                                                                                                     <  2:364654/149‑1 (MQ=255)
                              gCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATAGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATc                                                                                                                     >  8:94238/1‑140 (MQ=255)
                              gCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATAGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATc                                                                                                                     <  7:94238/140‑1 (MQ=255)
                                aaCATTTCTACCCATTCGGACCCGCTCGCAGGCGTCCATTTACCGTTATCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACACGGCAAGCCGTTCAATGCTGGTGGTGGAATCCATCGAGCCGGTAAg                                                                                                          >  8:63041/1‑149 (MQ=255)
                                 aCATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAg                                                                                                          >  4:170224/1‑148 (MQ=255)
                                        tACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCaa                                                                                                  >  7:175547/1‑149 (MQ=255)
                                         aCCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCGATGGTGGTgggg                                                                                                                              >  7:337616/1‑120 (MQ=255)
                                         aCCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCGATGGTGGTgggg                                                                                                                              <  8:337616/120‑1 (MQ=255)
                                         aCCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCGATGCTGGTgggg                                                                                                                              <  6:150230/120‑1 (MQ=255)
                                         aCCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCGATGCTGGTgggg                                                                                                                              >  5:150230/1‑120 (MQ=255)
                                         aCCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTTGTGGGGGAATCCATCGAGCCg                                                                                                               <  5:226134/135‑1 (MQ=255)
                                         aCCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCg                                                                                                               >  6:226134/1‑135 (MQ=255)
                                          ccccTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAAtc                                                                                                >  1:269189/1‑149 (MQ=255)
                                          ccccTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAAt                                                                                                 <  3:302957/148‑1 (MQ=255)
                                            ccTTCGGCCCAGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAAtctc                                                                                              <  2:328010/149‑1 (MQ=255)
                                             cTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTcc                                                                                             >  1:373738/1‑149 (MQ=255)
                                                cGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTg                                                                                           <  3:170224/148‑1 (MQ=255)
                                                  gCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATg                                                                                        >  2:325934/1‑149 (MQ=255)
                                                    cccGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGcc                                                                                      >  8:110248/1‑149 (MQ=255)
                                                    cccGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGcc                                                                                      >  6:237801/1‑149 (MQ=255)
                                                     ccGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGccc                                                                                     >  2:207614/1‑149 (MQ=255)
                                                            cAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGttt                                                                              <  1:207614/149‑1 (MQ=255)
                                                                 gTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCaa                                                                         >  7:176655/1‑149 (MQ=255)
                                                                    cATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGACCGGCGTTTATCAActc                                                                      >  1:218220/1‑149 (MQ=255)
                                                                       ttACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGc                                                                   >  5:364878/1‑149 (MQ=255)
                                                                           cGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTgggg                                                                                                                              >  8:331223/1‑86 (MQ=255)
                                                                           cGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTgggg                                                                                                                              <  7:331223/86‑1 (MQ=255)
                                                                             ttAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCg                                                                                    >  5:108448/1‑126 (MQ=255)
                                                                             ttAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCg                                                                                    <  6:108448/126‑1 (MQ=255)
                                                                                   aTAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTaa                                                        >  5:113653/1‑148 (MQ=255)
                                                                                    tAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTg                                                                                           <  5:15796/112‑1 (MQ=255)
                                                                                    tAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTg                                                                                           >  6:15796/1‑112 (MQ=255)
                                                                                    tAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTaaaa                                                      <  7:63041/149‑1 (MQ=255)
                                                                                      aGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACt                                                    <  7:179457/149‑1 (MQ=255)
                                                                                        tttttCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTgg                                                  <  1:223221/149‑1 (MQ=255)
                                                                                        tttttCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTgg                                                  >  2:223221/1‑149 (MQ=255)
                                                                                              aTTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGt                                                 <  7:350850/144‑1 (MQ=255)
                                                                                              aTTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGt                                                 >  8:350850/1‑144 (MQ=255)
                                                                                                  aGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATg                                        <  5:99858/149‑1 (MQ=255)
                                                                                                   gcgcAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGa                                       <  8:176655/149‑1 (MQ=255)
                                                                                                       aaTACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCg                                    >  5:215839/1‑148 (MQ=255)
                                                                                                       aaTACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGc                                   >  1:338330/1‑149 (MQ=255)
                                                                                                         tACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCAc                                 <  6:113653/149‑1 (MQ=255)
                                                                                                         tACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCAc                                 <  1:325934/149‑1 (MQ=255)
                                                                                                                   cgcCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGAGTTTATCAACTCTGCCACGCCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTAt                       >  1:281762/1‑149 (MQ=255)
                                                                                                                         cgcgTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGgctgc                 <  7:110248/149‑1 (MQ=255)
                                                                                                                           cgTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGACTCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGGctgctg               <  2:289901/149‑1 (MQ=255)
                                                                                                                                 tataCCCTGCAAGTCGTTCAATGCTTGTGGGGTAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGGCTGCTGttttt          <  2:281762/148‑1 (MQ=255)
                                                                                                                                 tataCCCGGCAAGACGTTCAATGCTGGTGTGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCGGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGGCTGCTGtttt           <  2:172758/147‑1 (MQ=255)
                                                                                                                                 tataCCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGGCTGCTGttttt          <  6:124757/148‑1 (MQ=255)
                                                                                                                                 tataCCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGGCTGCTGttttt          <  6:215839/148‑1 (MQ=255)
                                                                                                                                      ccGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGGCTGCTGTTTTTGCGCAg    <  7:318473/149‑1 (MQ=255)
                                                                                                                                      ccGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGGCTGCTGTTTTTGCGCAg    >  8:318473/1‑149 (MQ=255)
                                                                                                                                        ggCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGGCTGCTGTTTTTGCGCAggc  <  2:338330/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
CTAACGCGCTAAATTTAGCCCCACGTTCGAGCAACATTTCTACCCCTTCGGCCCCGCTGGCAGGCGTCCATTTACCGTTAGCGATAAGTTTTTCATTGAGCGCAATACCATGCTGCGCCAGCGCGTCTTTATACCCGGCAAGACGTTCAATGCTGGTGGGGGAATCCATCGAGCCGGTAAGGAAAGCAATCTCCTGATGCCCGGCGTTTATCAACTCTGCCACGGCGTTAAAACTGGTCTGTTTATGATCGCACCAGACGCTATGGCTGCTGTTTTTGCGCAGGC  >  NC_000913/2838532‑2838816

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