Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 735,021 C→T A50V (GCC→GTC)  ybfC → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913735,0210CT99.6% 864.0 / ‑5.2 244A50V (GCC→GTC) ybfChypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (111/132);  total (112/132)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

Only 100 of 248 total aligned reads displayed.

TGATTGCGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGCCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAATTCACACCAA  >  NC_000913/734929‑735110
                                                                                            |                                                                                         
tGATTGCGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGc                                                                                           >  1:243807/1‑93 (MQ=255)
  aTTGCGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCa                                                                                         <  1:4533737/93‑1 (MQ=255)
   ttGCGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCAt                                                                                        <  1:1203196/93‑1 (MQ=255)
   ttGCGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCAtg                                                                                       <  1:1569282/93‑1 (MQ=255)
    tGCGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCAtg                                                                                       <  1:2711317/93‑1 (MQ=255)
    tGCGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCAtg                                                                                       <  2:2189449/93‑1 (MQ=255)
     gCGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCAtgt                                                                                      <  2:1630516/93‑1 (MQ=255)
     gCGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCAtgt                                                                                      >  1:1877210/1‑93 (MQ=255)
      cGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCAtgtg                                                                                     <  1:2394619/93‑1 (MQ=255)
        aaaTACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGga                                                                                   <  1:1104830/93‑1 (MQ=255)
          aTACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGgaga                                                                                 <  2:1903812/93‑1 (MQ=255)
           tACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGAt                                                                                >  2:2775423/1‑93 (MQ=255)
           tACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGAt                                                                                <  2:151759/93‑1 (MQ=255)
            aCAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATg                                                                               <  2:1831573/93‑1 (MQ=255)
            aCAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATg                                                                               <  1:2194639/93‑1 (MQ=255)
            aCAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATg                                                                               <  2:1038357/93‑1 (MQ=255)
              aaGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGt                                                                             >  2:1585004/1‑93 (MQ=255)
               aGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTa                                                                            >  1:131115/1‑93 (MQ=255)
               aGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTa                                                                            <  1:799835/93‑1 (MQ=255)
                gCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTTAATAATAAAGTCATGTGGAGATGGTAt                                                                           >  2:1220189/1‑93 (MQ=255)
                gCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTAt                                                                           >  1:723159/1‑93 (MQ=255)
                gCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTAt                                                                           >  1:3136595/1‑93 (MQ=255)
                gCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTAt                                                                           >  1:2825199/1‑93 (MQ=255)
                  tAAAGATCATGATTGTTTTAGTATTTTCGCTAAGGCTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCa                                                                         <  1:961077/93‑1 (MQ=255)
                  tAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCa                                                                         >  1:3920185/1‑93 (MQ=255)
                  tAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCa                                                                         >  2:1482774/1‑93 (MQ=255)
                    aaGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGa                                                                       <  2:3034641/93‑1 (MQ=255)
                     aGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGaa                                                                      <  1:984543/93‑1 (MQ=255)
                       aTCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATg                                                                    >  1:1406253/1‑93 (MQ=255)
                       aTCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATg                                                                    >  2:3855282/1‑93 (MQ=255)
                        tCATGATTGTTTTACTCTTTTCGCTAAGGGGACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGa                                                                   >  2:4063240/1‑93 (MQ=255)
                        tCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGa                                                                   >  1:2371182/1‑93 (MQ=255)
                           tGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGa                                                                >  1:2677576/1‑93 (MQ=255)
                           tGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGa                                                                >  2:564559/1‑93 (MQ=255)
                           tGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGGGGAGAGGGTATAAAAATGAgca                                                                >  1:4067799/1‑91 (MQ=255)
                            gATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGAc                                                               >  2:3535750/1‑93 (MQ=255)
                            gATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGAc                                                               >  2:3942551/1‑93 (MQ=255)
                                gTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCg                                                           >  1:4013021/1‑93 (MQ=255)
                                gTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCg                                                           >  1:2580019/1‑93 (MQ=255)
                                gTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCg                                                           <  2:419448/93‑1 (MQ=255)
                                  tttATTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGt                                                         <  2:2669553/93‑1 (MQ=255)
                                  tttACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGt                                                         <  2:4153243/93‑1 (MQ=255)
                                   ttACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTg                                                        <  1:2256507/93‑1 (MQ=255)
                                    tACTATTTTCGCTAAGGGTACTATTTTTTTCGCGTTTCCTGTCTTAATTAATAAAGTCAAGTGGAGATGGTATCAGAATGAGGACATCggggg                                                       >  1:4199108/1‑90 (MQ=255)
                                    tACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGa                                                       >  1:3886650/1‑93 (MQ=255)
                                    tACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGa                                                       >  2:4079010/1‑93 (MQ=255)
                                    tACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGa                                                       >  2:4286970/1‑93 (MQ=255)
                                     aCTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTTAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAg                                                      >  1:1847828/1‑93 (MQ=255)
                                      cTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGt                                                     <  1:2479173/93‑1 (MQ=255)
                                       tATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGta                                                    <  2:2396244/93‑1 (MQ=255)
                                       tATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGta                                                    >  2:2928559/1‑93 (MQ=255)
                                        aTTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGtat                                                   <  1:4635636/93‑1 (MQ=255)
                                        aTTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGtat                                                   <  1:3361621/93‑1 (MQ=255)
                                         ttttCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGtatt                                                  >  2:1304953/1‑93 (MQ=255)
                                           ttCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGtattat                                                <  2:3776370/93‑1 (MQ=255)
                                            tCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTattatt                                               <  2:172451/93‑1 (MQ=255)
                                            tCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTattatt                                               <  2:1385566/93‑1 (MQ=255)
                                             cGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTg                                              >  1:3638345/1‑93 (MQ=255)
                                              gCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTgg                                             <  1:3822180/93‑1 (MQ=255)
                                                tAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCa                                           >  1:143061/1‑93 (MQ=255)
                                                 aaGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCaa                                          >  2:3370593/1‑93 (MQ=255)
                                                  aGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGATATCGGTGAGTATTATTGGCaaa                                         >  2:425011/1‑93 (MQ=255)
                                                  aGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCaaa                                         >  1:914075/1‑93 (MQ=255)
                                                  aGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCaaa                                         >  2:2265149/1‑93 (MQ=255)
                                                   gggTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAAc                                        <  1:3960389/93‑1 (MQ=255)
                                                    ggTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACa                                       <  2:4348685/93‑1 (MQ=255)
                                                     gTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAg                                      >  1:4262884/1‑93 (MQ=255)
                                                        cTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAAc                                   <  1:4090056/93‑1 (MQ=255)
                                                        cTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAAc                                   <  1:2854128/93‑1 (MQ=255)
                                                         tATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACt                                  >  1:2972849/1‑93 (MQ=255)
                                                          aTTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTg                                 <  2:742639/93‑1 (MQ=255)
                                                           ttttttCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTgg                                <  2:903975/93‑1 (MQ=255)
                                                            tttttCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTggg                               >  2:79528/1‑93 (MQ=255)
                                                             ttttCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGGCATCGGTGAGTATTATTGGCAAACAGAACTGGGt                              <  2:4539865/93‑1 (MQ=255)
                                                             ttttCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGt                              <  1:4145744/93‑1 (MQ=255)
                                                             ttttCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGt                              >  1:3171857/1‑93 (MQ=255)
                                                             ttttCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGt                              <  1:2010505/93‑1 (MQ=255)
                                                              tttCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTa                             <  2:2951578/93‑1 (MQ=255)
                                                                tCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTAct                           >  2:4628956/1‑92 (MQ=255)
                                                                tCCGCGTTTCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACa                           <  2:1168220/93‑1 (MQ=255)
                                                                      tttCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTaaa                     >  2:1762570/1‑93 (MQ=255)
                                                                      tttCCTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTaaa                     >  2:3713076/1‑93 (MQ=255)
                                                                          cTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAaata                 >  1:4577653/1‑93 (MQ=255)
                                                                          cTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAaata                 <  2:2888533/93‑1 (MQ=255)
                                                                          cTGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGGATTATTGGCAAACAGAACTGGGTACATGTAAAaata                 >  1:3467709/1‑93 (MQ=255)
                                                                           tGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAaataa                <  1:2948093/93‑1 (MQ=255)
                                                                           tGTCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAaataa                <  1:1161793/93‑1 (MQ=255)
                                                                             tCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAataata              <  2:3746180/93‑1 (MQ=255)
                                                                             tCTTAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAataata              >  1:3760544/1‑93 (MQ=255)
                                                                               ttAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGTGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATaaa            <  2:2297665/93‑1 (MQ=255)
                                                                               ttAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATaaa            <  2:3068466/93‑1 (MQ=255)
                                                                                tAAATAATAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAAt           <  1:1940410/93‑1 (MQ=255)
                                                                                  aataatAAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAATTc         >  2:3150611/1‑93 (MQ=255)
                                                                                   ataataAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAATTca        >  2:1116587/1‑93 (MQ=255)
                                                                                   ataataAAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAATTca        <  1:1244004/93‑1 (MQ=255)
                                                                                    taataaAGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAATTcac       >  2:1317399/1‑93 (MQ=255)
                                                                                        aaaGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAATTCACACCa   >  2:2068397/1‑93 (MQ=255)
                                                                                        aaaGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAATTCACACCa   <  1:3119358/93‑1 (MQ=255)
                                                                                         aaGTCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAATTCACACCaa  >  1:4100810/1‑93 (MQ=255)
                                                                                         aaGCCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAATTCACACCaa  >  1:3668669/1‑93 (MQ=255)
                                                                                            |                                                                                         
TGATTGCGAAATACAAGCTAAAGATCATGATTGTTTTACTATTTTCGCTAAGGGTACTATTTTTTCCGCGTTTCCTGTCTTAAATAATAAAGCCATGTGGAGATGGTATCAGAATGAGGACATCGGTGAGTATTATTGGCAAACAGAACTGGGTACATGTAAAAATAATAAATTCACACCAA  >  NC_000913/734929‑735110

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