Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,039,960 Δ1 bp intergenic (+42/‑36) yihD → / → yihE hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,039,9590A.99.6% 1006.7 / ‑2.6 222intergenic (+41/‑37)yihD/yihEhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (101/119);  minor base C (0/1);  total (102/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047
                                                                                            |                                                                                        
agaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                                          >  2:1868721/1‑93 (MQ=255)
agaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                                          >  2:1556319/1‑93 (MQ=255)
agaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                                          <  1:459208/93‑1 (MQ=255)
 gaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                                         <  1:32370/93‑2 (MQ=255)
 gaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                                         <  1:3113387/93‑2 (MQ=255)
  aagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTACAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGCAATGATatc                                                                                        >  2:350818/1‑91 (MQ=255)
  aagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatc                                                                                        <  1:1441843/93‑3 (MQ=255)
    gaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATAGTGAAATGATatcct                                                                                      >  2:672249/1‑89 (MQ=255)
      cTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTc                                                                                   <  2:436992/93‑1 (MQ=255)
       ttCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCg                                                                                  >  2:4605973/1‑93 (MQ=255)
       ttCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCg                                                                                  >  2:1072336/1‑93 (MQ=255)
        tCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGt                                                                                 >  1:1092677/1‑93 (MQ=255)
         cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTc                                                                                <  2:3650152/93‑1 (MQ=255)
           aaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCAt                                                                              <  1:1907269/93‑1 (MQ=255)
               ggCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGt                                                                          >  2:1178652/1‑93 (MQ=255)
                gcgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTa                                                                         <  1:2553844/93‑1 (MQ=255)
                gcgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTa                                                                         <  1:1010728/93‑1 (MQ=255)
                 cgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTaa                                                                        <  1:3666317/93‑1 (MQ=255)
                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAAt                                                                       >  2:3439597/1‑93 (MQ=255)
                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGATTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAAt                                                                       >  2:4642677/1‑93 (MQ=255)
                   cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATg                                                                      >  2:809726/1‑93 (MQ=255)
                   cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATg                                                                      >  2:2054386/1‑93 (MQ=255)
                    tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGt                                                                     >  2:3929273/1‑93 (MQ=255)
                    tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGt                                                                     <  1:3925958/93‑1 (MQ=255)
                     gTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGGTATCCTTCGTCATTCGTAATGtt                                                                    >  2:3378837/1‑93 (MQ=255)
                     gTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtt                                                                    <  1:4683657/93‑1 (MQ=255)
                     gTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtt                                                                    >  1:4072814/1‑93 (MQ=255)
                       tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                                  <  1:1772842/93‑1 (MQ=255)
                       tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                                  <  2:2242607/93‑1 (MQ=255)
                        aCGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTc                                                                 <  1:3636327/93‑1 (MQ=255)
                          gcgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCg                                                               <  2:3003803/93‑1 (MQ=255)
                          gcgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCg                                                               >  1:2498636/1‑93 (MQ=255)
                            gcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGga                                                             <  1:3866966/93‑1 (MQ=255)
                            gcTCGCGGAGTAATTAAAGAGTAAAAGCTTGTGAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCggt                                                             >  2:2472475/1‑92 (MQ=255)
                             cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            >  2:155781/1‑93 (MQ=255)
                             cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            <  1:2317660/93‑1 (MQ=255)
                              tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatg                                                           <  2:2110954/93‑1 (MQ=255)
                              tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatg                                                           <  2:582699/93‑1 (MQ=255)
                              tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatg                                                           >  1:3652679/1‑93 (MQ=255)
                              tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatg                                                           <  1:2602593/93‑1 (MQ=255)
                                gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatgat                                                         >  1:4017766/1‑93 (MQ=255)
                                gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatgat                                                         <  1:3751562/93‑1 (MQ=255)
                                 cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGatgatg                                                        <  2:2533760/93‑1 (MQ=255)
                                  ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATgg                                                       <  2:3587022/93‑1 (MQ=255)
                                   gtgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                                      <  2:4038406/91‑1 (MQ=255)
                                   gcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                                      >  1:2031230/1‑93 (MQ=255)
                                   gcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                                      <  1:950151/93‑1 (MQ=255)
                                     gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                                    >  2:206841/1‑93 (MQ=255)
                                     gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                                    >  2:714428/1‑93 (MQ=255)
                                       aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                                  >  2:2207741/1‑93 (MQ=255)
                                        aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                 >  2:3102094/1‑93 (MQ=255)
                                        aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                 >  2:2723486/1‑93 (MQ=255)
                                        aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                 <  1:2999822/93‑1 (MQ=255)
                                         atAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                <  1:1388457/92‑1 (MQ=255)
                                          tAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                               >  2:879629/1‑93 (MQ=255)
                                          tAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                               <  2:3264105/93‑1 (MQ=255)
                                          tAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                               >  2:942327/1‑93 (MQ=255)
                                           aaaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                                              <  1:4307694/93‑1 (MQ=255)
                                             agagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACa                                            >  2:232412/1‑93 (MQ=255)
                                             agagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACa                                            >  1:4218649/1‑93 (MQ=255)
                                              gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAg                                           >  2:3996646/1‑93 (MQ=255)
                                              gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAg                                           <  2:2888326/93‑1 (MQ=255)
                                              gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAg                                           <  2:2783910/93‑1 (MQ=255)
                                                 tAAAATCTTATAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                        <  1:3280193/93‑1 (MQ=255)
                                                 tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                        <  2:3324546/93‑1 (MQ=255)
                                                 tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                        <  1:3601065/93‑1 (MQ=255)
                                                 tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                        >  2:901738/1‑93 (MQ=255)
                                                    aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                                     <  2:1354328/93‑1 (MQ=255)
                                                       cTTTTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTAc                                  <  1:2962466/93‑1 (MQ=255)
                                                       cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTat                                  >  1:2766452/1‑92 (MQ=255)
                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAACGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 >  2:3938047/1‑93 (MQ=255)
                                                         tGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                <  1:1270511/93‑1 (MQ=255)
                                                          gTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACttt                               <  1:484052/93‑1 (MQ=255)
                                                             aGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCa                            >  2:1341995/1‑93 (MQ=255)
                                                               cggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGa                          <  2:2606633/93‑1 (MQ=255)
                                                                ggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGac                         <  1:2824530/93‑1 (MQ=255)
                                                                   gcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACt                      <  1:3439765/93‑1 (MQ=255)
                                                                    ccaccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTa                     <  1:1281213/93‑1 (MQ=255)
                                                                    ccaccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTa                     <  2:4712076/93‑1 (MQ=255)
                                                                      accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                   >  2:4323177/1‑93 (MQ=255)
                                                                      accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                   >  1:1656572/1‑93 (MQ=255)
                                                                      accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                   >  1:1993311/1‑93 (MQ=255)
                                                                       ccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTacac                  <  2:3791129/93‑1 (MQ=255)
                                                                        caAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAcc                 <  2:1218678/93‑1 (MQ=255)
                                                                        caAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAcc                 <  2:2181466/93‑1 (MQ=255)
                                                                         aaaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                >  2:490846/1‑93 (MQ=255)
                                                                         aaaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                <  1:620596/93‑1 (MQ=255)
                                                                         aaaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                <  1:4267080/93‑1 (MQ=255)
                                                                           aaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCgg              >  1:562510/1‑93 (MQ=255)
                                                                             tcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGAt            >  2:252181/1‑93 (MQ=255)
                                                                             tcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGAt            <  1:3334793/93‑1 (MQ=255)
                                                                              catcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATa           <  1:676412/93‑1 (MQ=255)
                                                                               atcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATAc          <  1:4505626/93‑1 (MQ=255)
                                                                                tcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATAcc         >  2:1089460/1‑93 (MQ=255)
                                                                                tcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATAcc         >  1:3421827/1‑93 (MQ=255)
                                                                                  gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcat       <  2:1280997/93‑1 (MQ=255)
                                                                                  gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcat       >  1:1053335/1‑93 (MQ=255)
                                                                                   tGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatc      <  1:864781/93‑1 (MQ=255)
                                                                                      aaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATg   <  2:3556070/93‑1 (MQ=255)
                                                                                      aaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATg   <  2:3350774/93‑1 (MQ=255)
                                                                                            |                                                                                        
AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

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