Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,049,079 T→C intergenic (‑94/‑25) ymcD ← / → b4514 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,049,0790TC99.3% 359.7 / ‑4.7 145intergenic (‑94/‑25)ymcD/b4514hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (67/77);  total (67/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 157 total aligned reads displayed.

ATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGC  >  NC_000913/1048989‑1049142
                                                                                          |                                                               
aTGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGccc                                                               >  2:1066812/1‑93 (MQ=255)
aTGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGccc                                                               <  2:1687098/93‑1 (MQ=255)
aTGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGccc                                                               <  2:2490140/93‑1 (MQ=255)
 tGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGcccc                                                              >  2:1208917/1‑93 (MQ=255)
 tGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGcccc                                                              >  2:5231258/1‑93 (MQ=255)
 tGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGcccc                                                              >  2:2474539/1‑93 (MQ=255)
  ggACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCcct                                                             <  1:6009911/93‑1 (MQ=255)
  ggACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCcct                                                             <  2:2835245/93‑1 (MQ=255)
   gACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCcctc                                                            <  2:3024616/93‑1 (MQ=255)
    gcAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCcctcc                                                           <  1:2012322/92‑1 (MQ=255)
    aCAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCcctcc                                                           <  2:459049/93‑1 (MQ=255)
    aCAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCcctcc                                                           <  2:3343322/93‑1 (MQ=255)
    aCAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCcctcc                                                           <  2:3127260/93‑1 (MQ=255)
     cAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCcctcct                                                          <  2:3929690/93‑1 (MQ=255)
      aaTACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTg                                                         >  1:2229709/1‑93 (MQ=255)
      aaTACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTg                                                         >  2:1665726/1‑93 (MQ=255)
       aTACGGGTGATGCTGCCAGCTTCCTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGt                                                        <  1:6190877/93‑1 (MQ=255)
       aTACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGt                                                        >  2:5160180/1‑93 (MQ=255)
        tACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGtt                                                       <  1:1488817/93‑1 (MQ=255)
          cGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCa                                                     <  1:3834882/93‑1 (MQ=25)
          cGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCa                                                     <  1:107749/93‑1 (MQ=21)
           gggAGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTATTTCTATCAGCTGCCCCTCCTGTTCAg                                                    >  1:2515084/1‑93 (MQ=255)
            ggTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGc                                                   <  1:6213005/93‑1 (MQ=18)
            ggTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGc                                                   >  1:4714654/1‑93 (MQ=18)
             gTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCt                                                  >  1:5932353/1‑93 (MQ=18)
             gTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCt                                                  >  1:2988998/1‑93 (MQ=18)
               gATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTAc                                                >  1:365886/1‑93 (MQ=18)
               gATGCTGCCAAATTACTGATTTCGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATTAGCTGCGCCTTCTGTTCAGCTAc                                                >  1:365887/1‑93 (MQ=11)
                aTGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACt                                               >  1:5844838/1‑93 (MQ=18)
                aTGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACt                                               >  1:2554093/1‑93 (MQ=18)
                 tgctgcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTg                                              <  1:37988/93‑1 (MQ=18)
                 tgctgcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTg                                              <  2:5898250/93‑1 (MQ=18)
                  gctgcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGa                                             >  1:1809428/1‑93 (MQ=18)
                  gctgcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGa                                             >  1:6255175/1‑93 (MQ=18)
                    tgcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACg                                           <  2:5089542/93‑1 (MQ=18)
                    tgcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGCTCTCCAGTGGCTTATGTTTCTATCAGCTGCACCTCCTGTTCAGCTCCTGACg                                           >  1:6454924/1‑93 (MQ=11)
                     gcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACgg                                          >  1:765624/1‑93 (MQ=18)
                     gcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACgg                                          <  2:4097250/93‑1 (MQ=18)
                      ccAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACggg                                         <  2:4101237/93‑1 (MQ=18)
                       cANCTTNCTGATTTAGTGTATNNTGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACgggg                                        <  2:2798332/93‑1 (MQ=255)
                       cAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACgggg                                        <  1:5391466/93‑1 (MQ=18)
                        aaCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACcgggg                                       >  1:5931023/1‑92 (MQ=255)
                        aaCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGggt                                       >  1:6284273/1‑93 (MQ=18)
                        aaCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGggt                                       <  2:5490742/93‑1 (MQ=18)
                        aaCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGggt                                       <  2:3344458/93‑1 (MQ=18)
                        aaCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGggt                                       <  2:6286981/93‑1 (MQ=18)
                        aaCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGggt                                       >  2:363926/1‑93 (MQ=18)
                         aCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGggtg                                      >  1:3209757/1‑93 (MQ=18)
                         aCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGggtg                                      >  1:1298404/1‑93 (MQ=18)
                          cTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGggtgg                                     <  1:2220226/93‑1 (MQ=18)
                           ttACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGggtggt                                    <  1:2300858/93‑1 (MQ=18)
                           tcaCTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGggtggt                                    <  2:5239760/91‑1 (MQ=16)
                            tACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGgtggtg                                   <  1:1208339/93‑1 (MQ=18)
                            tACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCACCTCCTGTTCAGCTACTGACGGGgtggtg                                   >  1:907722/1‑93 (MQ=16)
                               tGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGGGCGt                                >  1:4486129/1‑93 (MQ=255)
                                gATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTa                               >  2:917966/1‑93 (MQ=21)
                                  tttAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAAc                             <  1:2783032/93‑1 (MQ=18)
                                   ttAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACg                            <  2:5601839/93‑1 (MQ=18)
                                    tAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACgg                           >  1:613081/1‑93 (MQ=18)
                                    tAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACgg                           >  2:1328690/1‑93 (MQ=21)
                                    tAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACgg                           <  2:4632657/93‑1 (MQ=18)
                                    tAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGATGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACgg                           >  1:1536556/1‑93 (MQ=16)
                                     aGTGTATGATGGTTTTTTTGAGTTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGc                          <  2:5328335/93‑1 (MQ=255)
                                      gtgtATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCa                         >  2:4784994/1‑93 (MQ=18)
                                        gtATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCaaa                       <  1:5820482/93‑1 (MQ=18)
                                        gtATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCaaa                       >  1:2651817/1‑93 (MQ=18)
                                         tATGATGGTGTTTTTGAGTTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCaaaa                      <  2:2986424/93‑1 (MQ=16)
                                         tATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCaaaa                      <  1:3407254/93‑1 (MQ=18)
                                           tgatgGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGc                    <  1:4136198/93‑1 (MQ=18)
                                             atgGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAc                  >  1:2859941/1‑93 (MQ=18)
                                             atgGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAc                  >  1:5175189/1‑93 (MQ=18)
                                             atgGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGGAACGGCAAAAGCAc                  >  2:4410822/1‑93 (MQ=255)
                                             atgGTGTTTTTGAGGTTCTCCAGTGGCGTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAc                  <  2:1926332/93‑1 (MQ=16)
                                              tgGTTTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAcc                 <  1:4321716/93‑2 (MQ=17)
                                              tgGTTTTTGTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAcc                 <  1:4578745/93‑2 (MQ=11)
                                              tgGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCCGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACt                 >  2:1979564/1‑93 (MQ=16)
                                              tgGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACt                 >  1:4955162/1‑93 (MQ=18)
                                              tgGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACt                 >  1:1499644/1‑93 (MQ=18)
                                              tgGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACt                 <  2:798518/93‑1 (MQ=18)
                                               ggTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATTAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGctgtg                >  1:5532822/1‑89 (MQ=11)
                                                gtgtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGc               >  1:2723521/1‑93 (MQ=18)
                                                gtgtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGc               <  1:109714/93‑1 (MQ=18)
                                                gtgtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGc               <  2:3666209/93‑1 (MQ=18)
                                                gtgtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGc               >  2:4462572/1‑93 (MQ=18)
                                                 tgtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGcc              <  1:779911/93‑1 (MQ=18)
                                                 tgtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGcc              <  2:6139153/93‑1 (MQ=18)
                                                 tgtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGcc              <  1:4622045/93‑1 (MQ=18)
                                                   tttttGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCgg            <  2:6287964/93‑1 (MQ=18)
                                                    ttttGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGa           <  2:2283522/93‑1 (MQ=18)
                                                       tGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACAt        >  1:5640177/1‑93 (MQ=21)
                                                        gAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACTGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATc       >  1:3006029/1‑93 (MQ=16)
                                                        gAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTCACGGCAAAAGCACTGCCGGAAATc       >  2:686904/1‑93 (MQ=255)
                                                        gAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATc       >  1:4530596/1‑93 (MQ=18)
                                                        gAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATc       <  1:5511987/93‑1 (MQ=18)
                                                        gAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATc       <  1:4288264/93‑1 (MQ=18)
                                                        gAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATc       <  1:3437491/93‑1 (MQ=18)
                                                         aGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCa      >  1:1485851/1‑93 (MQ=18)
                                                          ggTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAg     <  1:2054317/93‑1 (MQ=18)
                                                          ggTTCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAg     <  2:100245/93‑1 (MQ=18)
                                                             tctcCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAgcgc  <  1:5338901/93‑1 (MQ=14)
                                                                                          |                                                               
ATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGC  >  NC_000913/1048989‑1049142

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