Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,043,237 A→G intergenic (+28/+13) b3861 → / ← yihG putative GTP‑binding protein/putative endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,043,2370AG99.8% 1330.1 / ‑5.4 420intergenic (+28/+13)b3861/yihGputative GTP‑binding protein/putative endonuclease
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (209/209);  total (211/209)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 438 total aligned reads displayed.

TTTAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCATT  >  NC_000913/4043147‑4043326
                                                                                          |                                                                                         
 ttAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGc                                                                                        <  1:1178393/93‑1 (MQ=255)
 ttAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGc                                                                                        <  2:5891646/93‑1 (MQ=255)
 ttAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGc                                                                                        <  2:4812384/93‑1 (MQ=255)
 ttAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGc                                                                                        <  1:4958808/93‑1 (MQ=255)
  tAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCg                                                                                       <  2:3770337/93‑1 (MQ=255)
  tAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCg                                                                                       <  1:7183904/93‑1 (MQ=255)
  tAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCg                                                                                       <  1:4342465/93‑1 (MQ=255)
      ggACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTa                                                                                   <  2:846281/93‑1 (MQ=255)
       gACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGGGCGTaa                                                                                  >  2:3126461/1‑93 (MQ=255)
       gACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTaa                                                                                  <  2:2967290/93‑1 (MQ=255)
        aCAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTaaa                                                                                 <  1:4629695/93‑1 (MQ=255)
          aGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTTATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAac                                                                               >  2:6892900/1‑93 (MQ=255)
            gagaTGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAactt                                                                             >  1:874371/1‑93 (MQ=255)
             agaTGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAactta                                                                            <  2:5395861/93‑1 (MQ=255)
               aTGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAacttact                                                                          <  1:1864506/93‑1 (MQ=255)
                tGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAacttactt                                                                         <  2:3705749/93‑1 (MQ=255)
                 gAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAActtactta                                                                        <  1:7505500/93‑1 (MQ=255)
                    cGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTc                                                                     <  1:3226714/93‑1 (MQ=255)
                    cGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTc                                                                     <  2:3336297/93‑1 (MQ=255)
                     gaagaagaaTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCt                                                                    >  1:6239429/1‑93 (MQ=255)
                      aagaagaaTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTg                                                                   >  2:7778573/1‑93 (MQ=255)
                      aagaagaaTTTATGTCTCGTGCCGGACGTTTTGGTCATTAAAATACTAACTCGCCTCTGCTCATTCAGGGGCGAGGAAAACTTACTTATTCTg                                                                   >  1:4894512/1‑93 (MQ=255)
                       agaagaaTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGa                                                                  <  2:4592808/93‑1 (MQ=255)
                       agaagaaTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGCGGCGAGTAAAACTTACTTATTCTgc                                                                  >  1:2162040/1‑92 (MQ=255)
                       agaagaaTTTATGTCTCGTGCCGGACGTTTTGGTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGa                                                                  >  2:2356337/1‑93 (MQ=255)
                        gaagaaTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGAc                                                                 <  2:3566572/93‑1 (MQ=255)
                        gaagaaTTTATGTCTCGTGCCAGACGTTTTGTTCATTAAAATACTCACCCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGAc                                                                 <  1:3714487/93‑1 (MQ=255)
                           gaaTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCt                                                              >  2:397784/1‑93 (MQ=255)
                           gaaTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCt                                                              <  1:5431588/93‑1 (MQ=255)
                             aTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCtgt                                                            >  2:1640792/1‑93 (MQ=255)
                               ttATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGtga                                                          <  1:4261135/93‑1 (MQ=255)
                               ttATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGtga                                                          <  2:5819075/93‑1 (MQ=255)
                               ttATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGtga                                                          <  1:1373261/93‑1 (MQ=255)
                               ttATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGtga                                                          <  1:1536008/93‑1 (MQ=255)
                                tATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAAATTACTTATTCTGACGCTGtgat                                                         >  1:5653223/1‑93 (MQ=255)
                                 aTGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGtgatg                                                        >  1:2836900/1‑93 (MQ=255)
                                 aTGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGtgatg                                                        >  1:5997662/1‑93 (MQ=255)
                                   gTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGAc                                                      >  1:6917803/1‑93 (MQ=255)
                                       cGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATga                                                  <  1:1726644/93‑1 (MQ=255)
                                         tGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATgaga                                                <  1:42941/93‑1 (MQ=255)
                                           ccGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATgagaga                                              <  1:1101026/93‑1 (MQ=255)
                                           ccGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATgagaga                                              <  2:1285478/93‑1 (MQ=255)
                                            cGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGAc                                             <  2:6936178/93‑1 (MQ=255)
                                             ggACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACg                                            <  1:4492732/93‑1 (MQ=255)
                                              gACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGt                                           <  2:2667804/93‑1 (MQ=255)
                                               aCGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTc                                          <  1:4167963/93‑1 (MQ=255)
                                                  ttttGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGa                                       >  1:403085/1‑93 (MQ=255)
                                                   tttGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGaa                                      <  2:4228908/93‑1 (MQ=255)
                                                   tttGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGaa                                      <  2:3203844/93‑1 (MQ=255)
                                                    ttGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAg                                     <  1:5749165/93‑1 (MQ=255)
                                                     tGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGc                                    <  1:7704787/93‑1 (MQ=255)
                                                     tGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGc                                    <  2:5589086/93‑1 (MQ=255)
                                                     tGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGc                                    >  2:7469300/1‑93 (MQ=255)
                                                      gTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGcc                                   >  1:5085499/1‑93 (MQ=255)
                                                       ttCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCg                                  <  2:3924376/93‑1 (MQ=255)
                                                       ttCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCg                                  >  2:7352841/1‑93 (MQ=255)
                                                        tCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCgg                                 <  1:7255952/93‑1 (MQ=255)
                                                        tCATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCgg                                 >  2:7134278/1‑93 (MQ=255)
                                                         cATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGt                                >  2:519815/1‑93 (MQ=255)
                                                         cATTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGt                                <  1:6098316/93‑1 (MQ=255)
                                                          aTTAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTc                               <  2:4515248/93‑1 (MQ=255)
                                                            tAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCtt                             <  1:6518093/93‑1 (MQ=255)
                                                            tAAAATACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCtt                             >  1:3188692/1‑93 (MQ=255)
                                                             aaaaTACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCttt                            <  2:6408660/93‑1 (MQ=255)
                                                                 tACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTcc                        <  2:4873439/93‑1 (MQ=255)
                                                                 tACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTcc                        <  1:2700765/93‑1 (MQ=255)
                                                                 tACTCACTCGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTcc                        <  2:1788568/93‑1 (MQ=255)
                                                                    tcactcGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGc                     >  2:6539210/1‑93 (MQ=255)
                                                                    tcactcGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGc                     >  1:351221/1‑93 (MQ=255)
                                                                    tcactcGACTCTGCTCATTCAGGGGCGAGGAAAACTTACTTATTCTGACGCTGTGATGACATGAGCGACGTCAGAAGCCGGTCTTTTTCCTgt                     >  1:6708723/1‑92 (MQ=255)
                                                                       ctcGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCcac                  <  1:5386018/93‑1 (MQ=255)
                                                                        tcGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCcaca                 >  2:3026422/1‑93 (MQ=255)
                                                                         cGCCTCTGTTAATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACaa                <  2:7220541/93‑1 (MQ=255)
                                                                         cGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGCCGTCAGAAGCCGGTCTTTTCCCTGCCACaa                <  2:2474566/93‑1 (MQ=255)
                                                                         cGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACaa                <  2:5268259/93‑1 (MQ=255)
                                                                         cGCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACaa                <  1:2587384/93‑1 (MQ=255)
                                                                          gCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAAt               >  1:7431249/1‑93 (MQ=255)
                                                                          gCCTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAAt               >  1:6810829/1‑93 (MQ=255)
                                                                           ccTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATg              >  1:2301381/1‑93 (MQ=255)
                                                                           ccTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATg              >  1:1487081/1‑93 (MQ=255)
                                                                           ccTCTGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATg              >  1:2067374/1‑93 (MQ=255)
                                                                            ctctGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGa             >  2:6203689/1‑93 (MQ=255)
                                                                            ctctGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGa             >  2:7612212/1‑93 (MQ=255)
                                                                             tctGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAg            >  2:2060361/1‑93 (MQ=255)
                                                                             tctGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAg            >  1:7756696/1‑93 (MQ=255)
                                                                             tctGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAg            >  2:6697263/1‑93 (MQ=255)
                                                                             cctGCTCATCCAGGGGCGAGTAAACCTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAg            <  1:64009/92‑1 (MQ=255)
                                                                              ctGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGt           <  2:5706462/93‑1 (MQ=255)
                                                                               tGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGCCGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGtt          >  1:4429810/1‑93 (MQ=255)
                                                                               tGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGtt          <  2:1420857/93‑1 (MQ=255)
                                                                               tGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGtt          >  1:3591555/1‑93 (MQ=255)
                                                                               tGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGtt          >  2:6048773/1‑93 (MQ=255)
                                                                               tGCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGACGGCCTTTTTCCTGCCACAATGAGtt          >  1:2373764/1‑93 (MQ=255)
                                                                                gCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTc         <  1:6325470/93‑1 (MQ=255)
                                                                                gCTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTc         <  1:5561066/93‑1 (MQ=255)
                                                                                 cTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCa        >  1:4712000/1‑93 (MQ=255)
                                                                                 cTCATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCa        <  1:5174993/93‑1 (MQ=255)
                                                                                   cATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGc      >  1:3843480/1‑93 (MQ=255)
                                                                                   cATTCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGc      <  2:1178813/93‑1 (MQ=255)
                                                                                      tCAGGGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGCCAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCAt   <  2:2155006/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTTAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCATT  >  NC_000913/4043147‑4043326

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