Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,362 Δ1 bp 23.6% intergenic (+156/+16) ydbL → / ← feaR DUF1318 family protein/transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,3620A.23.6% 279.5 / 90.5 110intergenic (+156/+16)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base A (45/39);  new base . (11/15);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

Only 100 of 112 total aligned reads displayed.

CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAA  >  NC_000913/1446263‑1446460
                                                                                                    |                                                                                                  
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACa                                                                                                    >  2:652604/1‑100 (MQ=255)
 tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCATTTGCCCTGAAAAATCCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACa                                                                                                    >  1:913664/1‑99 (MQ=255)
  gAAACAACGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTTAAAAATGCATAGTAAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATg                                                                                                  <  1:1461613/100‑1 (MQ=255)
    aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCa                                                                                                >  1:2623984/1‑100 (MQ=255)
    aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATATGTTTTAACTAATTGGCGTTGCAGTACATGCa                                                                                                >  2:22470/1‑100 (MQ=255)
     aCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCCGTTTTAACTAATTGGCGTTGCAGTACATGCaa                                                                                               <  2:1836987/100‑1 (MQ=255)
       aaGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACg                                                                                             <  1:271216/100‑1 (MQ=38)
       aaGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACg                                                                                             <  2:2707648/100‑1 (MQ=255)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  1:2063970/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  1:1674379/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  2:288602/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  2:1818160/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  2:1500246/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  2:1241657/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  1:699050/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGcc                                                                                           >  1:1719470/1‑100 (MQ=255)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCa                                                                                          <  2:1129508/100‑1 (MQ=38)
          ccgatctttCCAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAACATGCATAGTCAAATATCTGTTTTAGCTAATTGGCGTTGCATGTACTGCAACGCCa                                                                                          <  1:538278/93‑1 (MQ=25)
           tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCaa                                                                                         <  2:2389772/100‑1 (MQ=38)
           tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                         <  1:1924958/100‑1 (MQ=255)
           tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                         <  1:844651/100‑1 (MQ=255)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAt                                                                                        <  2:2220686/100‑1 (MQ=38)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAt                                                                                        <  2:2288746/100‑1 (MQ=38)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAt                                                                                        <  2:327966/100‑1 (MQ=38)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAt                                                                                        <  1:2573092/100‑1 (MQ=255)
             aaTATTCAAAAACTGCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAtt                                                                                       >  1:2239186/1‑100 (MQ=255)
             aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAtt                                                                                       <  2:738304/100‑1 (MQ=255)
             aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAAGGCATAGTAAAATATCTGGTTTAACTAATTGGCGTTGCAGTACATGCAACGGCAAtt                                                                                       <  1:1166691/100‑1 (MQ=255)
               tatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttag                                                                                     <  1:1465510/100‑1 (MQ=255)
               tatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttag                                                                                     <  2:317710/100‑1 (MQ=255)
                atTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAttagt                                                                                    <  1:1097158/100‑1 (MQ=37)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                                   <  1:1204066/100‑1 (MQ=255)
                   cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagttag                                                                                 <  2:1393288/100‑1 (MQ=255)
                   cAAAAACTCCTGTAAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagttag                                                                                 >  2:1979309/1‑100 (MQ=255)
                     aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCg                                                                               <  2:337281/100‑1 (MQ=255)
                          tCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCGGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAt                                                                           <  1:277742/99‑1 (MQ=255)
                           ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAttt                                                                         >  2:1845945/1‑100 (MQ=255)
                            cTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTa                                                                        >  2:923751/1‑100 (MQ=255)
                              gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                      >  2:2891785/1‑100 (MQ=255)
                                 aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCg                                                                   <  1:791655/100‑1 (MQ=255)
                                 aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCg                                                                   >  2:378344/1‑100 (MQ=255)
                                    ttGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                                <  2:1053524/100‑1 (MQ=255)
                                     tGCCTTTTGCCCAGAAAAATGCATAGTCAAGTATCTGGTTTAACTAATTGGCGTTACAGTACATGCAACGCAAATGAGTTAGCGTAATTTACATAGATa                                                                >  1:2811678/1‑99 (MQ=255)
                                      gCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACt                                                              >  2:2496428/1‑100 (MQ=255)
                                        cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                            >  2:2906098/1‑100 (MQ=255)
                                             gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTgg                                                       >  2:2550698/1‑100 (MQ=255)
                                             gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTgg                                                       >  2:1702466/1‑100 (MQ=255)
                                             gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTgg                                                       <  2:167830/100‑1 (MQ=255)
                                              cccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGc                                                      >  2:338817/1‑100 (MQ=255)
                                              cccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAATACATGCAACGCCAAGTAGGTAGCGGAGTTTACGTCGATACGCGCCTGGc                                                      >  1:2442367/1‑100 (MQ=255)
                                               ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                     >  1:2294814/1‑100 (MQ=255)
                                               ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                     >  1:1190817/1‑100 (MQ=255)
                                               ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAATTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                     <  2:2764080/100‑1 (MQ=255)
                                                 tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                     <  1:2118331/98‑1 (MQ=255)
                                                  gAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                                   >  2:821274/1‑99 (MQ=255)
                                                    aaaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATc                                                >  1:2606128/1‑100 (MQ=255)
                                                       aggCATAGTCAAATATCTGTGTAAACGATTTGGCGTTGCAGTACAGGCAACGCCAAGTAGTTAGCGGAATTTACGTCGATACTCGCCTTGACTCATcc                                               <  1:2146452/96‑1 (MQ=255)
                                                           aTAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                         >  2:704099/1‑100 (MQ=255)
                                                            tAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGACATCACAAAGCg                                        <  1:332931/100‑1 (MQ=255)
                                                              gTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                      >  2:1492250/1‑100 (MQ=255)
                                                              gTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                      <  2:1393877/100‑1 (MQ=255)
                                                              gTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                      >  2:2000800/1‑100 (MQ=255)
                                                               tCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                     <  2:1489906/100‑1 (MQ=255)
                                                                   atatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                    <  1:983275/97‑1 (MQ=255)
                                                                    tatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                                >  2:2107298/1‑100 (MQ=255)
                                                                     atCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaaa                               >  1:2354578/1‑100 (MQ=255)
                                                                      tCTGTTTTAACTAATTCGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                              <  1:1702459/100‑1 (MQ=255)
                                                                        tGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAc                            >  1:1167766/1‑100 (MQ=255)
                                                                            ttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                        <  2:437083/100‑1 (MQ=255)
                                                                            ttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCACAAAGCGTTGCTTAAATACCGtt                        >  1:2138654/1‑100 (MQ=255)
                                                                             tAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                       <  2:347916/100‑1 (MQ=255)
                                                                             tAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                       <  2:2624000/100‑1 (MQ=255)
                                                                             tAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                       <  2:595313/100‑1 (MQ=255)
                                                                               aCTAATTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                     >  1:632699/1‑100 (MQ=37)
                                                                                cTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                    <  2:1022855/100‑1 (MQ=255)
                                                                                  aaTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                  >  2:300746/1‑100 (MQ=38)
                                                                                  aaTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                  >  1:1488491/1‑100 (MQ=255)
                                                                                  aaTTGGCGTTGAATGTACTGCGAGGCCAATTAGTTCACGGAATTTACAACGATACCAACCTGGAGTCATAACAAACCGCTGCATAAATACCGTTGaaa                    >  1:1005298/1‑98 (MQ=255)
                                                                                   aTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                 >  1:1249351/1‑100 (MQ=38)
                                                                                    ttGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg                <  1:2529849/100‑1 (MQ=255)
                                                                                      ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc              >  2:2435473/1‑100 (MQ=38)
                                                                                      ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc              >  1:137396/1‑100 (MQ=38)
                                                                                      ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc              >  1:2077918/1‑100 (MQ=255)
                                                                                       gCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAct             >  2:10323/1‑100 (MQ=255)
                                                                                        cGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc            >  2:2481381/1‑100 (MQ=38)
                                                                                        cGTTGCAGTACATGCAACGCCAAATAGTTAACGGAATTTACGTCGATACACGCCTGGCGTAATCACAAAGCGTGGCTTAAATACAGTTGAAAAATGAc              >  1:2768112/1‑98 (MQ=255)
                                                                                         cgtttctatACATTCAACGCTAATTAGTTAGCGGAATTTACGTCGATACTCACCTGGCGTCATCCCAAATCGTTGCTTAAATACCGTTGAAAAATGActc            >  1:2851606/9‑100 (MQ=255)
                                                                                           tGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg         <  2:1321266/100‑1 (MQ=255)
                                                                                           tGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg         <  2:1354062/100‑1 (MQ=255)
                                                                                            gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg         >  1:1438382/1‑99 (MQ=255)
                                                                                            gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt        >  1:340192/1‑100 (MQ=255)
                                                                                            gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt        >  2:177815/1‑100 (MQ=255)
                                                                                              aGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa      <  2:462093/100‑1 (MQ=255)
                                                                                              aGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa      >  1:1718893/1‑100 (MQ=255)
                                                                                              aGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAACGTTGAAAAATGACTCTGGTCa      <  1:776648/100‑1 (MQ=255)
                                                                                                 tACATGCAACGCGAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGa    >  2:1921232/1‑100 (MQ=255)
                                                                                                 tACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGa    <  2:2310332/100‑1 (MQ=255)
                                                                                                 tACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGa    <  2:840410/100‑1 (MQ=255)
                                                                                                  aCATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaa   >  2:2541837/1‑100 (MQ=255)
                                                                                                   cATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaa  >  2:582069/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAA  >  NC_000913/1446263‑1446460

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