Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,007,020 T→C C241R (TGC→CGC)  yedE → putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,007,0200TC98.5% 393.8 / ‑4.4 129C241R (TGC→CGC) yedEputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (57/70);  total (59/70)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

GCGCTGCTCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCTGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGCCATC  >  NC_000913/2006930‑2007109
                                                                                          |                                                                                         
gcgcTGCTCACAGCGATGAACCAACCCAAACTGGGGCTGGCAGTGCTGGTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGc                                                                                         >  2:1601119/1‑93 (MQ=255)
gcgcTGCTCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGc                                                                                         <  2:412120/93‑1 (MQ=255)
  gctgctCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCtt                                                                                       >  2:2514038/1‑93 (MQ=255)
  gctgctCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCtt                                                                                       <  2:2162995/93‑1 (MQ=255)
     gctCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTAc                                                                                    >  2:2065190/1‑93 (MQ=255)
     gctCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTAc                                                                                    >  1:4551335/1‑93 (MQ=255)
     gctCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTAc                                                                                    <  1:3142808/93‑1 (MQ=255)
      ctCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACt                                                                                   >  2:4676898/1‑93 (MQ=255)
      ctCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACt                                                                                   <  1:3327950/93‑1 (MQ=255)
      ctCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACt                                                                                   <  1:2330213/93‑1 (MQ=255)
        cacaGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTc                                                                                 >  1:1290969/1‑93 (MQ=255)
          caGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAg                                                                               >  2:1789167/1‑93 (MQ=255)
          caGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAg                                                                               <  1:526406/93‑1 (MQ=255)
            gctaTGAGCCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCg                                                                             <  2:2133724/90‑1 (MQ=255)
            gCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCg                                                                             <  2:2879005/93‑1 (MQ=255)
              gATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGtt                                                                           >  2:4067716/1‑93 (MQ=255)
               aTGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCTGCTTTACTTCAGCGTTc                                                                          >  1:3836262/1‑93 (MQ=255)
                tGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTcc                                                                         <  2:2416597/93‑1 (MQ=255)
                tGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTcc                                                                         <  1:1882330/93‑1 (MQ=255)
                  aaccaaccCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCcgc                                                                       >  2:1437399/1‑93 (MQ=255)
                    ccaaccCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGa                                                                     >  2:4536287/1‑93 (MQ=255)
                      aaccCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGata                                                                   <  1:1617033/93‑1 (MQ=255)
                       accCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGatat                                                                  >  2:914845/1‑93 (MQ=255)
                         ccAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATAtgt                                                                <  2:348098/93‑1 (MQ=255)
                         ccAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATAtgt                                                                <  2:4281169/93‑1 (MQ=255)
                          cAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATAtgtg                                                               <  2:1469406/93‑1 (MQ=255)
                          cAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATAtgtg                                                               <  2:1932292/93‑1 (MQ=255)
                          cAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATAtgtg                                                               <  2:1469403/93‑1 (MQ=255)
                           aaaCTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGCgg                                                              >  1:4326646/1‑93 (MQ=255)
                            aaCTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGa                                                             <  1:3370942/93‑1 (MQ=255)
                             aCTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGAt                                                            <  2:1622051/93‑1 (MQ=255)
                             aCTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGAt                                                            <  2:1240889/93‑1 (MQ=255)
                             aCTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGAt                                                            <  2:4396888/93‑1 (MQ=255)
                              cTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATc                                                           <  2:2782297/93‑1 (MQ=255)
                              cTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATc                                                           >  1:654729/1‑93 (MQ=255)
                              cTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATc                                                           <  1:1293326/93‑1 (MQ=255)
                                 gggCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCAcc                                                        >  2:2417257/1‑93 (MQ=255)
                                   gCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCgg                                                      <  1:3336568/93‑1 (MQ=255)
                                     tGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGAc                                                    >  2:4635986/1‑93 (MQ=255)
                                       gCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGt                                                  >  2:1986019/1‑93 (MQ=255)
                                        cAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTa                                                 >  1:2710103/1‑93 (MQ=255)
                                        cAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTa                                                 <  1:1569801/93‑1 (MQ=255)
                                         aaTGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTAc                                                >  1:2578328/1‑93 (MQ=255)
                                          aTGCTGTTTGGGGTCGGCTTTGGTTTAATGCTTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGCTCACAGGACGTAcc                                               >  2:1219657/1‑93 (MQ=255)
                                          aTGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTAcc                                               >  2:522867/1‑93 (MQ=255)
                                          aTGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTAcc                                               >  2:3495822/1‑93 (MQ=255)
                                          aTGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTAcc                                               >  1:3384628/1‑93 (MQ=255)
                                          aTGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTAcc                                               >  1:3384627/1‑93 (MQ=255)
                                            gCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATTACCGGACGTACCCa                                             <  2:2862986/93‑1 (MQ=255)
                                            gCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCa                                             <  1:2199467/93‑1 (MQ=255)
                                             gtGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCat                                            <  1:2467726/92‑1 (MQ=255)
                                             cTGTTTGTCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCat                                            <  2:4216681/93‑1 (MQ=255)
                                              tGTTTGGNGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCata                                           >  1:2446784/1‑93 (MQ=255)
                                              tGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCata                                           <  1:2986182/93‑1 (MQ=255)
                                               gTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCatat                                          <  1:3083457/93‑1 (MQ=255)
                                                tttGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATg                                         >  2:2647670/1‑93 (MQ=255)
                                                tttGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATg                                         <  1:1683477/93‑1 (MQ=255)
                                                tgtGGCGTCGTCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATg                                         <  2:4297182/91‑1 (MQ=255)
                                                  tGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGc                                       <  1:225919/93‑1 (MQ=255)
                                                      gTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGaaa                                   >  2:885637/1‑93 (MQ=255)
                                                      gTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGaaa                                   <  1:1903732/93‑1 (MQ=255)
                                                       tCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAg                                  >  1:925927/1‑93 (MQ=255)
                                                       tCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAg                                  <  2:2560300/93‑1 (MQ=255)
                                                        cGGCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGc                                 >  2:2967399/1‑93 (MQ=255)
                                                         ggCTTTGGTTTACTGATTGAACGTGTGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCa                                >  1:4199284/1‑93 (MQ=255)
                                                         ggCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCa                                >  2:1783026/1‑93 (MQ=255)
                                                          gCTTTGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCaa                               >  1:1252274/1‑93 (MQ=255)
                                                            tttGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATc                             <  1:3973877/93‑1 (MQ=255)
                                                             ttGGTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCa                            >  1:2885766/1‑93 (MQ=255)
                                                               ggTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCAtt                          >  1:2601928/1‑93 (MQ=255)
                                                                gTTTACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTa                         <  1:3995292/93‑1 (MQ=255)
                                                                 tttACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTAt                        >  1:2967394/1‑93 (MQ=255)
                                                                 tttACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTAt                        <  2:2987196/93‑1 (MQ=255)
                                                                  ttACTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATc                       <  1:2387575/93‑1 (MQ=255)
                                                                    aCTGATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATCGCGAAAGCAATCATTATCgg                     >  2:1874108/1‑93 (MQ=255)
                                                                       gATTGAACGTGCGCATATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTAt                  <  2:1584968/93‑1 (MQ=255)
                                                                       gATTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTAt                  >  1:674857/1‑93 (MQ=255)
                                                                        aTTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATg                 <  2:385621/93‑1 (MQ=255)
                                                                        aTTGAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATg                 <  2:2308409/93‑1 (MQ=255)
                                                                           gAACGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCg              <  2:1846320/93‑1 (MQ=255)
                                                                            aaCGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCgg             >  2:256852/1‑93 (MQ=255)
                                                                            aaCGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCgg             <  2:830601/93‑1 (MQ=255)
                                                                            aaCGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCgg             <  2:4449376/93‑1 (MQ=255)
                                                                             aCGTGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCggg            >  1:1416087/1‑92 (MQ=255)
                                                                                tGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAg         >  2:1253464/1‑93 (MQ=255)
                                                                                tGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAg         <  2:1990252/93‑1 (MQ=255)
                                                                                tGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAg         >  1:1602902/1‑93 (MQ=255)
                                                                                tGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAg         <  2:2177002/93‑1 (MQ=255)
                                                                                tGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAg         <  1:2963695/93‑1 (MQ=255)
                                                                                tGCGCAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAg         >  1:4133815/1‑93 (MQ=255)
                                                                                   gcAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGc      >  1:578559/1‑93 (MQ=255)
                                                                                   gcAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGc      >  2:596914/1‑93 (MQ=255)
                                                                                   gcAAATCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGc      >  2:896683/1‑93 (MQ=255)
                                                                                     aaaTCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGCCa    <  2:3440127/93‑1 (MQ=255)
                                                                                      aaTCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGCCAt   <  2:3793591/93‑1 (MQ=255)
                                                                                      aaTCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGCCAt   <  2:3672482/93‑1 (MQ=255)
                                                                                      aaTCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGCCAt   >  1:3212334/1‑93 (MQ=255)
                                                                                       aTCTGCTTTACTTCAGCGTTCCCCGATATGTGCAACCCCGGACGTACCCATTTGGCGAAAGCAATGATTATGGGTATAGCGGTGAGTGCCATc  >  2:4486090/1‑93 (MQ=255)
                                                                                       aTCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGCCATc  <  2:835544/93‑1 (MQ=255)
                                                                                       aTCCGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGCCATc  <  2:1754310/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGCTGCTCACAGCGATGAACCAACCCAAACTGGGGCTGGCAATGCTGTTTGGCGTCGGCTTTGGTTTACTGATTGAACGTGCGCAAATCTGCTTTACTTCAGCGTTCCGCGATATGTGGATCACCGGACGTACCCATATGGCGAAAGCAATCATTATCGGTATGGCGGTGAGTGCCATC  >  NC_000913/2006930‑2007109

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