Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 706,394 (T)9→8 intergenic (+82/‑121) nagE → / → glnS fused N‑acetyl glucosamine specific PTS enzyme IICBA components/glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_000091706,3860T.99.2% 603.1 / ‑2.0 117intergenic (+74/‑129)nagE/glnSfused N‑acetyl glucosamine specific PTS enzyme IICBA components/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (58/58);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 141 total aligned reads displayed.

TGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAG  >  AC_000091/706297‑706471
                                                                                         |                                                                                     
 gAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtt                                                                                      >  2:1795040/1‑90 (MQ=255)
 gAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtt                                                                                      >  1:1284163/1‑90 (MQ=255)
 gAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtt                                                                                      <  2:1619443/90‑1 (MQ=255)
 gAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtt                                                                                      <  1:3025377/90‑1 (MQ=255)
  aaaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCttt                                                                                     <  1:139122/90‑1 (MQ=255)
   aaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtttt                                                                                    <  2:333033/90‑1 (MQ=255)
    aTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCttttt                                                                                   <  2:1799564/90‑1 (MQ=255)
      cAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCttttttt                                                                                 <  2:1450906/90‑1 (MQ=255)
      cAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCttttttt                                                                                 <  2:862067/90‑1 (MQ=255)
       aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtttttttt                                                                                >  2:2852848/1‑90 (MQ=255)
       aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTCGGGCGTCGCATCCAACAATGACAAGCGGTGGAGAACTTCTCTCCCGCtttttttt                                                                                >  2:713071/1‑90 (MQ=255)
       aaaaaGCAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtttttttt                                                                                >  2:1537427/1‑90 (MQ=255)
         aaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCttttttttca                                                                              >  2:1619469/1‑88 (MQ=255)
           aGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCttttttttcatc                                                                            >  1:3334384/1‑86 (MQ=255)
              aaTCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCatcaat                                                                        >  1:1017289/1‑90 (MQ=255)
              aaTCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCatcaat                                                                        >  2:2669532/1‑90 (MQ=255)
               aTCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCatcaatc                                                                       >  1:1607675/1‑90 (MQ=255)
               aTCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCatcaatc                                                                       <  2:163048/90‑1 (MQ=255)
               aTCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCatcaatc                                                                       >  1:2146449/1‑90 (MQ=255)
                 cTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAatcat                                                                     >  1:1981081/1‑90 (MQ=255)
                   gCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATcc                                                                   <  1:3052639/90‑1 (MQ=255)
                   gCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATcc                                                                   >  2:2749864/1‑90 (MQ=255)
                   gCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATcc                                                                   >  1:1475112/1‑90 (MQ=255)
                    cTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATccc                                                                  >  2:3295998/1‑90 (MQ=255)
                     tttATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCTGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATcccc                                                                 <  2:1332248/90‑1 (MQ=255)
                     tttATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATcccc                                                                 <  1:2971738/90‑1 (MQ=255)
                     tttATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATcccc                                                                 >  2:2988993/1‑90 (MQ=255)
                       tATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCAt                                                               <  1:3219492/90‑1 (MQ=255)
                       tATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCAt                                                               >  2:1686466/1‑90 (MQ=255)
                        aTGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCTTAATTATTCCCAtt                                                              >  1:3405266/1‑89 (MQ=255)
                        aTGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATa                                                              >  1:3136183/1‑90 (MQ=255)
                         tGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATaa                                                             >  1:471745/1‑90 (MQ=255)
                          gCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACACGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAAt                                                            <  1:180909/90‑1 (MQ=255)
                          gCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAAt                                                            <  1:35509/90‑1 (MQ=255)
                           ccTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATc                                                           <  2:1663362/90‑1 (MQ=255)
                           ccTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATc                                                           <  1:1391058/90‑1 (MQ=255)
                           ccTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATc                                                           >  2:1928312/1‑90 (MQ=255)
                           ccTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATc                                                           >  1:1851629/1‑90 (MQ=255)
                              gATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCtt                                                        <  1:1247932/90‑1 (MQ=255)
                              gATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCtt                                                        >  2:2596543/1‑90 (MQ=255)
                               aTGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTg                                                       <  1:224190/90‑1 (MQ=255)
                                tGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGt                                                      <  1:329587/90‑1 (MQ=255)
                                 gCGTCGCATGAGCGTCGCATCCAACAATGACGAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGtt                                                     <  2:1266704/90‑1 (MQ=255)
                                  cGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTa                                                    >  2:3288711/1‑90 (MQ=255)
                                  cGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTa                                                    <  2:1004095/90‑1 (MQ=255)
                                    aCGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGa                                                  <  2:1506055/90‑1 (MQ=255)
                                     cGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAt                                                 <  2:3198849/90‑1 (MQ=255)
                                      gCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAtt                                                <  2:452163/90‑1 (MQ=255)
                                      gCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAtt                                                >  1:1866311/1‑90 (MQ=255)
                                      gCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAtt                                                <  2:1441336/90‑1 (MQ=255)
                                      gCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAtt                                                >  1:1471829/1‑90 (MQ=255)
                                        ttGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTAt                                              >  1:84438/1‑90 (MQ=255)
                                        ttGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTAt                                              >  2:2454544/1‑90 (MQ=255)
                                         tGAGCGGCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATc                                             >  2:3019323/1‑90 (MQ=255)
                                           aGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCaa                                           <  1:1035792/90‑1 (MQ=255)
                                           aGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCaa                                           >  1:2456866/1‑90 (MQ=255)
                                            gCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAt                                          <  1:2318662/90‑1 (MQ=255)
                                            gCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAt                                          >  2:3059381/1‑90 (MQ=255)
                                            gCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAt                                          >  1:1801030/1‑90 (MQ=255)
                                             cGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAtt                                         <  2:640940/90‑1 (MQ=255)
                                                cGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTa                                      <  2:1711956/90‑1 (MQ=255)
                                                cGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTa                                      <  2:2057650/90‑1 (MQ=255)
                                                 gCAGCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaa                                     <  2:2223589/90‑1 (MQ=255)
                                                   aTCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaa                                   <  1:1705171/90‑1 (MQ=255)
                                                   aTCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaa                                   >  2:2758206/1‑90 (MQ=255)
                                                   aTCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaa                                   <  1:2588576/90‑1 (MQ=255)
                                                    tCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaaa                                  >  1:2143597/1‑90 (MQ=255)
                                                    tCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaaa                                  >  2:836560/1‑90 (MQ=255)
                                                     ccAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaaaa                                 >  2:1963494/1‑90 (MQ=255)
                                                      caacaaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAAc                                >  1:2783902/1‑90 (MQ=255)
                                                        acaaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTa                              <  1:2821161/90‑1 (MQ=255)
                                                          aaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAAc                            >  2:1537471/1‑90 (MQ=255)
                                                           aTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACa                           >  1:2026486/1‑90 (MQ=255)
                                                            tGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAg                          <  1:2795233/90‑1 (MQ=255)
                                                            tGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAg                          <  2:3285680/90‑1 (MQ=255)
                                                            tGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAg                          >  2:2312348/1‑90 (MQ=255)
                                                              aCAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGtt                        <  2:2438636/90‑1 (MQ=255)
                                                              aCAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGtt                        <  2:1794888/90‑1 (MQ=255)
                                                              aCAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGtt                        >  2:577260/1‑90 (MQ=255)
                                                               cAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTg                       >  2:3134063/1‑90 (MQ=255)
                                                                aaGCGGTTGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGt                      >  1:1672977/1‑90 (MQ=255)
                                                                 aGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTc                     >  2:1106941/1‑90 (MQ=255)
                                                                 aGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTc                     <  1:2803417/90‑1 (MQ=255)
                                                                 aGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTc                     <  1:1095705/90‑1 (MQ=255)
                                                                   cGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAg                   <  2:567785/90‑1 (MQ=255)
                                                                   cGGTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAg                   >  2:2087012/1‑90 (MQ=255)
                                                                    ggTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGc                  >  2:78425/1‑90 (MQ=255)
                                                                     gTGGAGATCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGcc                 <  2:240011/90‑1 (MQ=255)
                                                                        gagaTCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGt              >  1:1878330/1‑90 (MQ=255)
                                                                         agaTCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTc             <  1:1148167/90‑1 (MQ=255)
                                                                         agaTCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTc             >  2:3443853/1‑90 (MQ=255)
                                                                           aTCTTCTCTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTccc           <  1:2377933/90‑1 (MQ=255)
                                                                               tctcTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCtt       >  2:1119145/1‑90 (MQ=255)
                                                                               tctcTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCtt       <  1:2516318/90‑1 (MQ=255)
                                                                               tctcTGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCtt       <  2:2337573/90‑1 (MQ=255)
                                                                                ctctGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTta      <  2:1124686/90‑1 (MQ=255)
                                                                                  ctGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTtata    >  2:3159819/1‑90 (MQ=255)
                                                                                  ctGCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTtata    <  2:1750233/90‑1 (MQ=255)
                                                                                    gCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAg  <  2:3330366/90‑1 (MQ=255)
                                                                                    gCCGCTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTCGTCAGCCTGTCCCGCTTATAAg  <  2:1916774/90‑1 (MQ=255)
                                                                                         |                                                                                     
TGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAG  >  AC_000091/706297‑706471

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