Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,357:1 +T 17.8% intergenic (+151/+21) 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,3571.T17.8% 490.1 / 120.7 157intergenic (+151/+21)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base . (73/56);  new base T (13/15);  total (86/71)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

Only 100 of 160 total aligned reads displayed.

TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTC  >  NC_000913/1446258‑1446455
                                                                                                    |                                                                                                  
  ttCCTGAAACGAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑Gt                                                                                                  <  1:1003986/100‑1 (MQ=255)
  ttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑Gt                                                                                                  <  1:4401920/100‑1 (MQ=255)
    ccTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTAc                                                                                                >  2:3264281/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACa                                                                                               >  1:2058095/1‑100 (MQ=255)
       gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATg                                                                                             >  1:3447508/1‑100 (MQ=255)
       gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATg                                                                                             >  1:1665298/1‑100 (MQ=255)
        aaaGAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCACTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGc                                                                                            <  1:739/100‑1 (MQ=255)
        aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGc                                                                                            <  1:1763876/100‑1 (MQ=255)
         aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCa                                                                                           <  1:1310166/100‑1 (MQ=255)
         aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCa                                                                                           >  1:2953861/1‑100 (MQ=255)
          aCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCaa                                                                                          >  2:134638/1‑100 (MQ=255)
            aaGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACg                                                                                        <  1:1647756/100‑1 (MQ=255)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                      <  2:3796301/100‑1 (MQ=38)
               gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCa                                                                                     <  1:3467998/100‑1 (MQ=38)
               gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCa                                                                                     >  1:4401406/1‑100 (MQ=255)
                tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCaa                                                                                    >  2:2397110/1‑100 (MQ=255)
                tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCaa                                                                                    >  1:1770452/1‑100 (MQ=255)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAt                                                                                   <  2:2473039/100‑1 (MQ=38)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAt                                                                                   <  2:3008718/100‑1 (MQ=38)
                 gAATATTAAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGcc                                                                                      <  1:3289985/97‑1 (MQ=255)
                  aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAtt                                                                                  <  2:152134/100‑1 (MQ=38)
                   atatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtta                                                                                 <  2:3899633/100‑1 (MQ=255)
                    tatGCAAAAACTCCTTTTAAATTGTTTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAttag                                                                                <  1:4042245/100‑1 (MQ=21)
                     atTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagt                                                                               <  2:2440815/100‑1 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATTCATCGTCAAtt                                                                                  >  1:3608452/1‑96 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                              >  2:385236/1‑100 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                              >  2:1085403/1‑100 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                              <  2:1684162/100‑1 (MQ=255)
                      ttCAAAAACACCTGTAAAATTGCCTTTTGCCCAGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCATTGCA‑GTACATGCAACGCGAAttaat                                                                               >  1:597439/1‑97 (MQ=255)
                       tCAAAAACTCCTGTCAAATTGCGTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagt                                                                               >  1:2547293/1‑98 (MQ=255)
                        cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCTAAttagttag                                                                            >  1:3502583/1‑100 (MQ=255)
                        cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagttag                                                                            >  2:2603318/1‑100 (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)
                                     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‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTg                                                   <  2:1128441/100‑1 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                >  1:3942786/1‑100 (MQ=255)
                                                     cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                               >  1:1482703/1‑100 (MQ=255)
                                                     cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                               >  1:2274675/1‑100 (MQ=255)
                                                     cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                               >  1:2290582/1‑100 (MQ=255)
                                                         aaaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATc                                           <  2:2724068/100‑1 (MQ=255)
                                                          aaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATcc                                          <  1:2841033/100‑1 (MQ=255)
                                                           aaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTTATccc                                         >  1:1345040/1‑100 (MQ=255)
                                                              gCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaaa                                      >  2:436807/1‑100 (MQ=255)
                                                               cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAg                                     >  2:70984/1‑100 (MQ=255)
                                                                  aGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGt                                  <  2:3221984/100‑1 (MQ=255)
                                                                   gTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                 >  2:291800/1‑100 (MQ=255)
                                                                    tCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGTCGTAATCCCAAAGc                                    >  1:287140/1‑96 (MQ=255)
                                                                    tCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                <  2:2443690/100‑1 (MQ=255)
                                                                    tCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                >  2:210835/1‑100 (MQ=255)
                                                                     cAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                               <  2:3942811/100‑1 (MQ=255)
                                                                     cAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                               <  2:4269912/100‑1 (MQ=255)
                                                                      aagtatCTGTTGTAACGAATTGGCGTTGCA‑GTACGTGCCACGCCAATTAGTTAGCGGAATTTACGTCGATACGCGCCTGGCGTCATCCCAAAGCGTTGCt                              <  1:3436312/97‑1 (MQ=255)
                                                                       aaTATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                             <  2:1004322/100‑1 (MQ=255)
                                                                         tatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                           <  2:3927497/100‑1 (MQ=255)
                                                                         tatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                           >  2:2689890/1‑100 (MQ=255)
                                                                           tCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                         >  2:1080104/1‑100 (MQ=255)
                                                                             tGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGGTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATa                        >  1:2469247/1‑99 (MQ=255)
                                                                              gTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                      >  2:2336145/1‑100 (MQ=255)
                                                                              gTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                      >  1:3811610/1‑100 (MQ=255)
                                                                              gTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAACTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                      >  1:910480/1‑100 (MQ=255)
                                                                                 ttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                   <  2:1775870/100‑1 (MQ=255)
                                                                                 ttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                   <  2:3464031/100‑1 (MQ=255)
                                                                                  tAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATCCCGTTg                  >  1:2962649/1‑100 (MQ=255)
                                                                                  tAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                  <  1:2348478/100‑1 (MQ=255)
                                                                                     cTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa               >  2:1276613/1‑100 (MQ=255)
                                                                                     cTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa               >  2:482633/1‑100 (MQ=255)
                                                                                      tAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaa              <  2:1956481/100‑1 (MQ=255)
                                                                                       aaTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa             >  2:2505036/1‑100 (MQ=255)
                                                                                       aaTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa             >  2:2833732/1‑100 (MQ=255)
                                                                                        aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  1:187246/1‑100 (MQ=38)
                                                                                        aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  2:2816898/1‑100 (MQ=38)
                                                                                        aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  2:3655282/1‑100 (MQ=38)
                                                                                        aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  2:3718804/1‑100 (MQ=38)
                                                                                        aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  1:3321457/1‑100 (MQ=38)
                                                                                        aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTCAGCGGAATTTACGTCGATACTCGCCTGGCATCATACCAAAGCGTGGCTTAAATACCGTTCaaaa              >  1:3034891/1‑98 (MQ=25)
                                                                                        aTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            <  2:386995/100‑1 (MQ=255)
                                                                                        aTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  2:3025366/1‑100 (MQ=255)
                                                                                         ttGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg           <  2:1731864/100‑1 (MQ=255)
                                                                                          tGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAGATACCGTTGAAAAATGa          >  2:1893004/1‑100 (MQ=255)
                                                                                           ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  1:2024804/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         >  1:939505/1‑100 (MQ=255)
                                                                                             cGTTGCA‑GTACATGCAACGCAAATAAGTTAGCGGAATTTCCGTCGATACTCGCCTGGCGTCATCCCAAAGCATTCCTTAAATACCGTTGAAAAATGActc       >  1:1192078/1‑100 (MQ=255)
                                                                                              gTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActct      >  2:3745930/1‑100 (MQ=255)
                                                                                              gTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc       <  1:2166771/99‑1 (MQ=255)
                                                                                                tGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg    <  2:4232055/100‑1 (MQ=255)
                                                                                                tGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg    <  2:4311876/100‑1 (MQ=255)
                                                                                                 gCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt   <  2:528241/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTC  >  NC_000913/1446258‑1446455

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