Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,357:1 +T 17.6% 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.6% 309.1 / 75.2 102intergenic (+151/+21)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base . (43/41);  new base T (10/8);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 110 total aligned reads displayed.

TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTC  >  NC_000913/1446258‑1446455
                                                                                                    |                                                                                                  
ttttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCa                                                                                                     <  1:1622062/100‑1 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                   >  1:3600912/1‑100 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                   >  1:3466745/1‑100 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                   <  2:1010814/100‑1 (MQ=255)
   tCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTa                                                                                                 <  2:1509207/100‑1 (MQ=255)
    ccTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTAc                                                                                                >  2:114935/1‑100 (MQ=255)
    ccTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTAc                                                                                                <  2:3277429/100‑1 (MQ=255)
      tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACan                                                                                              <  1:2328079/100‑2 (MQ=255)
      tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACAt                                                                                              <  2:3175660/100‑1 (MQ=255)
      tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTATGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACAt                                                                                              >  1:2587718/1‑100 (MQ=255)
       gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCatgtactg                                                                                              <  2:605960/100‑8 (MQ=255)
        aaaCAAGGTGAATATTCAAAATCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACAt                                                                                              <  1:3105909/98‑1 (MQ=255)
        aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCatgtactgc                                                                                             <  2:1158466/100‑9 (MQ=255)
         aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCatgtactgca                                                                                            <  1:2748353/100‑10 (MQ=255)
          aCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCaa                                                                                          <  1:1035547/100‑1 (MQ=255)
           cAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAAc                                                                                         <  2:25771/100‑1 (MQ=39)
            aaGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACg                                                                                        <  1:59030/100‑1 (MQ=38)
             aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGc                                                                                       <  1:2345498/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                      <  1:2191964/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGcc                                                                                      <  2:3679079/100‑1 (MQ=255)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTACCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                      <  2:2659699/100‑1 (MQ=38)
               gTGAATATTCAAAAACTCCTGTCATATTGCATTTTGCCCTGATAAATTCATATTCATATATCTTTTTTAACTACTTGGCGTTGCAGTACATGCatc                                                                                         <  2:3090401/96‑3 (MQ=255)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAt                                                                                   <  2:1744432/100‑1 (MQ=38)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAt                                                                                   <  1:2640291/100‑1 (MQ=255)
                 gAATATAAAAAAATTCTTGTCAAATTGCTTTTTGACCTAAAAAATGCATCCTCAAATATCTCTTTAAACTAATTGGCGTTGCAGTACTTGCAACGCCa                                                                                     <  2:1088691/98‑1 (MQ=255)
                  aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAtt                                                                                  <  2:34420/100‑1 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTTCCCTGAAAAATTCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGGACTGCAACGCCAAttagt                                                                               <  2:3790533/99‑1 (MQ=25)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                              >  2:3362742/1‑100 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                              >  2:2565789/1‑100 (MQ=255)
                        cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagttag                                                                            >  1:3657348/1‑100 (MQ=255)
                         aaaagCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAttagttaa                                                                            <  2:1341723/99‑2 (MQ=25)
                         aaaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGc                                                                           >  1:151965/1‑100 (MQ=255)
                          aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCg                                                                          <  1:3198414/100‑1 (MQ=255)
                           aaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCCGTTTTAACTAATTGGCGTTGCAATACATGCAAAGCCCATTAATTAGCgg                                                                         >  2:2159595/1‑100 (MQ=255)
                                ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAttt                                                                    <  2:399221/100‑1 (MQ=255)
                                    tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGt                                                                >  1:1835554/1‑100 (MQ=255)
                                     cAAATTTCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                               <  2:2937453/100‑1 (MQ=255)
                                     cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                               >  1:260704/1‑100 (MQ=255)
                                       aaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGa                                                             >  1:2564904/1‑100 (MQ=255)
                                        aTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGAt                                                            <  2:1182944/100‑1 (MQ=255)
                                         ttGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                           >  2:742685/1‑100 (MQ=255)
                                              ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGc                                                      >  2:456132/1‑100 (MQ=255)
                                              ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGc                                                      <  1:3362768/100‑1 (MQ=255)
                                              ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGc                                                      <  1:1906781/100‑1 (MQ=255)
                                               tttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGcc                                                     >  1:2814346/1‑100 (MQ=255)
                                                  gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCt                                                    >  2:3380917/1‑98 (MQ=255)
                                                   cccTGAAAAATGCATAGTAAAATATCTGTTTTAACTAATTGGCGTAGCAGTACATGCAACGCCTATTAGTTAGCGGAATTTACGTCGAGACTCTcc                                                     >  2:1265809/1‑96 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                <  2:1954369/100‑1 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                >  2:1768541/1‑100 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                <  1:2307984/100‑1 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                <  1:1692292/100‑1 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                >  1:1909312/1‑100 (MQ=255)
                                                       gAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                              >  2:2418760/1‑99 (MQ=255)
                                                         aaaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTAGCAGTACATGCAACGCCAATCAGTTAGCGGAATTTACGTCGATACTCGCCTGGGGt                                               >  2:1880771/1‑96 (MQ=255)
                                                           aaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATccc                                         <  2:679305/100‑1 (MQ=255)
                                                               cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAg                                     >  1:2899238/1‑100 (MQ=255)
                                                                aTAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                    >  2:1485099/1‑100 (MQ=255)
                                                                aTAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                    <  1:2061105/100‑1 (MQ=255)
                                                                 tAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCg                                   <  1:59482/100‑1 (MQ=255)
                                                                   gTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                 <  1:3065087/100‑1 (MQ=255)
                                                                    tCAAATATCTGTTTTGACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                <  1:3764253/100‑1 (MQ=255)
                                                                    tCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                >  1:2555654/1‑100 (MQ=255)
                                                                     cAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                               <  2:1567191/100‑1 (MQ=255)
                                                                     cAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                               >  1:3703894/1‑100 (MQ=255)
                                                                     cAAATATCTGTTTTAACTAATTGGCGTTGCAGTACAGGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                               <  1:1270057/100‑1 (MQ=255)
                                                                       aaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                             <  2:2564887/100‑1 (MQ=255)
                                                                         tatCTGTTTTAACTAATTGGCGTTGCAGTACATAAAACGCCTATTAGTTAGCAGTATTTACATCGATACCCCCCAGGCGTCAACCCAAAGCGt                                  >  2:1687317/1‑93 (MQ=255)
                                                                           tCTGTTTTAACTAATTGGCGTTGCAGTACATGAAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                         <  1:2565806/100‑1 (MQ=255)
                                                                             tGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATa                        >  2:479484/1‑99 (MQ=255)
                                                                              gTTTTAACTGATGGGCGTTGCAGTAAATGCAACGCCCATGAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGGCATCCCAAAGCGTTGCTTAAATa                        <  2:3414192/98‑1 (MQ=255)
                                                                              gTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                      >  1:3441897/1‑100 (MQ=255)
                                                                                 ttAACTAATTGGCTTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                   <  1:3065446/100‑1 (MQ=255)
                                                                                 ttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                   >  2:3404551/1‑100 (MQ=255)
                                                                                  tAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                  >  1:2191767/1‑100 (MQ=255)
                                                                                   aaCTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                 >  2:2380562/1‑100 (MQ=255)
                                                                                      tAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaa              <  2:2530711/100‑1 (MQ=255)
                                                                                      tAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaa              >  1:3654457/1‑100 (MQ=255)
                                                                                      tAATTGGCGTTGCAGTACATGAAACCCCAATTAGTTCGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa               <  2:2517715/99‑1 (MQ=255)
                                                                                       aaTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa             >  2:3079679/1‑100 (MQ=255)
                                                                                       aaTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa             <  1:2879625/100‑1 (MQ=255)
                                                                                        aTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  2:69705/1‑100 (MQ=38)
                                                                                        aTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  1:1885326/1‑100 (MQ=38)
                                                                                        aTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  2:1958267/1‑100 (MQ=255)
                                                                                        aTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  2:2025366/1‑100 (MQ=255)
                                                                                        aTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGACGATACTCGCCTGGCGTCATCCCAAAGCGTTGATTACATACAGTTGAAAAAt            >  1:414771/1‑100 (MQ=255)
                                                                                         ttGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg           <  1:2173268/100‑1 (MQ=255)
                                                                                         ttGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg           <  2:3098103/100‑1 (MQ=255)
                                                                                          tGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa          >  1:2090322/1‑100 (MQ=38)
                                                                                           ggCGTTGCATGTACTGCAACGCCTATTAGTTAGCGGAATTTACTTCGATACTCGCCTTGCGTCATCCCAAAGCTTTGCTTAAATACCGATGAAAAATGAc         >  2:1669320/1‑100 (MQ=37)
                                                                                           ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  2:3733127/1‑100 (MQ=38)
                                                                                           ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  1:1347847/1‑100 (MQ=38)
                                                                                           ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  2:1760565/1‑100 (MQ=38)
                                                                                           ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  1:2437951/1‑100 (MQ=255)
                                                                                           ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  1:728621/1‑100 (MQ=255)
                                                                                           ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  1:2436925/1‑100 (MQ=255)
                                                                                            gCGTTGCAGTACATGCAACGCCAATTAGTTAGCGCAATTTACGTCGATACTCGCCTGGCGTCATCCCAAACCGTTGCTTAAATACCGTTGAAAAATGAc         <  2:1270824/99‑1 (MQ=255)
                                                                                                tGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg    <  2:2845547/100‑1 (MQ=255)
                                                                                                tGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg    >  1:3161749/1‑100 (MQ=255)
                                                                                                 gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt   <  1:3404577/100‑1 (MQ=255)
                                                                                                  cAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc  <  1:3351679/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTC  >  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: