Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 2,651,283 G→A G45E (GGA→GAA)  sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_0000912,651,2830GA99.2% 522.7 / NA 131G45E (GGA→GAA) sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (48/82);  total (49/82)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

TTGATGACCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAA  >  AC_000091/2651196‑2651369
                                                                                       |                                                                                      
ttGATGACCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAAc                                                                                      <  1:350101/90‑1 (MQ=255)
ttGATGACCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAAc                                                                                      >  2:1850584/1‑90 (MQ=255)
 tgatgaCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACa                                                                                     <  2:751587/90‑1 (MQ=255)
 tgatgaCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACa                                                                                     >  1:3226189/1‑90 (MQ=255)
  gatgaCCCGGGAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACat                                                                                    >  2:630392/1‑90 (MQ=255)
  gatgaCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACat                                                                                    <  2:3041196/90‑1 (MQ=255)
  gatgaCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACat                                                                                    <  2:1008996/90‑1 (MQ=255)
   atgaCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACata                                                                                   <  1:853924/90‑1 (MQ=255)
    tgaCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACatat                                                                                  <  2:1370618/90‑1 (MQ=255)
    tgaCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACatat                                                                                  <  2:2681583/90‑1 (MQ=255)
    tgaCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACatat                                                                                  <  1:1693562/90‑1 (MQ=255)
    tgaCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACatat                                                                                  >  2:2237579/1‑90 (MQ=255)
      aCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTc                                                                                <  2:1598747/90‑1 (MQ=255)
      aCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTc                                                                                <  1:3021679/90‑1 (MQ=255)
      aCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTc                                                                                >  2:2197879/1‑90 (MQ=255)
      aCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTc                                                                                >  2:2936642/1‑90 (MQ=255)
       cccGGAAATTCAGATTATCCATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTcc                                                                               <  1:2981265/90‑1 (MQ=255)
        ccGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTccc                                                                              <  2:1415715/90‑1 (MQ=255)
         cGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCg                                                                             >  1:3254790/1‑90 (MQ=255)
          ggAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGGTCAGGAGTATCTGAATGAACATATTCCCgg                                                                            >  1:895845/1‑90 (MQ=255)
          ggAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCgg                                                                            <  2:1777585/90‑1 (MQ=255)
          ggAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCgg                                                                            <  1:1330901/90‑1 (MQ=255)
           gAAATTCAGATTATCGATGCCCGCCTGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGgc                                                                           <  1:2392918/90‑1 (MQ=255)
           gAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGgc                                                                           >  2:2036465/1‑90 (MQ=255)
            aaaTTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGgcg                                                                          >  1:987510/1‑90 (MQ=255)
            aaaTTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGgcg                                                                          <  2:488373/90‑1 (MQ=255)
             aaTTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGgcgc                                                                         <  2:1792762/90‑1 (MQ=255)
              aTTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCa                                                                        <  1:3451999/90‑1 (MQ=255)
                tCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAgt                                                                      <  1:559953/90‑1 (MQ=255)
                tCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAgt                                                                      <  1:1084287/90‑1 (MQ=255)
                tCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAgt                                                                      <  1:1509343/90‑1 (MQ=255)
                 cAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAgtg                                                                     <  1:3215645/90‑1 (MQ=255)
                  aGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAgtgt                                                                    <  1:2987511/90‑1 (MQ=255)
                  aGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAgtgt                                                                    >  1:1165034/1‑90 (MQ=255)
                  aGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAgtgt                                                                    <  1:1104103/90‑1 (MQ=255)
                     ttATCGATGTCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGtttt                                                                 <  2:1967200/90‑1 (MQ=255)
                       aTCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGtttttt                                                               <  1:2093839/90‑1 (MQ=255)
                       aTCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGtttttt                                                               <  1:201121/90‑1 (MQ=255)
                        tCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTg                                                              <  2:1442779/90‑1 (MQ=255)
                         cGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGa                                                             >  2:2446538/1‑90 (MQ=255)
                         cGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGa                                                             <  2:2755656/90‑1 (MQ=255)
                          gATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGat                                                            <  1:1009967/90‑1 (MQ=255)
                          gATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGat                                                            >  2:2721119/1‑90 (MQ=255)
                               ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGa                                                       >  2:715129/1‑90 (MQ=255)
                                cGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGaa                                                      <  1:604928/90‑1 (MQ=255)
                                cGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGaa                                                      >  2:2761862/1‑90 (MQ=255)
                                cGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGaa                                                      >  2:440235/1‑90 (MQ=255)
                                 gCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAg                                                     <  1:523878/90‑1 (MQ=255)
                                  cATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAgc                                                    >  2:2479440/1‑90 (MQ=255)
                                    tGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAgcgc                                                  >  1:1794834/1‑90 (MQ=255)
                                     ggCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCt                                                 <  1:2581202/90‑1 (MQ=255)
                                     ggCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCt                                                 >  1:1612792/1‑90 (MQ=255)
                                      gCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCtt                                                <  2:3010119/90‑1 (MQ=255)
                                       cGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCttt                                               <  2:3183051/90‑1 (MQ=255)
                                       cGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCttt                                               <  1:102415/90‑1 (MQ=255)
                                        ggcGCCTGGACAGGAGGATCGTAACGTTGCCCAGGAGGATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTc                                              >  2:3400811/3‑90 (MQ=255)
                                        gTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTc                                              <  1:1562162/90‑1 (MQ=255)
                                         tCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCt                                             >  1:1058579/1‑90 (MQ=255)
                                         tCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCt                                             <  1:402226/90‑1 (MQ=255)
                                         tCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCt                                             >  2:1545129/1‑90 (MQ=255)
                                          cGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTg                                            >  2:1366622/1‑90 (MQ=255)
                                           gCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGa                                           <  1:74527/90‑1 (MQ=255)
                                           gCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGa                                           <  2:1173990/90‑1 (MQ=255)
                                            ccTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGAt                                          <  1:1931444/90‑1 (MQ=255)
                                             cTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATc                                         >  1:343419/1‑90 (MQ=255)
                                              tGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATca                                        >  1:1839885/1‑90 (MQ=255)
                                               ggACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATcac                                       >  1:1289299/1‑90 (MQ=255)
                                                gACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATcaca                                      >  2:2943101/1‑90 (MQ=255)
                                                  cAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACt                                    <  1:2468435/90‑1 (MQ=255)
                                                  cAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACt                                    >  1:2829315/1‑90 (MQ=255)
                                                   aggaagATCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACtt                                   >  2:1940765/1‑90 (MQ=255)
                                                     gaggaTCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTTTTTTTTGATATCGAAGCGCTTTCTGATCACACTTcc                                 <  1:2676942/90‑1 (MQ=255)
                                                       ggaTCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTcccc                               <  2:3451211/90‑1 (MQ=255)
                                                       ggaTCGTAACGTTGCGCAAGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTcccc                               >  2:2736797/1‑90 (MQ=255)
                                                        gaTCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCg                              >  2:2643069/1‑90 (MQ=255)
                                                         aTCGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGc                             >  1:401514/1‑90 (MQ=255)
                                                           cGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACTCTTCCCCGCtt                           <  2:343341/90‑1 (MQ=255)
                                                           cGTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCtt                           >  2:1074247/1‑90 (MQ=255)
                                                            gTAACGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTc                          <  2:610987/90‑1 (MQ=255)
                                                              aaCGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCg                        <  2:3289664/90‑1 (MQ=255)
                                                              aaCGTTGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGAAATCGAAGCGCTTTCTGATCACCCTTCCCCTCTTCCg                        <  1:2678413/90‑1 (MQ=255)
                                                                  ttGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGcaca                    <  2:1162905/90‑1 (MQ=255)
                                                                   tGCTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACAt                   <  2:2701778/90‑1 (MQ=255)
                                                                     cTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACAtgc                 >  2:61327/1‑90 (MQ=255)
                                                                     cTCAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACAtgc                 >  1:1564519/1‑90 (MQ=255)
                                                                       cAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACAtgctg               <  2:1264879/90‑1 (MQ=255)
                                                                       cAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACAtgctg               >  1:721017/1‑90 (MQ=255)
                                                                       cAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACAtgctg               <  2:2906096/90‑1 (MQ=255)
                                                                       cAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACAtgctg               >  1:2196520/1‑90 (MQ=255)
                                                                       cAGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACAtgctg               <  1:3429808/90‑1 (MQ=255)
                                                                        aGGAGTATCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACAtgctgc              <  2:200564/90‑1 (MQ=255)
                                                                              aTCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGcc        >  1:1904426/1‑90 (MQ=255)
                                                                               tCTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCGGATCCCACTTCCCCGCTTCCGCACATGCTGCCGCGccc       <  2:454792/90‑1 (MQ=255)
                                                                                cTGAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCg      >  1:1194645/1‑90 (MQ=255)
                                                                                 tgaatgaaCATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCgg     <  1:877186/90‑1 (MQ=255)
                                                                                 tgaatgaaCATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCgg     >  2:2828333/1‑90 (MQ=255)
                                                                                  gAATGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGa    <  2:3007304/90‑1 (MQ=255)
                                                                                   aaTGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGaa   <  2:1377081/90‑1 (MQ=255)
                                                                                    aTGAACATATTGCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGaaa  >  1:1635874/1‑90 (MQ=255)
                                                                                    aTGAACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGaaa  >  2:1602320/1‑90 (MQ=255)
                                                                                       |                                                                                      
TTGATGACCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAA  >  AC_000091/2651196‑2651369

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