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.

CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAA  >  NC_000913/1446263‑1446461
                                                                                                    |                                                                                                   
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACa                                                                                                     >  1:2058095/1‑100 (MQ=255)
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACa                                                                                                     >  1:2929403/1‑100 (MQ=255)
  gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATg                                                                                                   >  1:1665298/1‑100 (MQ=255)
   aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGc                                                                                                  <  1:1763876/100‑1 (MQ=255)
    aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCa                                                                                                 <  1:1310166/100‑1 (MQ=255)
        aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGc                                                                                             <  1:3455709/100‑1 (MQ=38)
         ggcgAATATTCAAAAACTAATGTCAAATTGCCTTTTGCCGTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGc                                                                                             >  1:1658883/4‑99 (MQ=255)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                            <  2:3796301/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                            <  2:128647/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCATGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGACGTTGCA‑GTACATGCAACGc                                                                                             >  1:3695371/1‑99 (MQ=255)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCa                                                                                           <  1:3467998/100‑1 (MQ=38)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATATGTTTTAACTAATTGGCGTTGCA‑GAACATGCAACGCCa                                                                                           <  1:2714010/100‑1 (MQ=255)
           tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCaa                                                                                          >  1:2542700/1‑100 (MQ=255)
           tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCaa                                                                                          >  2:2397110/1‑100 (MQ=255)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAt                                                                                         <  2:3008718/100‑1 (MQ=38)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAt                                                                                         <  2:2473039/100‑1 (MQ=38)
             actatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAACTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtt                                                                                        >  2:2977671/3‑100 (MQ=255)
             aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAtt                                                                                        <  2:152134/100‑1 (MQ=38)
              atatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtta                                                                                       <  2:3899633/100‑1 (MQ=255)
               tatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttag                                                                                      <  2:2786556/100‑1 (MQ=255)
                atTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagt                                                                                     <  2:2440815/100‑1 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                                    >  2:1085403/1‑100 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                                    >  2:385236/1‑100 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTATTAACATAAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                                    <  1:4153885/100‑1 (MQ=38)
                 ttCAAAAACACCTGTAAAATTGCCTTTTGCCCAGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCATTGCA‑GTACATGCAACGCGAAttaat                                                                                     >  1:597439/1‑97 (MQ=255)
                   cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCTAAttagttag                                                                                  >  1:3502583/1‑100 (MQ=255)
                   cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagttag                                                                                  <  2:2060032/100‑1 (MQ=255)
                   cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagttag                                                                                  >  2:2603318/1‑100 (MQ=255)
                   cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagttag                                                                                  <  1:3926045/100‑1 (MQ=255)
                    aaaaaCTCCTGTCAAATTGCCTTTTTCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAttagttaaa                                                                                 <  1:1137551/100‑3 (MQ=25)
                     aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCg                                                                                <  2:940857/100‑1 (MQ=255)
                     aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCg                                                                                <  1:2686108/100‑1 (MQ=255)
                       aaCTCCTGTCAAATTGCCTTTTGCTCTGAAAAATGCAGAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGa                                                                              <  1:212345/100‑1 (MQ=255)
                          tCCTGTCAAATTGCATTTTAAACTAAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAAtt                                                                           <  1:4089732/100‑1 (MQ=38)
                           ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAAttt                                                                          <  1:2214397/100‑1 (MQ=255)
                             tGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTAc                                                                        <  2:4307465/100‑1 (MQ=255)
                                cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                     <  1:2922343/100‑1 (MQ=255)
                                  aaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCAAATTAGTTAGCGGAATTTACGTCGa                                                                   >  1:1093268/1‑100 (MQ=255)
                                  aaTTGCCTTTTGCACTGAAAAATTAATAGTCAAATATCTGTTTTAACTAATTGGCGTTGTA‑GTACATGCAACGCCTATTAGTTAGCGGAATTTACGt                                                                      >  1:134833/1‑97 (MQ=255)
                                   aTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGAt                                                                  >  2:591320/1‑100 (MQ=255)
                                       ccTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTc                                                              >  2:1976168/1‑100 (MQ=255)
                                          tttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGcc                                                           <  1:3813187/100‑1 (MQ=255)
                                            tGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCCACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCt                                                          >  2:3632042/1‑99 (MQ=255)
                                                cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                                     >  1:2196662/1‑100 (MQ=255)
                                                cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                                     >  1:1482703/1‑100 (MQ=255)
                                                 tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                                    >  1:4012496/1‑100 (MQ=255)
                                                      aaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTTATccc                                               >  1:1345040/1‑100 (MQ=255)
                                                        tGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaa                                             >  1:3018608/1‑100 (MQ=255)
                                                         gCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaaa                                            >  2:436807/1‑100 (MQ=255)
                                                         gCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaaa                                            <  2:3887673/100‑1 (MQ=255)
                                                          cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAg                                           >  2:70984/1‑100 (MQ=255)
                                                            tAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCg                                         <  1:2281979/100‑1 (MQ=255)
                                                             aGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGt                                        <  1:2865497/100‑1 (MQ=255)
                                                             aGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGt                                        <  2:3221984/100‑1 (MQ=255)
                                                               tCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                      <  2:2443690/100‑1 (MQ=255)
                                                                cAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                     <  2:3942811/100‑1 (MQ=255)
                                                                 aagtatCTGTTGTAACGAATTGGCGTTGCA‑GTACGTGCCACGCCAATTAGTTAGCGGAATTTACGTCGATACGCGCCTGGCGTCATCCCAAAGCGTTGCt                                    <  1:3436312/97‑1 (MQ=255)
                                                                  aaTATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                                   <  2:1004322/100‑1 (MQ=255)
                                                                   atatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTa                                  <  2:3768279/100‑1 (MQ=255)
                                                                    tatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                                 >  2:2689890/1‑100 (MQ=255)
                                                                      tCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                               >  2:1080104/1‑100 (MQ=255)
                                                                         gTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                            <  2:1133709/100‑1 (MQ=255)
                                                                         gTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                            >  2:2336145/1‑100 (MQ=255)
                                                                         gTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAACTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                            >  1:910480/1‑100 (MQ=255)
                                                                          ttttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCg                           >  2:1427557/1‑100 (MQ=255)
                                                                            ttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                         <  2:1775870/100‑1 (MQ=255)
                                                                             tAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATCCCGTTg                        >  1:2962649/1‑100 (MQ=255)
                                                                             tAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                        <  1:2348478/100‑1 (MQ=255)
                                                                              aaCTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                       <  2:491020/100‑1 (MQ=255)
                                                                               aCTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                      >  2:1918978/1‑100 (MQ=255)
                                                                                cTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                     >  2:482633/1‑100 (MQ=255)
                                                                                  aaTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                   >  2:394110/1‑100 (MQ=255)
                                                                                   aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  2:2816898/1‑100 (MQ=38)
                                                                                   aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  2:3718804/1‑100 (MQ=38)
                                                                                   aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  1:3321457/1‑100 (MQ=38)
                                                                                   aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  1:187246/1‑100 (MQ=38)
                                                                                    ttGGCGTTGCA‑GTACATGCAACGTCAATTAGTTAGCTGAATTTACGTCGATAATCGCCTGGCCTCATCCCAAAGCGTTGCTTAAATACAGTTGAAAAAt                  <  1:1816772/99‑1 (MQ=255)
                                                                                     tGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                >  2:287110/1‑100 (MQ=38)
                                                                                     tGGCGTTGCATCTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                >  1:2301226/1‑100 (MQ=38)
                                                                                      ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:1162366/1‑100 (MQ=38)
                                                                                      ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:3634811/1‑100 (MQ=38)
                                                                                      ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:757142/1‑100 (MQ=38)
                                                                                      ggCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:3053527/1‑100 (MQ=255)
                                                                                      ggCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  1:939505/1‑100 (MQ=255)
                                                                                        cGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc             <  2:2514792/100‑1 (MQ=255)
                                                                                        cGTTGCA‑GTACATGCAACGCAAATAAGTTAGCGGAATTTCCGTCGATACTCGCCTGGCGTCATCCCAAAGCATTCCTTAAATACCGTTGAAAAATGActc             >  1:1192078/1‑100 (MQ=255)
                                                                                         gTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActct            <  2:1831912/100‑1 (MQ=255)
                                                                                          tttta‑GTACATGCAACGCCAATTAGTTATCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTg           <  1:2040547/96‑1 (MQ=255)
                                                                                           tGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg          <  2:4232055/100‑1 (MQ=255)
                                                                                            gCA‑GTACATGCAACGCCAATTAGTTCGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCGAAAGCGTTGCTTTAATACGGTGGAAAAATGCCTCTGGt         >  1:3279877/1‑100 (MQ=255)
                                                                                            gCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt         >  1:645736/1‑100 (MQ=255)
                                                                                             cA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc        >  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)
                                                                                                    |                                                                                                   
CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAA  >  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: ‑