Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 318,936 C→T G106D (GGC→GAC)  ykgC ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913318,9360CT87.6% 877.5 / 74.0 307G106D (GGC→GAC) ykgCputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (20/18);  new base T (125/144);  total (145/162)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

Only 100 of 313 total aligned reads displayed.

CGGTATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGTATG  >  NC_000913/318844‑319025
                                                                                            |                                                                                         
  gTATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCg                                                                                         <  1:111535/93‑1 (MQ=255)
  gTATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCg                                                                                         <  1:5838316/93‑1 (MQ=255)
  gTATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCg                                                                                         >  2:5469999/1‑93 (MQ=255)
   tATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGt                                                                                        >  1:4538308/1‑93 (MQ=255)
   tATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGt                                                                                        <  2:6159368/93‑1 (MQ=255)
     ttAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGtcgtcg                                                                                      >  1:5431391/1‑93 (MQ=255)
     ttAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGtcgtcg                                                                                      >  2:5353182/1‑93 (MQ=255)
        ataaAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATc                                                                                   >  2:2329239/1‑93 (MQ=255)
         taaAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCa                                                                                  <  1:5104478/93‑1 (MQ=255)
         taaAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCa                                                                                  >  1:6050275/1‑93 (MQ=255)
          aaaaaTTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCAc                                                                                 >  2:249078/1‑93 (MQ=255)
           aaaaTTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACg                                                                                >  1:4214729/1‑93 (MQ=255)
            aaaTTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGt                                                                               <  2:6776383/93‑1 (MQ=255)
               ttttCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGa                                                                            >  2:4669077/1‑93 (MQ=255)
                tttCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGat                                                                           >  1:2060238/1‑93 (MQ=255)
                 ttCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGata                                                                          >  1:6489879/1‑93 (MQ=255)
                  tctcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGatat                                                                         <  2:609484/93‑1 (MQ=255)
                   ctcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATAtt                                                                        <  1:2857524/93‑1 (MQ=255)
                   ctcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATAtt                                                                        <  1:5658195/93‑1 (MQ=255)
                   ctcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATAtt                                                                        >  1:7457415/1‑93 (MQ=255)
                   ctcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATAtt                                                                        <  2:4137953/93‑1 (MQ=255)
                    tcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTg                                                                       >  1:1672635/1‑93 (MQ=255)
                    tcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTg                                                                       <  1:4293300/93‑1 (MQ=255)
                    nngCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTg                                                                       <  2:2741332/91‑1 (MQ=255)
                      gCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTggg                                                                     <  2:3263447/93‑1 (MQ=255)
                      gCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTggg                                                                     >  1:5289577/1‑93 (MQ=255)
                       ccATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGc                                                                    >  2:7070307/1‑93 (MQ=255)
                         aTGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCat                                                                  >  1:3397955/1‑93 (MQ=255)
                         aTGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCat                                                                  >  2:2565894/1‑93 (MQ=255)
                           gAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCatat                                                                <  1:4226328/93‑1 (MQ=255)
                           gAATTTCCAGATTGCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCatat                                                                <  2:4434370/93‑1 (MQ=255)
                             aTTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATcc                                                              >  2:3084458/1‑93 (MQ=255)
                              tttCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCg                                                             <  2:2215974/93‑1 (MQ=255)
                              tttCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCg                                                             >  1:6317735/1‑93 (MQ=255)
                                 ccAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCaa                                                          <  2:558571/93‑1 (MQ=255)
                                 ccAGATTGCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCaa                                                          <  2:5715841/93‑1 (MQ=255)
                                   aGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGa                                                        <  1:1010555/93‑1 (MQ=255)
                                   aGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGa                                                        <  1:6712713/93‑1 (MQ=255)
                                   aGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGa                                                        <  1:4397221/93‑1 (MQ=255)
                                   aGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGa                                                        >  2:4238462/1‑93 (MQ=255)
                                   aGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGa                                                        <  2:4203851/93‑1 (MQ=255)
                                      tttCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTa                                                     <  1:4159428/93‑1 (MQ=255)
                                      tttCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTa                                                     <  2:4616804/93‑1 (MQ=255)
                                      tttCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTa                                                     <  1:816064/93‑1 (MQ=255)
                                      tttCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTa                                                     <  2:4304540/93‑1 (MQ=255)
                                      tttCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTa                                                     <  2:2052347/93‑1 (MQ=255)
                                      tttCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGACCACGTCGATATTGGGCATATCCGCAAGATTa                                                     <  1:4486153/93‑1 (MQ=255)
                                        tCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATg                                                   >  1:2198644/1‑93 (MQ=255)
                                         cccTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGa                                                  <  1:6516598/93‑1 (MQ=255)
                                           cTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGaaa                                                >  1:3792184/1‑93 (MQ=255)
                                             cAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAAt                                              >  2:1442299/1‑93 (MQ=255)
                                             cAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAAt                                              <  1:3285379/93‑1 (MQ=255)
                                              aGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAAtt                                             >  1:5058109/1‑93 (MQ=255)
                                               ggCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTc                                            >  1:5256095/1‑93 (MQ=255)
                                                gCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCt                                           >  1:3991804/1‑93 (MQ=255)
                                                gCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCt                                           >  2:748193/1‑93 (MQ=255)
                                                  cGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCtta                                         <  2:4801919/93‑1 (MQ=255)
                                                   gATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCttat                                        <  2:7192259/93‑1 (MQ=255)
                                                    aTGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCttatt                                       <  2:2316709/93‑1 (MQ=255)
                                                    aTGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCttatt                                       >  2:404492/1‑93 (MQ=255)
                                                    aTGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCttatt                                       >  1:6934236/1‑93 (MQ=255)
                                                      gAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTAc                                     <  2:1666764/93‑1 (MQ=255)
                                                       aaCACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACg                                    <  2:2367206/93‑1 (MQ=255)
                                                         cacGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTa                                  <  1:1534661/93‑1 (MQ=255)
                                                            gCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTaaaa                               >  2:1236576/1‑93 (MQ=255)
                                                              aGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTaaaaaa                             >  1:7021615/1‑93 (MQ=255)
                                                              aGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTaaaaaa                             >  2:5943914/1‑93 (MQ=255)
                                                              aGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTaaaaaa                             <  1:6258621/93‑1 (MQ=255)
                                                               ggCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAAt                            >  2:1093887/1‑93 (MQ=255)
                                                                  tATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTaa                         >  1:1304289/1‑93 (MQ=255)
                                                                   aTGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAAc                        >  2:6472247/1‑93 (MQ=255)
                                                                    tGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAAcc                       >  2:5972225/1‑93 (MQ=255)
                                                                    tGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAAcc                       <  2:5140852/93‑1 (MQ=255)
                                                                    tGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAAcc                       <  2:4929973/93‑1 (MQ=255)
                                                                    tGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAAcc                       >  2:7343390/1‑93 (MQ=255)
                                                                    tGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAAcc                       <  2:862024/93‑1 (MQ=255)
                                                                    tGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAAcc                       >  2:2507353/1‑93 (MQ=255)
                                                                     gATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCa                      >  1:3618407/1‑93 (MQ=255)
                                                                      attattGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCAc                     <  2:2902399/93‑1 (MQ=255)
                                                                       ttattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACt                    <  1:4834853/93‑1 (MQ=255)
                                                                       ttattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACt                    <  2:6616475/93‑1 (MQ=255)
                                                                       ttattGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACt                    >  1:551152/1‑93 (MQ=255)
                                                                        tattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACtt                   >  2:7454590/1‑93 (MQ=255)
                                                                        tattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACtt                   >  1:1938201/1‑93 (MQ=255)
                                                                        tattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACtt                   <  2:3402810/93‑1 (MQ=255)
                                                                         attGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTc                  <  1:5351706/93‑1 (MQ=255)
                                                                          ttGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCa                 >  2:3977557/1‑93 (MQ=255)
                                                                          ttGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCa                 <  2:362481/93‑1 (MQ=255)
                                                                             aTAAACTACGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCAttt              <  2:1875471/93‑1 (MQ=255)
                                                                                aaCTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTa           <  1:1623464/93‑1 (MQ=255)
                                                                                aaCTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTa           >  1:993567/1‑93 (MQ=255)
                                                                                 aCTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTAc          <  1:2635377/93‑1 (MQ=255)
                                                                                 aCTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTAc          <  1:3053345/93‑1 (MQ=255)
                                                                                  cTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACg         >  1:7378097/1‑93 (MQ=255)
                                                                                     cGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTg      >  1:1723071/1‑93 (MQ=255)
                                                                                      gCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGt     <  2:2812241/93‑1 (MQ=255)
                                                                                      gCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGt     >  1:7270513/1‑93 (MQ=255)
                                                                                      gCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGt     >  2:3748748/1‑93 (MQ=255)
                                                                                        cTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGTAt   <  2:6365720/93‑1 (MQ=255)
                                                                                         tGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGTATg  >  2:3482161/1‑93 (MQ=255)
                                                                                            |                                                                                         
CGGTATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGTATG  >  NC_000913/318844‑319025

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