Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,008,480 +A coding (637/654 nt) yedK → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,008,4761.A99.0% 708.3 / ‑1.7 196coding (633/654 nt)yedKhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (90/104);  total (91/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

AGCGGTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTATTA  >  NC_000913/2008384‑2008563
                                                                                              |                                                                                       
aGCGGTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGtt                                                                                          <  2:377657/93‑1 (MQ=255)
 gCGGTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTa                                                                                        >  2:1597086/1‑93 (MQ=255)
 gCGGTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGGGGGTAATGTTa                                                                                        >  2:2309133/1‑93 (MQ=255)
  cGGTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTaa                                                                                       <  1:1686473/93‑1 (MQ=255)
    gTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTaaaa                                                                                     >  1:4243489/1‑93 (MQ=255)
    gTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTaaaa                                                                                     >  1:7877913/1‑93 (MQ=255)
    gTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTaaaa                                                                                     <  1:4487410/93‑1 (MQ=255)
     tAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTaaaaa                                                                                    <  1:7367234/93‑2 (MQ=255)
     tAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTaaaaa                                                                                    <  1:4090806/93‑2 (MQ=255)
     tAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTaaaaa                                                                                    <  1:5829730/93‑2 (MQ=255)
       aagaagCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTaaaaacc                                                                                  <  1:4140573/93‑4 (MQ=255)
          aagCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAgg                                                                                >  2:5695586/1‑93 (MQ=255)
           agCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAggg                                                                               >  1:975423/1‑93 (MQ=255)
            gCCTCAGAAATCGCGGCCAGCGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGgc                                                                              >  2:7605669/1‑93 (MQ=255)
              cTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGgcgc                                                                            <  2:5086405/93‑1 (MQ=255)
              cTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGgcgc                                                                            >  1:306027/1‑93 (MQ=255)
              cTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGgcgc                                                                            <  1:6863016/93‑1 (MQ=255)
               tCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGcgcg                                                                           >  1:3433497/1‑93 (MQ=255)
               tCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGcgcg                                                                           <  2:5993254/93‑1 (MQ=255)
                 aGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGa                                                                         >  2:3798504/1‑93 (MQ=255)
                 aGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGa                                                                         <  1:6008020/93‑1 (MQ=255)
                 aGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGa                                                                         <  1:6510277/93‑1 (MQ=255)
                  gAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTACCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAg                                                                        <  2:5386575/93‑1 (MQ=255)
                   aaaTCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGt                                                                       <  1:4769348/93‑1 (MQ=255)
                   aaaTCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGt                                                                       <  1:7559997/93‑1 (MQ=255)
                     aTCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTa                                                                     >  1:5703048/1‑93 (MQ=255)
                       cgcgGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAAt                                                                   <  2:4254540/93‑1 (MQ=255)
                       cgcgGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAAt                                                                   <  1:3247640/93‑1 (MQ=255)
                        gcgGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAAtt                                                                  >  1:1093314/1‑93 (MQ=255)
                        gcgGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAAtt                                                                  >  2:2064289/1‑93 (MQ=255)
                         cgGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTc                                                                 <  1:6156893/93‑1 (MQ=255)
                          ggCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCa                                                                >  1:7205098/1‑93 (MQ=255)
                          ggCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCa                                                                >  2:6831698/1‑93 (MQ=255)
                           gCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCaa                                                               >  2:1394311/1‑93 (MQ=255)
                           gCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCaa                                                               <  2:1326220/93‑1 (MQ=255)
                            ccAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAAc                                                              <  1:3958686/93‑1 (MQ=255)
                             cAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAAcc                                                             <  2:2177440/93‑1 (MQ=255)
                             cAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAAcc                                                             >  2:7680228/1‑93 (MQ=255)
                              aGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCt                                                            >  1:2435795/1‑93 (MQ=255)
                              aGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCt                                                            <  1:59189/93‑1 (MQ=255)
                              aGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGCATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAAccg                                                            >  2:4546254/1‑92 (MQ=255)
                               gTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTg                                                           >  2:1640723/1‑93 (MQ=255)
                               gTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTg                                                           >  2:7369032/1‑93 (MQ=255)
                                  gCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGttt                                                        >  2:4751464/1‑93 (MQ=255)
                                  gCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGttt                                                        >  1:1068340/1‑93 (MQ=255)
                                   cTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTg                                                       >  2:5919993/1‑93 (MQ=255)
                                    tgtgTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGa                                                      >  2:6948094/1‑93 (MQ=255)
                                    tgtgTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGa                                                      <  1:4986424/93‑1 (MQ=255)
                                    tgtgTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGa                                                      >  2:2427616/1‑93 (MQ=255)
                                     gtgtTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGAt                                                     <  2:1527513/93‑1 (MQ=255)
                                      tgtTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATa                                                    <  2:1470354/93‑1 (MQ=255)
                                       gtTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATaa                                                   <  1:6304972/93‑1 (MQ=255)
                                          tcAGGAACCCAGGTCTCCTGGCCGCCGGTAACGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGcc                                                <  1:1208818/92‑1 (MQ=255)
                                          ccAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGcc                                                <  1:2932070/93‑1 (MQ=255)
                                           cAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCt                                               <  1:754122/93‑1 (MQ=255)
                                           cAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGAGCGGAGTTAATTCAACCTGTTTGATAAGCCt                                               <  2:3753787/93‑1 (MQ=255)
                                             gCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTgg                                             <  2:1226513/93‑1 (MQ=255)
                                             gCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTgg                                             >  2:541979/1‑93 (MQ=255)
                                              gaaaCCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGga                                            <  1:5770771/92‑1 (MQ=255)
                                               aaaCCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGgag                                           >  2:7201753/1‑93 (MQ=255)
                                                aaCCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGgaga                                          >  1:7807090/1‑93 (MQ=255)
                                                 aCCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGat                                         >  1:6676510/1‑93 (MQ=255)
                                                    aGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGattat                                      <  1:4382348/93‑1 (MQ=255)
                                                       tctcCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGa                                   <  2:6970607/93‑1 (MQ=255)
                                                        ctcCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGAt                                  <  1:4345477/93‑1 (MQ=255)
                                                        ctcCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGAt                                  <  1:5873814/93‑1 (MQ=255)
                                                         tcCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATc                                 <  2:4846739/93‑1 (MQ=255)
                                                          ccTGGCACGCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATca                                <  2:3088733/93‑1 (MQ=255)
                                                             ggCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATcaaca                             >  2:7504990/1‑93 (MQ=255)
                                                              gCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATcaacaa                            >  1:6265383/1‑93 (MQ=255)
                                                                ccccGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCACCCTGTTTGATAAGCCTGGAGATTATTGATCAACAATa                          <  2:649071/92‑1 (MQ=255)
                                                                aCCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATa                          <  1:3033910/93‑1 (MQ=255)
                                                                aCCCGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATAAACAATa                          <  1:944470/93‑1 (MQ=255)
                                                                   cGGTATCGCGCGCGGTGGGTAATGGTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTg                       >  2:3734380/1‑93 (MQ=255)
                                                                   cGGTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTg                       >  1:5114315/1‑93 (MQ=255)
                                                                    ggTATCGCGCGCGGTGTGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGc                      <  2:1393184/93‑1 (MQ=255)
                                                                    ggTATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGc                      <  1:1897584/93‑1 (MQ=255)
                                                                      tATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGt                    <  1:673993/93‑1 (MQ=255)
                                                                      tATCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGt                    <  2:4684258/93‑1 (MQ=255)
                                                                       aTCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTc                   >  2:5333230/1‑93 (MQ=255)
                                                                       aTCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTc                   <  1:2131675/93‑1 (MQ=255)
                                                                        tCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCa                  >  2:5158428/1‑93 (MQ=255)
                                                                        tCGCGCGCGGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCa                  <  1:5574642/93‑1 (MQ=255)
                                                                         cgcgcgcgGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCAt                 <  2:3482365/93‑1 (MQ=255)
                                                                              gcgGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaa            <  1:7647718/93‑1 (MQ=255)
                                                                               cgGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTGAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaaa           <  1:6331426/93‑1 (MQ=255)
                                                                               cgGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaaa           >  2:3909558/1‑93 (MQ=255)
                                                                               cgGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaaa           <  1:6051861/93‑1 (MQ=255)
                                                                               cgGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaaa           <  2:1890617/93‑1 (MQ=255)
                                                                               cgGTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaaa           <  1:5535248/93‑1 (MQ=255)
                                                                               cgGTGGATAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaaa           <  1:5261252/93‑1 (MQ=255)
                                                                                 gTGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAAtc         <  2:6226784/93‑1 (MQ=255)
                                                                                  tGGGTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAAtct        >  1:3167439/1‑93 (MQ=255)
                                                                                   gggTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAAtctc       <  1:3808344/93‑1 (MQ=255)
                                                                                   gggTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAAtctc       <  1:4519293/93‑1 (MQ=255)
                                                                                    ggTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATctct      >  1:1254368/1‑93 (MQ=255)
                                                                                     gTAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAACAAATctctt     >  1:2627946/1‑92 (MQ=255)
                                                                                      tAATGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTAt    >  1:1008697/1‑93 (MQ=255)
                                                                                       aaTGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTAtt   <  1:1802262/93‑1 (MQ=255)
                                                                                        aTGTTAAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTATTa  >  1:4327118/1‑93 (MQ=255)
                                                                                              |                                                                                       
AGCGGTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAACCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTATTA  >  NC_000913/2008384‑2008563

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