Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,797 (C)6→7 coding (1314/3978 nt) ydeK ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,7911.C97.4% 477.8 / 5.2 151coding (1320/3978 nt)ydeKhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (79/68);  total (81/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 184 total aligned reads displayed.

GCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGATT  >  NC_000913/1594700‑1594877
                                                                                            |                                                                                       
gCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTc                                                                                        <  2:4250325/93‑1 (MQ=255)
  gATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTccc                                                                                      >  2:3625860/1‑93 (MQ=255)
  gATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTccc                                                                                      >  1:2415468/1‑93 (MQ=255)
  gATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTccc                                                                                      >  1:2070178/1‑93 (MQ=255)
  gATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTccc                                                                                      >  2:2399596/1‑93 (MQ=255)
  gATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTccc                                                                                      >  2:2816221/1‑93 (MQ=255)
  gATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcac                                                                                      >  2:3127029/1‑91 (MQ=255)
    tACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTccccc                                                                                   >  1:3578888/1‑93 (MQ=255)
      cGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTccccccc                                                                                 <  2:6660428/93‑2 (MQ=255)
       gAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcccccccg                                                                                <  2:2901425/93‑3 (MQ=255)
        aaCTTCCCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcccccccgt                                                                               <  2:2480937/93‑4 (MQ=255)
        aaCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcccccccgt                                                                               <  2:4952673/93‑4 (MQ=255)
        aaCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcccccccgt                                                                               >  2:6315237/1‑90 (MQ=255)
        aaCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcccccccgt                                                                               <  1:4985454/93‑4 (MQ=255)
        aaCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcccccccgg                                                                               >  2:1772242/1‑90 (MQ=255)
         aCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcccccccgtg                                                                              <  2:2853663/93‑5 (MQ=255)
         aCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcccccccgtg                                                                              >  1:1002243/1‑89 (MQ=255)
         aCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcccccccgtg                                                                              >  1:5627399/1‑89 (MQ=255)
         aCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTcccccccgtg                                                                              <  1:3465592/93‑5 (MQ=255)
            tACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCgtggtg                                                                            >  1:1242022/1‑93 (MQ=255)
             aCCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGa                                                                           >  2:1921186/1‑93 (MQ=255)
               cGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATg                                                                         >  2:1820947/1‑93 (MQ=255)
               cGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATg                                                                         >  1:4325689/1‑93 (MQ=255)
                gCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGc                                                                        <  1:1968/93‑1 (MQ=255)
                 ccATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCt                                                                       <  2:336416/93‑1 (MQ=255)
                 ccATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCt                                                                       >  2:4517445/1‑93 (MQ=255)
                  cATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTg                                                                      >  2:4044197/1‑93 (MQ=255)
                  cATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTg                                                                      >  1:894714/1‑93 (MQ=255)
                  cATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTg                                                                      >  1:4551453/1‑93 (MQ=255)
                  cATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCGCCCACATTCCCCCCCGTGGTGATGCTg                                                                      <  1:690312/93‑1 (MQ=255)
                    tCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGAc                                                                    <  2:4272045/93‑1 (MQ=255)
                    tCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGAc                                                                    >  1:5942068/1‑93 (MQ=255)
                      gCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACAt                                                                  >  1:4955124/1‑93 (MQ=255)
                       cTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACAtt                                                                 >  2:381671/1‑93 (MQ=255)
                       cTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACAtt                                                                 >  1:4402452/1‑93 (MQ=255)
                       cTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACAtt                                                                 >  1:2718727/1‑93 (MQ=255)
                        tAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTc                                                                <  1:987869/93‑1 (MQ=255)
                        tAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTc                                                                >  1:1743544/1‑93 (MQ=255)
                        tAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTc                                                                >  1:2945001/1‑93 (MQ=255)
                         tatatTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTcc                                                               <  2:1681585/92‑1 (MQ=255)
                         aaTATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTcc                                                               <  2:6164241/93‑1 (MQ=255)
                           tatTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCg                                                             >  2:348608/1‑93 (MQ=255)
                           tatTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGATGACATTCACg                                                             >  1:135512/1‑93 (MQ=255)
                               gAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCcgccgc                                                         >  1:1202547/1‑93 (MQ=255)
                               gAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGc                                                         <  2:1394871/93‑1 (MQ=255)
                                aGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGcc                                                        <  1:1576254/93‑1 (MQ=255)
                                 ggTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCg                                                       >  1:3813634/1‑93 (MQ=255)
                                 ggTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCg                                                       <  2:3781099/93‑1 (MQ=255)
                                 ggTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCGCCGTGGTGATGCTGACATTCCCGACGCCg                                                       >  2:4085463/1‑93 (MQ=255)
                                   tttCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCCCATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGt                                                     <  1:4536204/93‑1 (MQ=255)
                                    ttCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGtt                                                    <  2:6402402/93‑1 (MQ=255)
                                     tCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGttt                                                   <  1:3593661/93‑1 (MQ=255)
                                     tCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGttt                                                   <  2:2445860/93‑1 (MQ=255)
                                      ccGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTtc                                                  <  1:2948276/93‑1 (MQ=255)
                                      ccGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTtc                                                  <  1:2800842/93‑1 (MQ=255)
                                       cGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTtct                                                 <  2:6001113/93‑1 (MQ=255)
                                        gACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTtctc                                                >  1:581094/1‑93 (MQ=255)
                                        gACGCCCGCAACACCAATGGTGGGATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTtctc                                                >  2:3948714/1‑93 (MQ=255)
                                         aCGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTcc                                               >  1:3856101/1‑93 (MQ=255)
                                          cGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCa                                              <  1:1444671/93‑1 (MQ=255)
                                             ccGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATaa                                           >  2:2197133/1‑93 (MQ=255)
                                              cGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAAc                                          <  2:1731537/93‑1 (MQ=255)
                                                  acacCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGATTCTCCATAAACGAg                                      >  1:3938049/1‑93 (MQ=255)
                                                   cacCAATGGTGGGATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGAGATTCCCGACGCGGGGTTCTTCATAAGCGAGt                                     >  1:3790699/1‑93 (MQ=255)
                                                    acCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTg                                    >  1:3974254/1‑93 (MQ=255)
                                                     ccAATGGTGGTATAGACATTCTCATTACCCTCCCACATTCCCCCCCGTGGTGATGCTGAAATTCACGACGCCGGTTTCTCCATAACCGAGTGa                                   <  1:3695014/93‑1 (MQ=255)
                                                       aaTGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAAc                                 <  2:710583/93‑1 (MQ=255)
                                                       aaTGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAAc                                 <  2:1469807/93‑1 (MQ=255)
                                                        aTGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACt                                <  1:5506170/93‑1 (MQ=255)
                                                        aTGGTGCTATAGACATTCTTATTGCCCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACt                                <  1:4873637/93‑1 (MQ=255)
                                                          ggtggtATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTtc                              <  1:542247/93‑1 (MQ=255)
                                                            tggtATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTtcct                            <  1:4720532/93‑1 (MQ=255)
                                                             ggtATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTtcctc                           <  1:927040/93‑1 (MQ=255)
                                                             ggtATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTtcctc                           >  1:915506/1‑93 (MQ=255)
                                                               tataGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCg                         <  2:5228091/93‑1 (MQ=255)
                                                               tataGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCg                         >  2:6715036/1‑93 (MQ=255)
                                                                ataGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGtt                       <  1:4511562/93‑1 (MQ=255)
                                                                 taGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGCCGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGtt                       >  2:267402/1‑93 (MQ=255)
                                                                 taGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGtt                       >  1:1073254/1‑93 (MQ=255)
                                                                  aGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTg                      >  1:4224835/1‑93 (MQ=255)
                                                                   gACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGc                     <  2:2750495/93‑1 (MQ=255)
                                                                   gACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGc                     <  2:5047941/93‑1 (MQ=255)
                                                                    aCATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCt                    >  2:5497510/1‑93 (MQ=255)
                                                                    aCATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCt                    <  1:2502715/93‑1 (MQ=255)
                                                                      aTTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGac                 >  2:4071074/1‑93 (MQ=255)
                                                                       ttCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGac                 >  1:6461963/1‑93 (MQ=255)
                                                                         cTTATTGACCTCCCCCATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGacaa               <  2:628984/93‑1 (MQ=255)
                                                                         cTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGacaa               >  1:5920910/1‑93 (MQ=255)
                                                                           tATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGacaaca             <  1:3314847/93‑1 (MQ=255)
                                                                            aTTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGAcaacaa            >  1:3038043/1‑93 (MQ=255)
                                                                             ttGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAg           >  2:5353153/1‑93 (MQ=255)
                                                                             ttGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAg           >  2:1349957/1‑93 (MQ=255)
                                                                              tGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGc          <  2:5715437/93‑1 (MQ=255)
                                                                                 ccTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGcccc       >  1:610857/1‑93 (MQ=255)
                                                                                    cccACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGa    >  2:5441846/1‑93 (MQ=255)
                                                                                     ccACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGAt   >  2:3686824/1‑93 (MQ=255)
                                                                                     ccACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGAt   <  1:5851909/93‑1 (MQ=255)
                                                                                      cacaTTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGAtt  <  1:5583844/93‑1 (MQ=255)
                                                                                      cacaTTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGAtt  >  1:671122/1‑93 (MQ=255)
                                                                                      cacaTTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGAtt  <  2:822785/93‑1 (MQ=255)
                                                                                            |                                                                                       
GCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGATT  >  NC_000913/1594700‑1594877

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