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.

CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAA  >  NC_000913/1446263‑1446460
                                                                                                    |                                                                                                  
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACa                                                                                                    >  2:2095808/1‑100 (MQ=255)
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGACAAATATCTGTTTTAACTAATAGGCGATGCA‑GTACa                                                                                                    >  1:1569169/1‑100 (MQ=255)
 tGAAACAAGGTGCATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATACATACTCAACTCTCTGTTTTAACTAATTAGCGTTGCA‑GTACa                                                                                                    >  2:3488023/1‑99 (MQ=255)
 tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACAt                                                                                                   >  2:1476137/1‑100 (MQ=255)
  gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATg                                                                                                  <  1:1408271/100‑1 (MQ=255)
   aaaCAAGGTGAATATTCAAAATCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGc                                                                                                 <  1:2547387/100‑1 (MQ=255)
   aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACAt                                                                                                   >  1:3574442/1‑98 (MQ=255)
   aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGc                                                                                                 <  1:1121226/100‑1 (MQ=255)
       atcttGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACg                                                                                             <  1:1369684/96‑1 (MQ=38)
        aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGc                                                                                            <  2:677341/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                           <  2:3599464/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                           <  2:3590474/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                           <  1:1015870/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                           <  2:1663515/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGcc                                                                                           <  2:2600484/100‑1 (MQ=255)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTAAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                           <  1:1061496/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTGTGCCCTGAAAAATGCATAATCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACTcc                                                                                           <  2:2870803/100‑1 (MQ=255)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCa                                                                                          <  1:2141528/100‑1 (MQ=38)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCa                                                                                          <  2:1548571/100‑1 (MQ=38)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCa                                                                                          <  2:1705968/100‑1 (MQ=255)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCa                                                                                          <  1:3239482/100‑1 (MQ=255)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCa                                                                                          <  1:174606/100‑1 (MQ=255)
           tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTATCTAATTGGCGTTGCA‑GTTCATGCAACGccta                                                                                         <  1:2641120/100‑3 (MQ=255)
           tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCaa                                                                                         <  1:654450/100‑1 (MQ=255)
           tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCaa                                                                                         >  1:470004/1‑100 (MQ=255)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAt                                                                                        <  1:628195/100‑1 (MQ=38)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAt                                                                                        <  1:172933/100‑1 (MQ=255)
               tatTCAAAAACTCCTGTCAAATTGCCTTTTGTCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATCGGCGTTGCATGTAC‑TGCAACGCCAAttag                                                                                     <  1:3424727/100‑1 (MQ=37)
               tatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAttag                                                                                     <  2:3356380/100‑1 (MQ=37)
                atTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAttagt                                                                                    <  1:3438192/100‑1 (MQ=37)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                                   >  2:1064313/1‑100 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                                   >  1:824978/1‑100 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                                   >  2:2883170/1‑100 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAATCAAATCTCTGTTTTAACTAATTGAAGTTGCA‑GTACATACAACACCAAtaagt                                                                                    >  1:1390315/1‑99 (MQ=255)
                  tCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtta                                                                                  >  2:2942735/1‑100 (MQ=255)
                      aaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCgg                                                                              <  1:2342295/100‑1 (MQ=255)
                           ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAAttt                                                                         <  2:2553216/100‑1 (MQ=255)
                            cTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTa                                                                        <  1:1911530/100‑1 (MQ=255)
                             tGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTAc                                                                       >  2:1050091/1‑100 (MQ=255)
                              gTCAAATTGCCTTTTGCCCTGAAAAATGTATAGTCAAATATCTGTTTTAACTAATTGGCGCTGCA‑GTATATGCAAAGCCAATTAGTTAGCGCAATTTa                                                                        >  2:1794042/1‑98 (MQ=255)
                              gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                      >  1:2094821/1‑100 (MQ=255)
                               tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTAGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                      >  1:1829356/1‑99 (MQ=255)
                                cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATTCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                    >  1:2512889/1‑100 (MQ=255)
                                cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                    <  1:3510473/100‑1 (MQ=255)
                                 aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCg                                                                   <  1:1973327/100‑1 (MQ=255)
                                        cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                            <  1:2980687/100‑1 (MQ=255)
                                        cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                            >  1:2475308/1‑100 (MQ=255)
                                         ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGc                                                           <  2:807744/100‑1 (MQ=255)
                                                 tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                                   <  1:3169979/100‑1 (MQ=255)
                                                  taaaaaTGCATAGTTAAATATCTGGTTTAACTAATTGTCGTTGCA‑TTACATCCAACGCCAATTAGTTAGCGGAATTTACGTCGATAATCGCCt                                                         <  2:3226760/92‑1 (MQ=38)
                                                  gAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCa                                                  <  1:2365832/100‑1 (MQ=255)
                                                  gAAAAATACATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATTCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTg                                                        <  2:2605814/94‑1 (MQ=255)
                                                        tGCATAGTCAAATATCTGTTTTAAGTAATTGGCGTTTCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaa                                            <  2:3828335/100‑1 (MQ=255)
                                                          cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGTCGTCATCCCaacg                                          <  1:2657373/100‑3 (MQ=255)
                                                          cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAg                                          >  1:390961/1‑100 (MQ=255)
                                                               tCAAATATCTGTTTTAACTAATTGGCGTTGAA‑GTACATGCAACTCCAATTAGTTAGCGGAATTTACGTCTATACTCGCCTGGCGTCATCCCAAAGCGTTg                                     <  2:3840129/100‑1 (MQ=255)
                                                                cAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                    >  2:1483359/1‑100 (MQ=255)
                                                                cAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                    <  2:1920135/100‑1 (MQ=255)
                                                                 aaaTATCTGTTTTAACTAATTTGATTTGCA‑TTACATTCAACTCTAATTAGTTAGCGATATCTACCTCGTTATTCGCCTTGCGTCATCCCAAAGCGTTGc                                    <  2:2435904/99‑1 (MQ=255)
                                                                  aaTATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                                  >  2:14198/1‑100 (MQ=255)
                                                                    tatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                                >  1:3310404/1‑100 (MQ=255)
                                                                    tatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                                <  1:2648255/100‑1 (MQ=255)
                                                                     atCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaaa                               <  2:62664/100‑1 (MQ=255)
                                                                     atCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCCTaaa                               <  2:1438722/100‑1 (MQ=255)
                                                                      tCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                              <  2:1590406/100‑1 (MQ=255)
                                                                      tCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                              >  2:1787346/1‑100 (MQ=255)
                                                                       cTGTTTTAACTAATTGGCGTTGCA‑GTACATGCCACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaaa                               <  2:1709627/98‑1 (MQ=255)
                                                                       cTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCTAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATa                             <  1:3206634/100‑1 (MQ=255)
                                                                       cTGTTTTAACTAATTGGCCTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTAGCGTCATCCCAAAGCGTTGCTTAAATa                             >  2:2923542/1‑100 (MQ=255)
                                                                           gctAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGt                         >  1:647368/3‑100 (MQ=25)
                                                                            ttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                        <  1:348664/100‑1 (MQ=255)
                                                                            ttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCATTGCTTAAATACCGtt                        >  2:1299792/1‑100 (MQ=255)
                                                                             tAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGAGCCATCCCAAAGCGTTGCTTAAATACCGTTg                       >  2:2351434/1‑100 (MQ=255)
                                                                              aaCTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGAGTTGCTTACATACCGTTGa                      >  2:2801816/1‑100 (MQ=255)
                                                                              aaCTAATTGGCGTTGCA‑GTACATGCAACGCCAATAAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                      >  2:513812/1‑100 (MQ=255)
                                                                               aCTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                     >  1:2298724/1‑100 (MQ=255)
                                                                                  aaTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                  >  2:1043437/1‑100 (MQ=255)
                                                                                   aTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                 <  1:3514410/100‑1 (MQ=255)
                                                                                    ttGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg                >  2:1605941/1‑100 (MQ=255)
                                                                                    ttGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg                <  1:1882292/100‑1 (MQ=255)
                                                                                    ttGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg                <  2:1520361/100‑1 (MQ=255)
                                                                                     tGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa               >  1:609161/1‑100 (MQ=38)
                                                                                     tGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa               >  2:3372533/1‑100 (MQ=255)
                                                                                      ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGACTTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACAGTTGAAAAATGa               >  1:1426553/1‑99 (MQ=37)
                                                                                      ggCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc              >  2:2685382/1‑100 (MQ=255)
                                                                                      ggCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc              >  1:1123029/1‑100 (MQ=255)
                                                                                       gCGTTGCA‑GTACATGCAACGCCTATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAct             >  2:1552850/1‑100 (MQ=255)
                                                                                       gCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc              <  1:1185357/99‑1 (MQ=255)
                                                                                        cGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc            >  1:3534961/1‑100 (MQ=255)
                                                                                          ttGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc            <  1:1857662/98‑1 (MQ=255)
                                                                                            gCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActctcat        >  1:3163908/1‑97 (MQ=255)
                                                                                            gCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt        >  2:2244040/1‑100 (MQ=255)
                                                                                             cA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc       <  2:803388/100‑1 (MQ=255)
                                                                                              a‑GTACGTGCAACGCGAATTAGTTAGCGGAATATACGTCGAGACTCGACTGACATCAGCCCAAAGCGTAGCATGAAGACCGTGGAAAAATGACTCTgatga      >  2:2372254/1‑96 (MQ=255)
                                                                                              a‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAATGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa      <  2:2094802/100‑1 (MQ=255)
                                                                                              a‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa      <  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)
                                                                                                    |                                                                                                  
CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAA  >  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: ‑