Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,525 +G intergenic (+92/‑141) b4555 → / → yigM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,5241.G99.6% 962.4 / ‑3.8 256intergenic (+91/‑142)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (126/129);  total (127/129)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01

Only 100 of 273 total aligned reads displayed.

CGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTT  >  NC_000913/4008434‑4008611
                                                                                             |                                                                                        
cGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt                                                                                           <  1:737491/93‑1 (MQ=255)
  ggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg                                                                                        >  2:172112/1‑92 (MQ=255)
         cgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                                                                                 <  1:3177658/93‑1 (MQ=255)
         cgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                                                                                 >  1:5113939/1‑93 (MQ=255)
           acgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                                                                               >  2:515563/1‑93 (MQ=255)
                gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                                                                          >  2:3294343/1‑93 (MQ=255)
                gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                                                                          >  2:608246/1‑93 (MQ=255)
                gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                                                                          >  1:1199901/1‑93 (MQ=255)
                   aTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt                                                                       <  2:2542082/93‑1 (MQ=255)
                   aTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt                                                                       >  2:6487699/1‑93 (MQ=255)
                    tGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                                                      >  2:4729384/1‑93 (MQ=255)
                    tGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                                                      <  2:5953053/93‑1 (MQ=255)
                      cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                                                                    >  1:1743840/1‑93 (MQ=255)
                      cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                                                                    <  1:4629335/93‑1 (MQ=255)
                       tGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg                                                                   >  1:6401371/1‑93 (MQ=255)
                       tGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg                                                                   >  1:3625383/1‑93 (MQ=255)
                       cgTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg                                                                   <  1:506992/92‑1 (MQ=255)
                        gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                                                  <  2:4865780/93‑1 (MQ=255)
                        gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                                                  >  2:1185427/1‑93 (MQ=255)
                         tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                                                                 <  2:3782441/93‑1 (MQ=255)
                         tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                                                                 <  2:3028737/93‑1 (MQ=255)
                          cGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacg                                                                >  1:523193/1‑93 (MQ=255)
                           gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                                                               <  2:5635097/93‑1 (MQ=255)
                            aTGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                                                              >  1:5701826/1‑93 (MQ=255)
                              ggcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc                                                            >  1:4753811/5‑93 (MQ=255)
                               gcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                                                           >  2:452360/4‑93 (MQ=255)
                               gcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                                                           >  1:1820823/4‑93 (MQ=255)
                               gcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                                                           >  1:558129/4‑93 (MQ=255)
                                cggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg                                                          <  1:4029892/91‑1 (MQ=255)
                                cggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg                                                          >  1:300252/3‑93 (MQ=255)
                                cggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg                                                          >  1:2382718/3‑93 (MQ=255)
                                 ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                                                         >  1:1705896/2‑93 (MQ=255)
                                 ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                                                         <  1:2670173/92‑1 (MQ=255)
                                 ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                                                         <  2:3844288/92‑1 (MQ=255)
                                   ggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGc                                                       >  1:2926756/1‑93 (MQ=255)
                                   ggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGc                                                       <  2:2055246/93‑1 (MQ=255)
                                   ggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGc                                                       <  2:3052035/93‑1 (MQ=255)
                                     aaaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTGCTAATGCCGACGACGCGGTGCCg                                                     <  1:4697058/89‑1 (MQ=255)
                                      aaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGc                                                    >  2:3675794/4‑93 (MQ=255)
                                       aggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCa                                                   <  2:3364328/91‑1 (MQ=255)
                                       aggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCa                                                   >  2:1312660/3‑93 (MQ=255)
                                           tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAt                                               <  2:2126216/93‑1 (MQ=255)
                                           tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAt                                               <  1:6465690/93‑1 (MQ=255)
                                            gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc                                              <  1:1378233/93‑1 (MQ=255)
                                            gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc                                              >  2:3730663/1‑93 (MQ=255)
                                             cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt                                             <  1:6020153/93‑1 (MQ=255)
                                             cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt                                             <  1:3570102/93‑1 (MQ=255)
                                             cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt                                             >  1:1503404/1‑93 (MQ=255)
                                                  tGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTa                                        <  2:5587924/93‑1 (MQ=255)
                                                   gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaa                                       >  2:1373889/1‑93 (MQ=255)
                                                    ggCAGTGCGCACGAGCGGCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa                                      <  1:979996/93‑1 (MQ=255)
                                                    ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa                                      >  1:6153694/1‑93 (MQ=255)
                                                     gCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAc                                     >  2:660382/1‑93 (MQ=255)
                                                      cAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAct                                    <  1:3048197/93‑1 (MQ=255)
                                                           cgcaccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTAtt                               <  2:5340615/91‑1 (MQ=255)
                                                           cgcaccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTAtt                               <  2:2328243/91‑1 (MQ=255)
                                                            gcaccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTAttt                              <  2:4233000/92‑1 (MQ=255)
                                                             caccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTa                             >  2:5687610/1‑93 (MQ=255)
                                                              accaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTAt                            <  1:100433/93‑1 (MQ=255)
                                                               ccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATc                           >  1:3834797/1‑93 (MQ=255)
                                                                 aGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGt                         >  1:577187/1‑93 (MQ=255)
                                                                 aGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGt                         <  1:687844/93‑1 (MQ=255)
                                                                 aGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGACGCATTATCTGCGTAAACTCTATTTATCGt                         >  1:2250569/1‑93 (MQ=255)
                                                                  gCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTa                        <  2:1855330/93‑1 (MQ=255)
                                                                   cGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTaa                       <  2:3500554/93‑1 (MQ=255)
                                                                   cGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTaa                       >  1:4310769/1‑93 (MQ=255)
                                                                     tCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                     <  1:4488876/93‑1 (MQ=255)
                                                                      gtGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCg                    <  1:1285537/92‑1 (MQ=255)
                                                                      cTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCg                    >  2:1951092/1‑93 (MQ=255)
                                                                       tGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGt                   <  1:422988/93‑1 (MQ=255)
                                                                        gAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGtt                  >  1:3517674/1‑93 (MQ=255)
                                                                         aaaGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTc                 >  2:1542161/1‑93 (MQ=255)
                                                                         aaaGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTc                 >  2:584949/1‑93 (MQ=255)
                                                                         aaaGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTc                 >  1:162501/1‑93 (MQ=255)
                                                                          aaaaCCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCt                <  2:5325532/90‑1 (MQ=255)
                                                                          aaGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCt                <  1:2129041/93‑1 (MQ=255)
                                                                           aGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCtt               <  2:3935201/93‑1 (MQ=255)
                                                                           aGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCtt               >  2:1511770/1‑93 (MQ=255)
                                                                            gACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCttt              <  1:629456/93‑1 (MQ=255)
                                                                            gACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCttt              >  2:2578301/1‑93 (MQ=255)
                                                                             aCCTTTATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTa             <  1:4923225/93‑1 (MQ=255)
                                                                              ccTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTAt            >  1:5424655/1‑93 (MQ=255)
                                                                              ccTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTAt            <  2:4409643/93‑1 (MQ=255)
                                                                               cTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTAtt           >  1:2760907/1‑93 (MQ=255)
                                                                               cTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTAtt           <  2:3238540/93‑1 (MQ=255)
                                                                                ttCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTAttt          <  2:3410116/93‑1 (MQ=255)
                                                                                ttCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTAttt          <  2:3142386/93‑1 (MQ=255)
                                                                                 tCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTg         >  1:2860219/1‑93 (MQ=255)
                                                                                 tCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTg         >  2:2440865/1‑93 (MQ=255)
                                                                                 tCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTg         <  1:3377039/93‑1 (MQ=255)
                                                                                   aTCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGt       >  2:1657541/1‑93 (MQ=255)
                                                                                   aTCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGt       >  2:6184634/1‑93 (MQ=255)
                                                                                    tCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTc      <  2:817646/93‑1 (MQ=255)
                                                                                    tCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTc      <  2:2171988/93‑1 (MQ=255)
                                                                                    tCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTc      >  2:2896225/1‑93 (MQ=255)
                                                                                      ccGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAg    >  2:876994/1‑93 (MQ=255)
                                                                                       cGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGt   <  2:5453072/93‑1 (MQ=255)
                                                                                       cGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGt   >  1:1104611/1‑93 (MQ=255)
                                                                                       cGAGCTGTAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGt   >  2:4815951/1‑93 (MQ=255)
                                                                                        gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGtt  >  2:6440628/1‑93 (MQ=255)
                                                                                             |                                                                                        
CGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTT  >  NC_000913/4008434‑4008611

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