Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,362 Δ1 bp 17.7% 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.17.7% 438.2 / 96.2 158intergenic (+156/+16)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base A (73/56);  new base . (13/15);  total (86/72)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01

Only 100 of 164 total aligned reads displayed.

CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAA  >  NC_000913/1446263‑1446461
                                                                                                    |                                                                                                   
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACa                                                                                                     >  1:2058095/1‑100 (MQ=255)
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACa                                                                                                     >  1:2929403/1‑100 (MQ=255)
  gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATg                                                                                                   >  1:1665298/1‑100 (MQ=255)
   aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGc                                                                                                  <  1:1763876/100‑1 (MQ=255)
    aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCa                                                                                                 <  1:1310166/100‑1 (MQ=255)
        aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGc                                                                                             <  1:3455709/100‑1 (MQ=38)
         ggcgAATATTCAAAAACTAATGTCAAATTGCCTTTTGCCGTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGc                                                                                             >  1:1658883/4‑99 (MQ=255)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                            <  2:3796301/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                            <  2:128647/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCATGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGACGTTGCAGTACATGCAACGc                                                                                             >  1:3695371/1‑99 (MQ=255)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCa                                                                                           <  1:3467998/100‑1 (MQ=38)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATATGTTTTAACTAATTGGCGTTGCAGAACATGCAACGCCa                                                                                           <  1:2714010/100‑1 (MQ=255)
           tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                          >  1:2542700/1‑100 (MQ=255)
           tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                          >  2:2397110/1‑100 (MQ=255)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAt                                                                                         <  2:3008718/100‑1 (MQ=38)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAt                                                                                         <  2:2473039/100‑1 (MQ=38)
             actatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAACTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAtt                                                                                        >  2:2977671/3‑100 (MQ=255)
             aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAtt                                                                                        <  2:152134/100‑1 (MQ=38)
              atatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAtta                                                                                       <  2:3899633/100‑1 (MQ=255)
               tatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttag                                                                                      <  2:2786556/100‑1 (MQ=255)
                atTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagt                                                                                     <  2:2440815/100‑1 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                                    >  2:1085403/1‑100 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                                    >  2:385236/1‑100 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTATTAACATAAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                                    <  1:4153885/100‑1 (MQ=38)
                 ttCAAAAACACCTGTAAAATTGCCTTTTGCCCAGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCATTGCAGTACATGCAACGCGAAttaat                                                                                     >  1:597439/1‑97 (MQ=255)
                   cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCTAAttagttag                                                                                  >  1:3502583/1‑100 (MQ=255)
                   cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagttag                                                                                  <  2:2060032/100‑1 (MQ=255)
                   cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagttag                                                                                  >  2:2603318/1‑100 (MQ=255)
                   cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagttag                                                                                  <  1:3926045/100‑1 (MQ=255)
                    aaaaaCTCCTGTCAAATTGCCTTTTTCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAttagttaaa                                                                                 <  1:1137551/100‑3 (MQ=25)
                     aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCg                                                                                <  2:940857/100‑1 (MQ=255)
                     aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCg                                                                                <  1:2686108/100‑1 (MQ=255)
                       aaCTCCTGTCAAATTGCCTTTTGCTCTGAAAAATGCAGAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGa                                                                              <  1:212345/100‑1 (MQ=255)
                          tCCTGTCAAATTGCATTTTAAACTAAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAtt                                                                           <  1:4089732/100‑1 (MQ=38)
                           ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAttt                                                                          <  1:2214397/100‑1 (MQ=255)
                             tGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTAc                                                                        <  2:4307465/100‑1 (MQ=255)
                                cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                     <  1:2922343/100‑1 (MQ=255)
                                  aaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCAAATTAGTTAGCGGAATTTACGTCGa                                                                   >  1:1093268/1‑100 (MQ=255)
                                  aaTTGCCTTTTGCACTGAAAAATTAATAGTCAAATATCTGTTTTAACTAATTGGCGTTGTAGTACATGCAACGCCTATTAGTTAGCGGAATTTACGt                                                                      >  1:134833/1‑97 (MQ=255)
                                   aTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGAt                                                                  >  2:591320/1‑100 (MQ=255)
                                       ccTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTc                                                              >  2:1976168/1‑100 (MQ=255)
                                          tttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGcc                                                           <  1:3813187/100‑1 (MQ=255)
                                            tGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCCACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCt                                                          >  2:3632042/1‑99 (MQ=255)
                                                cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                                     >  1:2196662/1‑100 (MQ=255)
                                                cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                                     >  1:1482703/1‑100 (MQ=255)
                                                 tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                                    >  1:4012496/1‑100 (MQ=255)
                                                      aaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTTATccc                                               >  1:1345040/1‑100 (MQ=255)
                                                        tGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaa                                             >  1:3018608/1‑100 (MQ=255)
                                                         gCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaaa                                            >  2:436807/1‑100 (MQ=255)
                                                         gCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaaa                                            <  2:3887673/100‑1 (MQ=255)
                                                          cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAg                                           >  2:70984/1‑100 (MQ=255)
                                                            tAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCg                                         <  1:2281979/100‑1 (MQ=255)
                                                             aGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGt                                        <  1:2865497/100‑1 (MQ=255)
                                                             aGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGt                                        <  2:3221984/100‑1 (MQ=255)
                                                               tCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                      <  2:2443690/100‑1 (MQ=255)
                                                                cAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                     <  2:3942811/100‑1 (MQ=255)
                                                                 aagtatCTGTTGTAACGAATTGGCGTTGCAGTACGTGCCACGCCAATTAGTTAGCGGAATTTACGTCGATACGCGCCTGGCGTCATCCCAAAGCGTTGCt                                    <  1:3436312/97‑1 (MQ=255)
                                                                  aaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                                   <  2:1004322/100‑1 (MQ=255)
                                                                   atatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTa                                  <  2:3768279/100‑1 (MQ=255)
                                                                    tatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                                 >  2:2689890/1‑100 (MQ=255)
                                                                      tCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                               >  2:1080104/1‑100 (MQ=255)
                                                                         gTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                            <  2:1133709/100‑1 (MQ=255)
                                                                         gTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                            >  2:2336145/1‑100 (MQ=255)
                                                                         gTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAACTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                            >  1:910480/1‑100 (MQ=255)
                                                                          ttttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCg                           >  2:1427557/1‑100 (MQ=255)
                                                                            ttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                         <  2:1775870/100‑1 (MQ=255)
                                                                             tAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATCCCGTTg                        >  1:2962649/1‑100 (MQ=255)
                                                                             tAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                        <  1:2348478/100‑1 (MQ=255)
                                                                              aaCTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                       <  2:491020/100‑1 (MQ=255)
                                                                               aCTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                      >  2:1918978/1‑100 (MQ=255)
                                                                                cTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                     >  2:482633/1‑100 (MQ=255)
                                                                                  aaTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                   >  2:394110/1‑100 (MQ=255)
                                                                                   aTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  2:2816898/1‑100 (MQ=38)
                                                                                   aTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  2:3718804/1‑100 (MQ=38)
                                                                                   aTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  1:3321457/1‑100 (MQ=38)
                                                                                   aTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  1:187246/1‑100 (MQ=38)
                                                                                    ttGGCGTTGCAGTACATGCAACGTCAATTAGTTAGCTGAATTTACGTCGATAATCGCCTGGCCTCATCCCAAAGCGTTGCTTAAATACAGTTGAAAAAt                  <  1:1816772/99‑1 (MQ=255)
                                                                                     tGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                >  2:287110/1‑100 (MQ=38)
                                                                                     tGGCGTTGCATCTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                >  1:2301226/1‑100 (MQ=38)
                                                                                      ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:1162366/1‑100 (MQ=38)
                                                                                      ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:3634811/1‑100 (MQ=38)
                                                                                      ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:757142/1‑100 (MQ=38)
                                                                                      ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:3053527/1‑100 (MQ=255)
                                                                                      ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  1:939505/1‑100 (MQ=255)
                                                                                        cGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc             <  2:2514792/100‑1 (MQ=255)
                                                                                        cGTTGCAGTACATGCAACGCAAATAAGTTAGCGGAATTTCCGTCGATACTCGCCTGGCGTCATCCCAAAGCATTCCTTAAATACCGTTGAAAAATGActc             >  1:1192078/1‑100 (MQ=255)
                                                                                         gTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActct            <  2:1831912/100‑1 (MQ=255)
                                                                                          ttttaGTACATGCAACGCCAATTAGTTATCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTg           <  1:2040547/96‑1 (MQ=255)
                                                                                           tGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg          <  2:4232055/100‑1 (MQ=255)
                                                                                            gCAGTACATGCAACGCCAATTAGTTCGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCGAAAGCGTTGCTTTAATACGGTGGAAAAATGCCTCTGGt         >  1:3279877/1‑100 (MQ=255)
                                                                                            gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt         >  1:645736/1‑100 (MQ=255)
                                                                                             cAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc        >  2:1715368/1‑100 (MQ=255)
                                                                                                gTACATGCAAGGGCATTTAGTTAGCGGAATTTACGTCGAGACTCGCCTGGCTTTATCCCAAGGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc        <  1:3334544/98‑1 (MQ=255)
                                                                                                gTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAATCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt         >  1:3669681/1‑97 (MQ=255)
                                                                                                gTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAg      >  2:4287102/1‑100 (MQ=255)
                                                                                                gTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAg      >  2:3268654/1‑100 (MQ=255)
                                                                                                  aCATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaa    >  1:1112155/1‑100 (MQ=255)
                                                                                                   cATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaa   <  2:2363703/100‑1 (MQ=255)
                                                                                                    aTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaan  <  1:482564/100‑2 (MQ=255)
                                                                                                    aTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaaa  >  2:1106991/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAA  >  NC_000913/1446263‑1446461

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: