Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,362 Δ1 bp 15.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.15.7% 297.6 / 55.4 108intergenic (+156/+16)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base A (43/47);  new base . (3/14);  total (47/61)
Fisher's exact test for biased strand distribution p-value = 3.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

Only 100 of 111 total aligned reads displayed.

CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAA  >  NC_000913/1446263‑1446460
                                                                                                    |                                                                                                  
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)
       atcttGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACg                                                                                             <  1:1369684/96‑1 (MQ=38)
        aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGc                                                                                            <  2:677341/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  2:3599464/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  2:3590474/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  1:1015870/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                           <  2:1663515/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:3239482/100‑1 (MQ=255)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCa                                                                                          <  1:174606/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)
               tatTCAAAAACTCCTGTCAAATTGCCTTTTGTCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATCGGCGTTGCATGTACTGCAACGCCAAttag                                                                                     <  1:3424727/100‑1 (MQ=37)
               tatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAttag                                                                                     <  2:3356380/100‑1 (MQ=37)
                atTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAttagt                                                                                    <  1:3438192/100‑1 (MQ=37)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                                   >  2:1064313/1‑100 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                                   >  1:824978/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)
                      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)
                                cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATTCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                    >  1:2512889/1‑100 (MQ=255)
                                cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                    <  1:3510473/100‑1 (MQ=255)
                                 aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCg                                                                   <  1:1973327/100‑1 (MQ=255)
                                        cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                            <  1:2980687/100‑1 (MQ=255)
                                        cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                            >  1:2475308/1‑100 (MQ=255)
                                         ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGc                                                           <  2:807744/100‑1 (MQ=255)
                                                 tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                                   <  1:3169979/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)
                                                        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)
                                                                  aaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                                  >  2:14198/1‑100 (MQ=255)
                                                                    tatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                                >  1:3310404/1‑100 (MQ=255)
                                                                    tatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                                <  1:2648255/100‑1 (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)
                                                                           gctAACTAATTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGt                         >  1:647368/3‑100 (MQ=25)
                                                                            ttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                        <  1:348664/100‑1 (MQ=255)
                                                                            ttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCATTGCTTAAATACCGtt                        >  2:1299792/1‑100 (MQ=255)
                                                                             tAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGAGCCATCCCAAAGCGTTGCTTAAATACCGTTg                       >  2:2351434/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:1605941/1‑100 (MQ=255)
                                                                                    ttGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg                <  1:1882292/100‑1 (MQ=255)
                                                                                    ttGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg                <  2:1520361/100‑1 (MQ=255)
                                                                                     tGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa               >  1:609161/1‑100 (MQ=38)
                                                                                     tGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa               >  2:3372533/1‑100 (MQ=255)
                                                                                      ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGACTTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACAGTTGAAAAATGa               >  1:1426553/1‑99 (MQ=37)
                                                                                      ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc              >  2:2685382/1‑100 (MQ=255)
                                                                                      ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc              >  1:1123029/1‑100 (MQ=255)
                                                                                       gCGTTGCAGTACATGCAACGCCTATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAct             >  2:1552850/1‑100 (MQ=255)
                                                                                       gCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc              <  1:1185357/99‑1 (MQ=255)
                                                                                        cGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc            >  1:3534961/1‑100 (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)
                                                                                                gTACATGCAGCGCCAATCAGTTAGCGGAATTTACGTTGATACTCGCCTGGCGTCATCCCAAAGCGTTGTTTAAATACCGTTGAAAAATGACTCTGGTCa      <  2:3014915/99‑1 (MQ=255)
                                                                                                atACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAg     >  1:2737330/2‑100 (MQ=255)
                                                                                                   cATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaa  >  2:2157186/1‑100 (MQ=255)
                                                                                                   cATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaa  >  1:2995255/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAA  >  NC_000913/1446263‑1446460

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