Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 502,653 C→T intergenic (‑191/+47) ybaL ← / ← fsr putative transport protein/fosmidomycin resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913502,6530CT99.7% 1109.4 / ‑4.2 318intergenic (‑191/+47)ybaL/fsrputative transport protein/fosmidomycin resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (170/147);  minor base . (1/0);  total (171/147)
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.94e-01

Only 100 of 322 total aligned reads displayed.

ATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATC  >  NC_000913/502563‑502742
                                                                                          |                                                                                         
 ttattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGAc                                                                                        >  1:4107962/1‑93 (MQ=255)
 ttattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGAc                                                                                        >  2:566514/1‑93 (MQ=255)
   attGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTg                                                                                      >  2:7733936/1‑93 (MQ=255)
   attGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTg                                                                                      >  1:2542258/1‑93 (MQ=255)
   attGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTg                                                                                      >  2:6970116/1‑93 (MQ=255)
     tGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACATAGGAAAGGCTTATGGTGACTGTa                                                                                    <  1:1226751/93‑1 (MQ=255)
      gAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaa                                                                                   >  2:8542038/1‑93 (MQ=255)
      gAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaa                                                                                   >  2:4494959/1‑93 (MQ=255)
       aaTTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaa                                                                                  <  1:308821/93‑1 (MQ=255)
        aTTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaaa                                                                                 >  1:1985179/1‑93 (MQ=255)
         ttttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAc                                                                                <  1:4980876/93‑1 (MQ=255)
         ttttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAc                                                                                >  1:6560357/1‑93 (MQ=255)
         ttttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAc                                                                                <  2:7209812/93‑1 (MQ=255)
           ttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGa                                                                              <  1:8535016/93‑1 (MQ=255)
           ttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGa                                                                              <  1:4332854/93‑1 (MQ=255)
           ttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAAGa                                                                              <  2:4765949/93‑1 (MQ=255)
              cAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGc                                                                           <  2:5070796/93‑1 (MQ=255)
              cAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGc                                                                           <  1:3502917/93‑1 (MQ=255)
              cAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGc                                                                           <  1:7401125/93‑1 (MQ=255)
              cAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGc                                                                           >  2:8461940/1‑93 (MQ=255)
               aGAAAAATGGAATTAGCTGCAAAACAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCa                                                                          >  1:3983571/1‑93 (MQ=255)
               aGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCTGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCa                                                                          <  1:8767943/93‑1 (MQ=255)
               aGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCa                                                                          >  2:1989974/1‑93 (MQ=255)
               aGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCa                                                                          <  2:2857218/93‑1 (MQ=255)
                 aaaaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGGTTAACAGAGGGAAGGCTTATGGTGACTGTAGAACGATGCAGc                                                                        >  1:8256173/1‑93 (MQ=255)
                   aaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCa                                                                      >  2:5577472/1‑93 (MQ=255)
                    aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaa                                                                     >  2:254610/1‑93 (MQ=255)
                     aTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaaa                                                                    >  1:8364110/1‑93 (MQ=255)
                      tGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAg                                                                   <  1:8793836/93‑1 (MQ=255)
                      aggAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAg                                                                   <  1:7508302/92‑1 (MQ=255)
                       ggAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGt                                                                  <  2:1276443/93‑1 (MQ=255)
                        gAATTAGCTGCAAAAAAAGCACGGATCGGAAAGCAGAAGGCTGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGtt                                                                 <  1:1785468/93‑1 (MQ=255)
                        gAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGtt                                                                 <  1:8052807/93‑1 (MQ=255)
                          aTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTg                                                               >  1:3785756/1‑93 (MQ=255)
                            tAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGcc                                                             <  2:6047502/93‑1 (MQ=255)
                             aGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCt                                                            <  2:2563371/93‑1 (MQ=255)
                              gCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCtt                                                           >  1:2582597/1‑93 (MQ=255)
                                tGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTg                                                         >  1:4550850/1‑93 (MQ=255)
                                  cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                                       <  2:2097028/93‑1 (MQ=255)
                                  cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                                       <  2:3642426/93‑1 (MQ=255)
                                  cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                                       <  2:620480/93‑1 (MQ=255)
                                   aaaaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATCAAGCCAAAGTTTGCCTTTGGCt                                                      >  2:2800511/1‑93 (MQ=255)
                                    aaaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTg                                                     <  1:705727/93‑1 (MQ=255)
                                     aaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGc                                                    >  1:5018230/1‑93 (MQ=255)
                                       aaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCaa                                                  >  1:156799/1‑93 (MQ=255)
                                        aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAAt                                                 <  1:1052556/93‑1 (MQ=255)
                                        aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAAt                                                 <  2:6139468/93‑1 (MQ=255)
                                         aGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATg                                                >  1:5439691/1‑93 (MQ=255)
                                         aGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATg                                                <  1:2883626/93‑1 (MQ=255)
                                         aGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATg                                                >  2:3546862/1‑93 (MQ=255)
                                            aCGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaaa                                             >  2:3323528/1‑93 (MQ=255)
                                            aCGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaaa                                             >  1:6582453/1‑93 (MQ=255)
                                             cGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAAt                                            >  2:3372481/1‑93 (MQ=255)
                                             cGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAAt                                            <  1:888722/93‑1 (MQ=255)
                                             cGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAAt                                            <  1:776401/93‑1 (MQ=255)
                                              ggATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATc                                           >  2:4311699/1‑93 (MQ=255)
                                              ggATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATc                                           <  1:8555954/93‑1 (MQ=255)
                                               gATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTgcaatgcaataa                                          >  2:4250862/1‑91 (MQ=255)
                                               gATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCa                                          >  2:8719410/1‑93 (MQ=255)
                                                aTCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAg                                         <  1:4808712/93‑1 (MQ=255)
                                                aTCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAg                                         >  2:6407735/1‑93 (MQ=255)
                                                 tCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGt                                        >  2:7933796/1‑93 (MQ=255)
                                                  cGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTc                                       >  1:6191745/1‑93 (MQ=255)
                                                   ggAAAGCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCt                                      >  1:7366574/1‑93 (MQ=255)
                                                   ggAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCt                                      >  1:277759/1‑93 (MQ=255)
                                                    gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCtt                                     >  2:8292716/1‑93 (MQ=255)
                                                    gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCtt                                     >  1:5209781/1‑93 (MQ=255)
                                                    gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCtt                                     >  2:3735564/1‑93 (MQ=255)
                                                     aaaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTCTGCCTTTGGCTGCAATGAAATCAGTCttt                                    <  2:8624086/93‑1 (MQ=255)
                                                      aaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTa                                   >  2:5382073/1‑93 (MQ=255)
                                                       aaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTAt                                  >  2:1491158/1‑93 (MQ=255)
                                                       aaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTAt                                  >  1:2981811/1‑93 (MQ=255)
                                                        aCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATg                                 >  2:3144334/1‑93 (MQ=255)
                                                         cAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTATTTATGc                                <  2:5717756/93‑1 (MQ=255)
                                                          aGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGcc                               >  1:3385774/1‑93 (MQ=255)
                                                           gAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCg                              >  1:4696644/1‑93 (MQ=255)
                                                           gAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCg                              >  2:8117514/1‑93 (MQ=255)
                                                           gAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCg                              >  1:793723/1‑93 (MQ=255)
                                                            aaGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCgg                             <  1:5094182/93‑1 (MQ=255)
                                                              ggCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGtt                           <  1:8206978/93‑1 (MQ=255)
                                                                cGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTAAGCCGGTTAt                         <  1:5711250/93‑1 (MQ=255)
                                                                 ggTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATc                        <  2:4081628/93‑1 (MQ=255)
                                                                  gTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATca                       <  2:5490414/93‑1 (MQ=255)
                                                                   tttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcag                      >  2:6764090/1‑93 (MQ=255)
                                                                         agagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGa                >  2:3067349/1‑93 (MQ=255)
                                                                         agagGAAAGGCTTATGGTGACTGTAAAACGAAGCAGCCAAAGGTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGa                >  2:1174169/1‑93 (MQ=255)
                                                                          gagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGaa               >  1:3904910/1‑93 (MQ=255)
                                                                          gagGAAAGGCTTATGGTGACTGTAAAACGACGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGaa               >  2:4530216/1‑93 (MQ=255)
                                                                            ggAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAata             >  1:3216928/1‑93 (MQ=255)
                                                                            ggAAAGCCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAata             >  2:7616627/1‑93 (MQ=255)
                                                                               aaGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATgg          >  2:1027353/1‑93 (MQ=255)
                                                                                aGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGt         <  1:4537937/93‑1 (MQ=255)
                                                                                aGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGt         >  1:528608/1‑93 (MQ=255)
                                                                                    ttATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAAc     <  1:2375092/93‑1 (MQ=255)
                                                                                    ttATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAAc     >  1:5928524/1‑93 (MQ=255)
                                                                                    ttATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAAc     >  1:6717157/1‑93 (MQ=255)
                                                                                     tATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACa    <  1:2160967/93‑1 (MQ=255)
                                                                                     tATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACa    <  2:3871314/93‑1 (MQ=255)
                                                                                      aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACAt   >  1:711701/1‑93 (MQ=255)
                                                                                       tGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATc  <  2:1764479/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATC  >  NC_000913/502563‑502742

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