Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,491 2 bp→CG intergenic (+58/‑174) 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,4910GC99.7% 1241.1 / NA 370intergenic (+58/‑175)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (171/198);  total (171/199)
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
*NC_0009134,008,4920CG99.5% 1224.9 / NA 365intergenic (+59/‑174)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base G (166/197);  total (166/199)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 370 total aligned reads displayed.

CTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACT  >  NC_000913/4008411‑4008577
                                                                                     ||                                                                                      
cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAg                                                                                  >  2:1382978/1‑93 (MQ=255)
 tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                 <  1:1202663/93‑1 (MQ=255)
 tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                 >  2:6209325/1‑93 (MQ=255)
 tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                 >  1:7078464/1‑93 (MQ=255)
 tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                 >  2:5455356/1‑93 (MQ=255)
 tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                 <  1:209109/93‑1 (MQ=255)
 tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                 <  2:1959311/93‑1 (MQ=255)
  gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCTTTATGGGCAGTGCGCACCAGCg                                                                                >  1:1112201/1‑93 (MQ=255)
  gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                                                                >  2:1138424/1‑93 (MQ=255)
    aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                                              >  2:3427937/1‑93 (MQ=255)
    aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                                              <  2:471445/93‑1 (MQ=255)
                  tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                                                >  2:6948155/5‑93 (MQ=255)
                  tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                                                >  2:2939129/5‑93 (MQ=255)
                    tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                                              <  1:2764591/91‑1 (MQ=255)
                       tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                                           <  1:5312574/93‑1 (MQ=255)
                        gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg                                                          <  2:1858739/93‑1 (MQ=255)
                           gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTg                                                      >  1:4515664/1‑93 (MQ=255)
                             gATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgga                                                    >  1:3510725/1‑91 (MQ=255)
                              aTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                                   <  2:4884585/93‑4 (MQ=255)
                              aTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                                   <  1:273441/93‑4 (MQ=255)
                              aTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                                   <  1:2399878/93‑4 (MQ=255)
                               tGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                                                  >  1:447193/1‑89 (MQ=255)
                                 aaCGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                                                 >  2:750328/1‑93 (MQ=255)
                                  acgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                                                <  2:2110599/93‑1 (MQ=255)
                                  acgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                                                >  1:7429522/1‑93 (MQ=255)
                                   cgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                                               <  2:3815664/93‑1 (MQ=255)
                                   cgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                                               <  2:2887079/93‑1 (MQ=255)
                                     acgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                                             >  1:1363953/1‑93 (MQ=255)
                                     acgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                                             >  2:3238232/1‑93 (MQ=255)
                                       gCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                                           <  1:719952/93‑1 (MQ=255)
                                        ccGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                                          <  2:5608195/93‑1 (MQ=255)
                                        ccGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                                          <  2:248902/93‑1 (MQ=255)
                                         cGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc                                         >  1:5644213/1‑93 (MQ=255)
                                         cGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc                                         <  2:7378296/93‑1 (MQ=255)
                                         cGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc                                         <  1:3942131/93‑1 (MQ=255)
                                          gagaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                                        <  1:5861038/90‑1 (MQ=255)
                                          gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                                        >  2:150043/1‑93 (MQ=255)
                                           aaaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTa                                       <  1:2436251/93‑1 (MQ=255)
                                            aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                      <  2:4088021/93‑1 (MQ=255)
                                            aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                      >  2:6796432/1‑93 (MQ=255)
                                            aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                      <  1:966786/93‑1 (MQ=255)
                                            aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                      <  1:414147/93‑1 (MQ=255)
                                             aTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt                                     <  1:3754984/93‑1 (MQ=255)
                                              tGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                    >  2:4941722/1‑93 (MQ=255)
                                              tGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                    <  1:6470282/93‑1 (MQ=255)
                                              tGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                    >  2:3168847/1‑93 (MQ=255)
                                               gCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGGGCTGGAAGTGATTGGTACTAATGc                                   <  1:659224/93‑1 (MQ=255)
                                               gCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGGGCTGGAAGTGATTGGTACTAATGc                                   <  2:726571/93‑1 (MQ=255)
                                               gCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                                   <  1:517814/93‑1 (MQ=255)
                                                 tGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg                                 >  2:6077482/1‑93 (MQ=255)
                                                  gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                <  1:3119614/93‑1 (MQ=255)
                                                  gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                <  2:3631821/93‑1 (MQ=255)
                                                  gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                <  2:6870522/93‑1 (MQ=255)
                                                   tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                               <  2:1086426/93‑1 (MQ=255)
                                                   tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                               >  1:6710388/1‑93 (MQ=255)
                                                   tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                               <  2:7208058/93‑1 (MQ=255)
                                                   tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                               <  2:6492572/93‑1 (MQ=255)
                                                   tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                               <  1:3602871/93‑1 (MQ=255)
                                                     gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                             >  1:1777695/1‑93 (MQ=255)
                                                     gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                             <  2:6032356/93‑1 (MQ=255)
                                                      aTGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                            >  2:6328395/1‑93 (MQ=255)
                                                      aTGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                            >  1:5502038/1‑93 (MQ=255)
                                                        ggcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc                          <  1:6057983/89‑1 (MQ=255)
                                                          cggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg                        <  1:4896224/91‑1 (MQ=255)
                                                           ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                       <  1:755291/92‑1 (MQ=255)
                                                           ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                       >  2:6118904/2‑93 (MQ=255)
                                                             ggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGc                     <  2:5218696/93‑1 (MQ=255)
                                                                aaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGc                  >  1:1685486/4‑93 (MQ=255)
                                                                  ggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAt                <  1:3907554/92‑1 (MQ=255)
                                                                  ggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAt                >  1:160993/2‑93 (MQ=255)
                                                                    cTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGGGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTa              <  1:902764/93‑1 (MQ=255)
                                                                    cTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTa              >  1:2938991/1‑93 (MQ=255)
                                                                     tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGGGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAt             <  1:6248073/93‑1 (MQ=255)
                                                                     tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAt             >  2:4176554/1‑93 (MQ=255)
                                                                     tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAt             <  2:76525/93‑1 (MQ=255)
                                                                     tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAt             <  1:427396/93‑1 (MQ=255)
                                                                     tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAt             <  1:5237510/93‑1 (MQ=255)
                                                                      gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc            <  2:3334585/93‑1 (MQ=255)
                                                                      gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc            <  2:4345596/93‑1 (MQ=255)
                                                                      gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc            >  1:1814614/1‑93 (MQ=255)
                                                                      gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc            <  2:2778616/93‑1 (MQ=255)
                                                                       cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt           <  2:5767344/93‑1 (MQ=255)
                                                                       cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt           <  1:931612/93‑1 (MQ=255)
                                                                        aTTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGGCGACGCGGTGCCGCATTATCTg          >  2:2427925/1‑93 (MQ=255)
                                                                        aTTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTg          >  1:2189860/1‑93 (MQ=255)
                                                                         ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc         >  1:5429835/1‑93 (MQ=255)
                                                                         ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc         >  1:3248900/1‑93 (MQ=255)
                                                                         ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc         <  1:3392610/93‑1 (MQ=255)
                                                                          tATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCg        <  2:3938331/93‑1 (MQ=255)
                                                                          tATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCg        <  1:5772837/93‑1 (MQ=255)
                                                                           aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGt       <  2:1821233/93‑1 (MQ=255)
                                                                            tGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTa      >  2:2171423/1‑93 (MQ=255)
                                                                             gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaa     <  1:6778698/93‑1 (MQ=255)
                                                                              ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa    <  1:2578335/93‑1 (MQ=255)
                                                                              ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa    <  1:809081/93‑1 (MQ=255)
                                                                              ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa    >  2:1339807/1‑93 (MQ=255)
                                                                               gCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAc   <  2:1169844/93‑1 (MQ=255)
                                                                               gCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAc   <  1:1574299/93‑1 (MQ=255)
                                                                               gCAGTGCGCACCAGCGGCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAc   <  1:4742218/93‑1 (MQ=255)
                                                                                cAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAct  <  2:4489685/93‑1 (MQ=255)
                                                                                     ||                                                                                      
CTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACT  >  NC_000913/4008411‑4008577

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