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.C95.3% 330.0 / 8.4 107coding (1320/3978 nt)ydeKhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base C (69/33);  total (71/36)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 117 total aligned reads displayed.

CGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑CCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGATTA  >  NC_000913/1594701‑1594878
                                                                                             |                                                                                        
cGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCACTGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑cc                                                                                       >  1:1101031/1‑93 (MQ=255)
cGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑cc                                                                                       <  1:3484533/93‑1 (MQ=255)
    aCGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑cccccc                                                                                   <  1:405965/93‑1 (MQ=255)
     cGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑ccccccc                                                                                  <  1:131605/93‑2 (MQ=255)
     cGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑ccccccc                                                                                  <  1:2782246/93‑2 (MQ=255)
      gAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑cccccccg                                                                                 >  2:420861/1‑91 (MQ=255)
       aaCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑cccccccgt                                                                                >  2:2409316/1‑90 (MQ=255)
       aaCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑cccccccgt                                                                                <  2:4660217/93‑4 (MQ=255)
        aCTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑cccccccgtg                                                                               >  2:4560891/1‑89 (MQ=255)
          ttACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCgtggt                                                                              >  1:3793823/1‑93 (MQ=255)
           tACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACTAATGGTGGTACAGACATTCTTATTGACCTCCCACATTCCCCCCCgtggtg                                                                             >  1:4546826/1‑93 (MQ=255)
           tACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTAGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCgtggtg                                                                             >  2:4095259/1‑93 (MQ=255)
            tccGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGa                                                                            <  1:3658976/92‑1 (MQ=255)
             ccGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGAt                                                                           >  1:2251899/1‑93 (MQ=255)
               gCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGc                                                                         >  2:3977346/1‑93 (MQ=255)
               gCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGc                                                                         >  2:3977345/1‑93 (MQ=255)
               gCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGc                                                                         >  2:2538895/1‑93 (MQ=255)
               gCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGc                                                                         <  1:1106184/93‑1 (MQ=255)
               gCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGc                                                                         >  2:1039160/1‑93 (MQ=255)
               gCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGc                                                                         >  1:3835740/1‑93 (MQ=255)
               gCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGc                                                                         >  1:3384680/1‑93 (MQ=255)
               gCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGc                                                                         >  1:245552/1‑93 (MQ=255)
               cccATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATC‑‑CCCCCCGTGGTGATGCt                                                                        <  1:534990/92‑1 (MQ=255)
                 cATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTg                                                                       >  1:3738877/1‑93 (MQ=255)
                   tCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑CCCCCCGTGGTGATGCTGACa                                                                    >  1:3135204/1‑93 (MQ=255)
                    cGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACa                                                                    <  2:2652326/93‑1 (MQ=255)
                     gCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACAt                                                                   <  1:4249415/93‑1 (MQ=255)
                     gCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACAt                                                                   >  1:1713222/1‑93 (MQ=255)
                       tAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTc                                                                 >  1:2556760/1‑93 (MQ=255)
                       tAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTc                                                                 <  1:1457095/93‑1 (MQ=255)
                        aaTATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTcc                                                                >  1:445088/1‑93 (MQ=255)
                        aaTATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTcc                                                                <  1:4005977/93‑1 (MQ=255)
                        aaTATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTcc                                                                <  1:2826552/93‑1 (MQ=255)
                         atatTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTTACCTCTCACATTT‑‑CCCCCCGTGGTGATGCTGACATTccc                                                               >  2:4218750/1‑93 (MQ=255)
                         atatTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTccc                                                               <  2:689587/93‑1 (MQ=255)
                         atatTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTccc                                                               >  2:2517192/1‑93 (MQ=255)
                         atatTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTccc                                                               >  2:2723786/1‑93 (MQ=255)
                         atatTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTccc                                                               >  2:2939075/1‑93 (MQ=255)
                          tatTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCg                                                              <  1:1841339/93‑1 (MQ=255)
                          tatTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCg                                                              <  1:675285/93‑1 (MQ=255)
                            ttGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGAc                                                            <  2:2236506/93‑1 (MQ=255)
                            ttGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGAc                                                            <  1:1875194/93‑1 (MQ=255)
                             tGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACg                                                           <  2:2010878/93‑1 (MQ=255)
                                ggTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCg                                                        <  2:3109295/93‑1 (MQ=255)
                                ggTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑CCCCCCGTGGTGATGCTGACATTCCCGACGCCgg                                                       <  1:3996049/93‑1 (MQ=255)
                                 gTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCgg                                                       >  2:3556316/1‑93 (MQ=255)
                                      cGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTtct                                                  >  2:2492878/1‑93 (MQ=255)
                                       gACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTtct                                                  <  1:4704221/93‑1 (MQ=255)
                                       gACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTtctc                                                 >  1:2193898/1‑93 (MQ=255)
                                          gCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCAt                                              >  2:565259/1‑93 (MQ=255)
                                          gCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCAt                                              >  2:3419137/1‑93 (MQ=255)
                                           cccGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATa                                             >  2:1850826/1‑93 (MQ=255)
                                            ccGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATaa                                            >  2:3051221/1‑93 (MQ=255)
                                            ccGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATaa                                            >  2:1729076/1‑93 (MQ=255)
                                             cGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAAc                                           >  1:4582063/1‑93 (MQ=255)
                                             cGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAAc                                           >  1:875643/1‑93 (MQ=255)
                                               cAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCg                                         >  1:618343/1‑93 (MQ=255)
                                                aaCACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGa                                        >  1:58105/1‑93 (MQ=255)
                                                 acacCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAg                                       >  2:3932105/1‑93 (MQ=255)
                                                 acacCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAg                                       >  2:594756/1‑93 (MQ=255)
                                                 acacCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAg                                       >  1:4625797/1‑93 (MQ=255)
                                                    ccAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGa                                    >  1:3299044/1‑93 (MQ=255)
                                                    ccAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGa                                    >  1:3299043/1‑93 (MQ=255)
                                                    ccAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGa                                    >  1:1284020/1‑93 (MQ=255)
                                                     cAATGGTGTTATAGACATTCTTAGTGTCCTCCCACCTTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGaa                                   <  1:65896/93‑1 (MQ=255)
                                                     cAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGaa                                   <  2:2391866/93‑1 (MQ=255)
                                                      aaTGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAAc                                  >  2:1298304/1‑93 (MQ=255)
                                                            ggtATAGACATTCTTATTGCCCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTtcctc                            <  2:1822105/93‑1 (MQ=255)
                                                            ggtATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTtcctc                            >  2:3851231/1‑93 (MQ=255)
                                                            ggtATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTtcctc                            >  2:764187/1‑93 (MQ=255)
                                                            ggtATAGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTtcctc                            >  2:1912261/1‑93 (MQ=255)
                                                              tataGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCg                          >  2:3501021/1‑93 (MQ=255)
                                                               ataGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGt                         <  1:1331908/93‑1 (MQ=255)
                                                               ataGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGt                         >  2:4215610/1‑93 (MQ=255)
                                                                taGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGtt                        >  1:3170796/1‑93 (MQ=255)
                                                                 aGACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTg                       >  2:2810636/1‑93 (MQ=255)
                                                                  gACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGc                      >  2:3201559/1‑93 (MQ=255)
                                                                  gACATTCTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGc                      >  2:3714787/1‑93 (MQ=255)
                                                                        cTTATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGacaa                >  1:306008/1‑93 (MQ=255)
                                                                         ttATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGacaac               <  1:2914461/93‑1 (MQ=255)
                                                                         ttATTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGacaac               >  2:3565595/1‑93 (MQ=255)
                                                                           tttGACCCCCCACAACC‑‑CCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCAAAAACCGAGTGACCTTCCTCCGTTGCTGAcaacaa             <  2:3931070/92‑1 (MQ=255)
                                                                           aTTGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGAcaacaa             <  2:3930913/93‑1 (MQ=255)
                                                                            ttGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAg            >  2:4035371/1‑93 (MQ=255)
                                                                             tGACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGc           <  2:3069867/93‑1 (MQ=255)
                                                                              gACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGcc          >  2:1175502/1‑93 (MQ=255)
                                                                              gACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGcc          <  2:3041384/93‑1 (MQ=255)
                                                                              gACCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGcc          >  1:3221090/1‑93 (MQ=255)
                                                                               aCCTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGccc         >  1:3365198/1‑93 (MQ=255)
                                                                                ccTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGcccc        >  1:2639183/1‑93 (MQ=255)
                                                                                ccTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGcccc        <  2:2709960/93‑1 (MQ=255)
                                                                                 cTCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCt       <  1:768685/93‑1 (MQ=255)
                                                                                  tCCCACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTg      >  1:3232887/1‑93 (MQ=255)
                                                                                   cccACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGa     >  2:3061908/1‑93 (MQ=255)
                                                                                   cccACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGa     >  2:2768317/1‑93 (MQ=255)
                                                                                   cccACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGa     <  2:3677816/93‑1 (MQ=255)
                                                                                    ccACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGAt    >  1:943680/1‑93 (MQ=255)
                                                                                    ccACATTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGAt    >  1:2416644/1‑93 (MQ=255)
                                                                                     cacaTTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGGTTCTCCATAACCGAATGAACTTCCTCCGTTGCTGACAACAAGCCCCTGCtt   >  2:1710747/1‑93 (MQ=255)
                                                                                      acaTTCCCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGATTa  >  1:1398673/1‑93 (MQ=255)
                                                                                             |                                                                                        
CGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATT‑‑CCCCCCGTGGTGATGCTGACATTCCCGACGCCGGTTTCTCCATAACCGAGTGAACTTCCTCCGTTGCTGACAACAAGCCCCTGATTA  >  NC_000913/1594701‑1594878

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