Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,854,389 G→A G116G (GGG→GGA yidI → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,854,3890GA99.2% 899.2 / ‑3.7 242G116G (GGG→GGAyidIhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (91/148);  total (94/148)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 248 total aligned reads displayed.

GGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGACCGGCG  >  NC_000913/3854299‑3854477
                                                                                          |                                                                                        
ggCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAAc                                                                                        <  1:3957180/93‑1 (MQ=255)
ggCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAAc                                                                                        <  2:417914/93‑1 (MQ=255)
ggCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAAc                                                                                        <  2:3583892/93‑1 (MQ=255)
ggCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAAc                                                                                        >  1:544217/1‑93 (MQ=255)
 gCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACg                                                                                       >  1:2807967/1‑93 (MQ=255)
  ccATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGc                                                                                      <  2:4447612/93‑1 (MQ=255)
  ccATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGc                                                                                      <  2:5697880/93‑1 (MQ=255)
    aTTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTg                                                                                    >  2:1541867/1‑93 (MQ=255)
      tgctgcGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGc                                                                                  <  1:119130/93‑1 (MQ=255)
      tgctgcGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGc                                                                                  <  1:682344/93‑1 (MQ=255)
       gctgcGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCt                                                                                 <  2:4239/93‑1 (MQ=255)
       gctgcGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCt                                                                                 <  1:2674639/93‑1 (MQ=255)
       gctgcGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCt                                                                                 <  2:2811346/93‑1 (MQ=255)
       gctgcGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCt                                                                                 >  1:5953302/1‑93 (MQ=255)
           cgcTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTcc                                                                             <  1:1265180/93‑1 (MQ=255)
            gcTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCa                                                                            >  2:3554485/1‑93 (MQ=255)
            gcTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCa                                                                            >  2:1750662/1‑93 (MQ=255)
              tACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGCACGCTGGCTTTCCATa                                                                          >  2:2508282/1‑93 (MQ=255)
               aCTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATAc                                                                         >  1:5718336/1‑93 (MQ=255)
               aCTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATAc                                                                         <  2:6393868/93‑1 (MQ=255)
                cTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATAcc                                                                        >  2:6337261/1‑93 (MQ=255)
                 tttGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCt                                                                       <  1:2915632/93‑1 (MQ=255)
                 tttGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCt                                                                       >  2:4639164/1‑93 (MQ=255)
                  ttGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCtt                                                                      <  1:6179178/93‑1 (MQ=255)
                    gcgcATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGc                                                                    <  1:6193327/93‑1 (MQ=255)
                      gcATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCtt                                                                  >  1:4989073/1‑93 (MQ=255)
                      gcATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCtt                                                                  <  2:1402906/93‑1 (MQ=255)
                      gcATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCtt                                                                  >  1:1471584/1‑93 (MQ=255)
                         tttATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCttttt                                                               <  2:5585039/93‑1 (MQ=255)
                         tttATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCttttt                                                               <  1:5626871/93‑1 (MQ=255)
                         attATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCttttt                                                               >  2:2384113/2‑93 (MQ=255)
                               ggCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATc                                                         <  2:1976482/93‑1 (MQ=255)
                               ggCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATc                                                         <  2:217284/93‑1 (MQ=255)
                               ggCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATc                                                         >  1:2718775/1‑93 (MQ=255)
                                  ggTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGc                                                      <  2:3189962/93‑1 (MQ=255)
                                  ggTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCCTTCCATACCTTGCTTTTTGGTATCGGc                                                      >  1:4841216/1‑93 (MQ=255)
                                   gTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCa                                                     <  1:3699081/93‑1 (MQ=255)
                                   gTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCa                                                     >  1:6408294/1‑93 (MQ=255)
                                   gTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCa                                                     <  2:4957598/93‑1 (MQ=255)
                                   gTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCa                                                     >  2:4685594/1‑93 (MQ=255)
                                    tATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCAt                                                    >  2:2977315/1‑93 (MQ=255)
                                    tATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCAt                                                    <  2:3967298/93‑1 (MQ=255)
                                     aTGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATc                                                   >  1:2077916/1‑93 (MQ=255)
                                       gCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATATCTTGCTTTTTGGTATCGGCATCgt                                                 >  1:4623142/1‑93 (MQ=255)
                                        cGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTCTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCgtg                                                <  1:3768527/93‑1 (MQ=255)
                                         gTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCgtgt                                               <  1:2312035/93‑1 (MQ=255)
                                          tAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGtgtg                                              >  1:2491790/1‑93 (MQ=255)
                                            aGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTgca                                            <  1:1141183/93‑1 (MQ=255)
                                            aGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTgca                                            >  1:2406199/1‑93 (MQ=255)
                                            aGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTgca                                            >  1:2246233/1‑93 (MQ=255)
                                            aGGAAAATCGCTGGGGCATTCGCGGCGCCCTGGTTTTTGGCGGCGGAACGCTGGCTTTTCATACCTTGCTTTTTGGTATCGGCATCGTGTgca                                            >  2:4477287/1‑93 (MQ=255)
                                              gAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTgcagc                                          <  1:1340320/93‑1 (MQ=255)
                                                  aTCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCAttt                                      >  1:879771/1‑93 (MQ=255)
                                                  aTCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCAttt                                      >  1:584291/1‑93 (MQ=255)
                                                   tCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCAtttt                                     >  2:5091005/1‑93 (MQ=255)
                                                    cGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTg                                    >  2:3038272/1‑93 (MQ=255)
                                                       tGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTg                                 >  2:5893962/1‑93 (MQ=255)
                                                       tGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTg                                 >  1:193607/1‑93 (MQ=255)
                                                        ggggCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGa                                <  2:5023864/93‑1 (MQ=255)
                                                        ggggCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGa                                <  1:910288/93‑1 (MQ=255)
                                                         gggCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGat                               >  1:4239984/1‑93 (MQ=255)
                                                          ggCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGata                              <  2:3819359/93‑1 (MQ=255)
                                                           gCATTCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGatat                             >  2:1007315/1‑93 (MQ=255)
                                                               tCGCGGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCt                         <  1:5909878/93‑1 (MQ=255)
                                                                cgcgGCGCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCtt                        <  1:3559451/93‑1 (MQ=255)
                                                                 gcggcgCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTa                       <  1:2635986/93‑1 (MQ=255)
                                                                 gcggcgCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTa                       >  2:5432981/1‑93 (MQ=255)
                                                                  cggcgCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTAt                      <  1:5239833/93‑1 (MQ=255)
                                                                   ggcgCGCTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATc                     <  1:516338/93‑1 (MQ=255)
                                                                     cgcgcTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCtt                   <  1:5937103/93‑1 (MQ=255)
                                                                     cgcgcTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCtt                   <  1:3857074/93‑1 (MQ=255)
                                                                      gcgcTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCttt                  <  1:2003234/93‑1 (MQ=255)
                                                                      gcgcTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCttt                  <  1:5783527/93‑1 (MQ=255)
                                                                      gcgcTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCttt                  <  2:5160182/93‑1 (MQ=255)
                                                                      gcgcTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCttt                  <  2:664800/93‑1 (MQ=255)
                                                                      gcgcTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCttt                  <  2:4097741/93‑1 (MQ=255)
                                                                       cgcTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCtttt                 >  2:2594999/1‑93 (MQ=255)
                                                                         cTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttcc               >  1:3322931/1‑93 (MQ=255)
                                                                         cTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttcc               <  1:1004159/93‑1 (MQ=255)
                                                                         cTGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttcc               <  2:6154095/93‑1 (MQ=255)
                                                                          tGGTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttcct              <  2:1384476/93‑1 (MQ=255)
                                                                           ggTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttcctt             <  1:4355083/93‑1 (MQ=255)
                                                                           ggTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttcctt             >  2:2699437/1‑93 (MQ=255)
                                                                            gTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttccttc            >  2:810504/1‑93 (MQ=255)
                                                                            gTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttccttc            <  1:4523457/93‑1 (MQ=255)
                                                                            gTTTTTGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttccttc            >  1:5549198/1‑93 (MQ=255)
                                                                             tttttGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttccttcc           <  2:3375987/93‑1 (MQ=255)
                                                                             tttttGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttccttcc           <  2:6004996/93‑1 (MQ=255)
                                                                             tttttGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTttccttcc           <  1:1712347/93‑1 (MQ=255)
                                                                              ttttGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTtccttcct          <  2:4835780/93‑1 (MQ=255)
                                                                              ttttGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTtccttcct          <  1:1576927/93‑1 (MQ=255)
                                                                                 tGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGAc       <  1:2995901/93‑1 (MQ=255)
                                                                                 tGGCGGCGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGAc       <  1:3272566/93‑1 (MQ=255)
                                                                                  ggcggcGGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGAcc      >  2:4818752/1‑93 (MQ=255)
                                                                                   gcggcgGAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGACcg     >  1:5611534/1‑93 (MQ=255)
                                                                                    cggcggAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGACcgg    <  1:905039/93‑1 (MQ=255)
                                                                                     ggcggAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGACcggc   <  2:2136623/93‑1 (MQ=255)
                                                                                     ggcggAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGACcggc   >  2:1173138/1‑93 (MQ=255)
                                                                                      gcggAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGACcggcg  <  2:1815309/93‑1 (MQ=255)
                                                                                      gcggAACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGACcggcg  <  2:1551038/93‑1 (MQ=255)
                                                                                          |                                                                                        
GGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCTTTTCCTTCCTGACCGGCG  >  NC_000913/3854299‑3854477

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