Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,854,698 T→C Y445C (TAC→TGC)  yidJ ← putative sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,854,6980TC84.6% 763.8 / 127.0 332Y445C (TAC→TGC) yidJputative sulfatase
Reads supporting (aligned to +/- strand):  ref base T (17/34);  new base C (122/159);  total (139/193)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 333 total aligned reads displayed.

ACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATC  >  NC_000913/3854608‑3854790
                                                                                          |                                                                                            
aCTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCaa                                                                                            >  1:5673292/1‑93 (MQ=255)
 cTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAAc                                                                                           <  1:7300261/93‑1 (MQ=255)
       aaTAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGa                                                                                     >  1:5276589/1‑93 (MQ=255)
        aTAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGaa                                                                                    <  1:7389831/93‑1 (MQ=255)
         tAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAAc                                                                                   >  2:3141813/1‑93 (MQ=255)
          aGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACg                                                                                  <  2:6465562/93‑1 (MQ=255)
           gCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACgg                                                                                 <  1:5376844/93‑1 (MQ=255)
            ccATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGa                                                                                <  2:360226/93‑1 (MQ=255)
             cATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGAt                                                                               <  2:2685210/93‑1 (MQ=255)
              aTCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATc                                                                              >  2:4600356/1‑93 (MQ=255)
              aTCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATc                                                                              <  1:1371078/93‑1 (MQ=255)
              aTCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATc                                                                              <  2:4336467/93‑1 (MQ=255)
               tCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATcg                                                                             >  1:933305/1‑93 (MQ=255)
               tCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATcg                                                                             >  1:3387607/1‑93 (MQ=255)
                cTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATcgc                                                                            <  2:6026349/93‑1 (MQ=255)
                cTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATcgc                                                                            <  1:4496256/93‑1 (MQ=255)
                cTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATcgc                                                                            >  2:5243811/1‑93 (MQ=255)
                    tggacgtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAAt                                                                        <  2:5750182/93‑1 (MQ=255)
                     ggacgtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAAtt                                                                       <  1:4727958/93‑1 (MQ=255)
                      gacgtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAAttt                                                                      <  1:5620662/93‑1 (MQ=255)
                       acgtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAAtttt                                                                     >  1:6540960/1‑93 (MQ=255)
                         gtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGt                                                                   <  1:3022915/93‑1 (MQ=255)
                          tggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTc                                                                  <  2:6358131/93‑1 (MQ=255)
                          tggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTc                                                                  >  2:169816/1‑93 (MQ=255)
                          tggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTc                                                                  >  1:5097662/1‑93 (MQ=255)
                           ggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTcc                                                                 >  2:7087245/1‑93 (MQ=255)
                            gacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCa                                                                <  2:3958332/93‑1 (MQ=255)
                            gacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCa                                                                >  2:4417745/1‑93 (MQ=255)
                             acgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCAt                                                               >  1:6774021/1‑93 (MQ=255)
                               gAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGt                                                             <  2:2185866/93‑1 (MQ=255)
                                aaaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTa                                                            >  1:1818993/1‑93 (MQ=255)
                                 aaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTaa                                                           <  2:1102216/93‑1 (MQ=255)
                                 aaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTaa                                                           <  1:1044557/93‑1 (MQ=255)
                                 aaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTaa                                                           <  2:4689238/93‑1 (MQ=255)
                                  aCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAAt                                                          <  2:2601068/93‑1 (MQ=255)
                                    gCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATcc                                                        <  1:4978053/93‑1 (MQ=255)
                                    gCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATcc                                                        >  1:3596988/1‑93 (MQ=255)
                                    gCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATcc                                                        >  1:379401/1‑93 (MQ=255)
                                      ccccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCaa                                                      >  1:6033563/1‑93 (MQ=255)
                                      ccccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCaa                                                      <  2:3387430/93‑1 (MQ=255)
                                       cccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAAt                                                     <  1:6142916/93‑1 (MQ=255)
                                       cccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAAt                                                     <  1:5429372/93‑1 (MQ=255)
                                       cccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAAt                                                     >  1:6914486/1‑93 (MQ=255)
                                         cATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATaa                                                   <  1:1983545/93‑1 (MQ=255)
                                         cATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATaa                                                   >  1:5948313/1‑93 (MQ=255)
                                         cATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATaa                                                   <  2:5114805/93‑1 (MQ=255)
                                         cATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATaa                                                   <  1:2853695/93‑1 (MQ=255)
                                          aTCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAg                                                  >  1:6671772/1‑93 (MQ=255)
                                            ccAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGc                                                <  2:692223/93‑1 (MQ=255)
                                             cAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCg                                               <  2:777136/93‑1 (MQ=255)
                                             cAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCg                                               <  2:5536878/93‑1 (MQ=255)
                                              aGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGt                                              >  2:1665794/1‑93 (MQ=255)
                                               gcgcGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTc                                             >  2:2282382/1‑93 (MQ=255)
                                               gcgcGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTc                                             >  1:1196075/1‑93 (MQ=255)
                                               gcgcGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTc                                             >  1:999931/1‑93 (MQ=255)
                                               gcgcGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTc                                             <  1:6632877/93‑1 (MQ=255)
                                                       gTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAttt                                     <  2:2774699/93‑1 (MQ=255)
                                                       gTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAttt                                     >  1:6375234/1‑93 (MQ=255)
                                                       gTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAttt                                     >  2:6207822/1‑93 (MQ=255)
                                                        tGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtttt                                    >  1:4838617/1‑93 (MQ=255)
                                                        tGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtttt                                    <  1:6965240/93‑1 (MQ=255)
                                                         gCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtg                                   >  1:1507570/1‑93 (MQ=255)
                                                          cATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgc                                  >  1:7170419/1‑93 (MQ=255)
                                                           aTCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgct                                 <  1:7318351/93‑1 (MQ=255)
                                                           aTCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgct                                 <  2:3798964/93‑1 (MQ=255)
                                                            tCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgctg                                >  2:3759658/1‑93 (MQ=255)
                                                            tCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgctg                                >  1:43325/1‑93 (MQ=255)
                                                            tCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgctg                                <  2:7249797/93‑1 (MQ=255)
                                                             cTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgctgc                               <  2:3579886/93‑1 (MQ=255)
                                                             cTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgctgc                               <  2:2521055/93‑1 (MQ=255)
                                                              tttACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCg                              >  1:7061771/1‑93 (MQ=255)
                                                               ttACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGa                             >  1:1738571/1‑93 (MQ=255)
                                                                 aCGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAAc                           >  1:3463292/1‑93 (MQ=255)
                                                                   gCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGt                         <  2:7388301/93‑1 (MQ=255)
                                                                    ccACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTc                        <  2:593232/93‑1 (MQ=255)
                                                                    ccACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTc                        <  1:2759404/93‑1 (MQ=255)
                                                                    ccACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTc                        <  1:3245271/93‑1 (MQ=255)
                                                                    ccACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTc                        <  1:4799327/93‑1 (MQ=255)
                                                                    ccACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTc                        <  2:427577/93‑1 (MQ=255)
                                                                      aCGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTg                      <  2:2341881/93‑1 (MQ=255)
                                                                       cGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGc                     <  2:1037770/93‑1 (MQ=255)
                                                                           cGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCaaaa                 >  2:1326307/1‑93 (MQ=255)
                                                                            gCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAAc                <  1:2282195/93‑1 (MQ=255)
                                                                              aGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACgg              >  1:3123301/1‑93 (MQ=255)
                                                                              aGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACgg              <  1:7469482/93‑1 (MQ=255)
                                                                               gACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGa             <  2:1993490/93‑1 (MQ=255)
                                                                               gACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGa             <  2:1517268/93‑1 (MQ=255)
                                                                                aCTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGat            >  1:6294300/1‑93 (MQ=255)
                                                                                 cTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGata           <  1:2462582/93‑1 (MQ=255)
                                                                                 cTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGata           >  1:1310260/1‑93 (MQ=255)
                                                                                   ccATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatc         <  1:5599179/93‑1 (MQ=255)
                                                                                   ccATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatc         >  2:5007283/1‑93 (MQ=255)
                                                                                   ccATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatc         <  1:811282/93‑1 (MQ=255)
                                                                                    cATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatca        <  2:5418561/93‑1 (MQ=255)
                                                                                    cATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatca        >  2:29034/1‑93 (MQ=255)
                                                                                    cATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatca        <  1:2371364/93‑1 (MQ=255)
                                                                                     aTTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatcat       <  2:842652/93‑1 (MQ=255)
                                                                                      ttGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatcatc      >  2:6526703/1‑93 (MQ=255)
                                                                                      ttGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatcatc      >  1:3902370/1‑93 (MQ=255)
                                                                                          tAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATc  >  2:4155866/1‑93 (MQ=255)
                                                                                          |                                                                                            
ACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATC  >  NC_000913/3854608‑3854790

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