Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,357:1 +T 15.5% 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.T15.5% 358.1 / 68.2 110intergenic (+151/+21)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base . (47/46);  new base T (3/14);  total (50/60)
Fisher's exact test for biased strand distribution p-value = 1.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-02

Only 100 of 122 total aligned reads displayed.

TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCA  >  NC_000913/1446258‑1446456
                                                                                                    |                                                                                                   
ttttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTAATAAAAAGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCa                                                                                                      <  2:2949300/100‑1 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                    <  1:407091/100‑1 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                    >  1:3567295/1‑100 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                    <  1:2925967/100‑1 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                    >  2:505418/1‑100 (MQ=255)
    ccTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGGAc                                                                                                 >  2:3108357/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGTTGTTGCAGTa                                                                                                  >  1:3324293/1‑98 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACa                                                                                                >  2:2095808/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGACAAATATCTGTTTTAACTAATAGGCGATGCAGTACa                                                                                                >  1:1569169/1‑100 (MQ=255)
      tGAAACAAGGTGCATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATACATACTCAACTCTCTGTTTTAACTAATTAGCGTTGCAGTACa                                                                                                >  2:3488023/1‑99 (MQ=255)
      tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACAt                                                                                               >  2:1476137/1‑100 (MQ=255)
       gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATg                                                                                              <  1:1408271/100‑1 (MQ=255)
        aaaCAAGGTGAATATTCAAAATCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGc                                                                                             <  1:2547387/100‑1 (MQ=255)
        aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACAt                                                                                               >  1:3574442/1‑98 (MQ=255)
        aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGc                                                                                             <  1:1121226/100‑1 (MQ=255)
             aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGc                                                                                        <  2:677341/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                       <  2:1663515/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                       <  1:1015870/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                       <  2:3599464/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGcc                                                                                       <  2:2600484/100‑1 (MQ=255)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTAAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                       <  1:1061496/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTGTGCCCTGAAAAATGCATAATCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACTcc                                                                                       <  2:2870803/100‑1 (MQ=255)
               gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCa                                                                                      <  1:2141528/100‑1 (MQ=38)
               gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCa                                                                                      <  2:1548571/100‑1 (MQ=38)
               gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCa                                                                                      <  2:1705968/100‑1 (MQ=255)
               gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCa                                                                                      <  1:174606/100‑1 (MQ=255)
               gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCa                                                                                      <  1:3239482/100‑1 (MQ=255)
                tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTATCTAATTGGCGTTGCAGTTCATGCAACGccta                                                                                     <  1:2641120/100‑3 (MQ=255)
                tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                     <  1:654450/100‑1 (MQ=255)
                tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                     >  1:470004/1‑100 (MQ=255)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAt                                                                                    <  1:628195/100‑1 (MQ=38)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAt                                                                                    <  1:172933/100‑1 (MQ=255)
                    tatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAttag                                                                                 <  2:3356380/100‑1 (MQ=37)
                     atTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAttagt                                                                                <  1:3438192/100‑1 (MQ=37)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                               >  2:1064313/1‑100 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                               >  2:2883170/1‑100 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAATCAAATCTCTGTTTTAACTAATTGAAGTTGCAGTACATACAACACCAAtaagt                                                                                >  1:1390315/1‑99 (MQ=255)
                       tCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtta                                                                              >  2:2942735/1‑100 (MQ=255)
                        cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagttag                                                                             >  2:2189274/1‑100 (MQ=255)
                           aaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCgg                                                                          <  1:2672036/100‑1 (MQ=255)
                           aaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCgg                                                                          <  1:2342295/100‑1 (MQ=255)
                                ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAttt                                                                     <  2:2553216/100‑1 (MQ=255)
                                 cTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTa                                                                    <  1:1911530/100‑1 (MQ=255)
                                  tGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTAc                                                                   >  2:1050091/1‑100 (MQ=255)
                                   gTCAAATTGCCTTTTGCCCTGAAAAATGTATAGTCAAATATCTGTTTTAACTAATTGGCGCTGCAGTATATGCAAAGCCAATTAGTTAGCGCAATTTa                                                                    >  2:1794042/1‑98 (MQ=255)
                                   gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                  >  1:2094821/1‑100 (MQ=255)
                                    tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTAGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                  >  1:1829356/1‑99 (MQ=255)
                                     cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                <  1:3510473/100‑1 (MQ=255)
                                      aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCg                                                               <  1:1973327/100‑1 (MQ=255)
                                             cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                        >  1:2475308/1‑100 (MQ=255)
                                             cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                        <  1:2980687/100‑1 (MQ=255)
                                             cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTAGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                        >  2:2826554/1‑100 (MQ=255)
                                              ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGc                                                       <  2:807744/100‑1 (MQ=255)
                                                       taaaaaTGCATAGTTAAATATCTGGTTTAACTAATTGTCGTTGCATTACATCCAACGCCAATTAGTTAGCGGAATTTACGTCGATAATCGCCt                                                     <  2:3226760/92‑1 (MQ=38)
                                                       gAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCa                                              <  1:2365832/100‑1 (MQ=255)
                                                       gAAAAATACATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATTCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTg                                                    <  2:2605814/94‑1 (MQ=255)
                                                           aaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATccc                                          <  2:1388980/100‑1 (MQ=255)
                                                             tGCATAGTCAAATATCTGTTTTAAGTAATTGGCGTTTCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaa                                        <  2:3828335/100‑1 (MQ=255)
                                                               cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGTCGTCATCCCaacg                                      <  1:2657373/100‑3 (MQ=255)
                                                               cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAg                                      >  1:390961/1‑100 (MQ=255)
                                                                    tCAAATATCTGTTTTAACTAATTGGCGTTGAAGTACATGCAACTCCAATTAGTTAGCGGAATTTACGTCTATACTCGCCTGGCGTCATCCCAAAGCGTTg                                 <  2:3840129/100‑1 (MQ=255)
                                                                     cAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                >  2:1483359/1‑100 (MQ=255)
                                                                     cAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                <  2:1920135/100‑1 (MQ=255)
                                                                      aaaTATCTGTTTTAACTAATTTGATTTGCATTACATTCAACTCTAATTAGTTAGCGATATCTACCTCGTTATTCGCCTTGCGTCATCCCAAAGCGTTGc                                <  2:2435904/99‑1 (MQ=255)
                                                                         tatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                            <  1:2648255/100‑1 (MQ=255)
                                                                         tatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                            >  1:3310404/1‑100 (MQ=255)
                                                                          atCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaaa                           <  2:62664/100‑1 (MQ=255)
                                                                          atCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCCTaaa                           <  2:1438722/100‑1 (MQ=255)
                                                                           tCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                          <  2:1590406/100‑1 (MQ=255)
                                                                           tCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                          >  2:1787346/1‑100 (MQ=255)
                                                                            cTGTTTTAACTAATTGGCGTTGCAGTACATGCCACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaaa                           <  2:1709627/98‑1 (MQ=255)
                                                                            cTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCTAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATa                         <  1:3206634/100‑1 (MQ=255)
                                                                            cTGTTTTAACTAATTGGCCTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTAGCGTCATCCCAAAGCGTTGCTTAAATa                         >  2:2923542/1‑100 (MQ=255)
                                                                             tGTTTTAACTAATTGGCGATGCAGTAAATGAAACGCCAAATAGTTAGCGGAATTTACGTCGATACTCGCCAGGCGTCATCCCAAAGCGTTACTTAAATAc                        >  2:3574125/1‑100 (MQ=255)
                                                                               ttttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCg                      >  2:165998/1‑100 (MQ=255)
                                                                                tttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGt                     <  1:3328400/100‑1 (MQ=255)
                                                                                gctAACTAATTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGt                     >  1:647368/3‑100 (MQ=25)
                                                                                 ttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                    <  1:348664/100‑1 (MQ=255)
                                                                                 ttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCATTGCTTAAATACCGtt                    >  2:1299792/1‑100 (MQ=255)
                                                                                   aaCTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGAGTTGCTTACATACCGTTGa                  >  2:2801816/1‑100 (MQ=255)
                                                                                   aaCTAATTGGCGTTGCAGTACATGCAACGCCAATAAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                  >  2:513812/1‑100 (MQ=255)
                                                                                    aCTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                 >  1:2298724/1‑100 (MQ=255)
                                                                                       aaTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa              >  2:1043437/1‑100 (MQ=255)
                                                                                        aTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt             <  1:3514410/100‑1 (MQ=255)
                                                                                         ttGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg            <  2:1520361/100‑1 (MQ=255)
                                                                                         ttGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg            >  2:1605941/1‑100 (MQ=255)
                                                                                         ttGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg            <  1:1882292/100‑1 (MQ=255)
                                                                                          tGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa           >  1:609161/1‑100 (MQ=38)
                                                                                          tGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGTAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa           <  1:166003/100‑1 (MQ=255)
                                                                                          tGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa           >  2:3372533/1‑100 (MQ=255)
                                                                                           ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGACTTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACAGTTGAAAAATGa           >  1:1426553/1‑99 (MQ=37)
                                                                                           ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc          >  1:1123029/1‑100 (MQ=255)
                                                                                            gCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc          <  1:1185357/99‑1 (MQ=255)
                                                                                               ttGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc        <  1:1857662/98‑1 (MQ=255)
                                                                                                 gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActctcat    >  1:3163908/1‑97 (MQ=255)
                                                                                                 gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt    >  2:2244040/1‑100 (MQ=255)
                                                                                                  cAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc   <  2:803388/100‑1 (MQ=255)
                                                                                                   aGTACGTGCAACGCGAATTAGTTAGCGGAATATACGTCGAGACTCGACTGACATCAGCCCAAAGCGTAGCATGAAGACCGTGGAAAAATGACTCTgatga  >  2:2372254/1‑96 (MQ=255)
                                                                                                   aGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAATGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa  <  2:2094802/100‑1 (MQ=255)
                                                                                                   aGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa  <  1:2942761/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCA  >  NC_000913/1446258‑1446456

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