Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 31,345 A→G S177G (AGC→GGC)  carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091331,3450AG80.2% 489.8 / 134.5 253S177G (AGC→GGC) carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base A (23/27);  new base G (105/98);  total (128/125)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 260 total aligned reads displayed.

CCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTGCTGATTGAT  >  NC_000913/31253‑31437
                                                                                            |                                                                                            
ccGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTa                                                                                              <  1:2577464/93‑1 (MQ=255)
   gTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGgcg                                                                                           >  2:4506532/1‑93 (MQ=255)
    tATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGgcgg                                                                                          <  2:1583059/93‑1 (MQ=255)
    tATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGgcgg                                                                                          >  2:2400849/1‑93 (MQ=255)
      tCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGgcggcg                                                                                        >  2:4071556/1‑93 (MQ=255)
       cGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGgcggcgg                                                                                       <  1:2739441/93‑1 (MQ=255)
        gCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTggcggcggc                                                                                      <  2:2617850/93‑1 (MQ=255)
        gCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTggcggcggc                                                                                      >  1:3531535/1‑93 (MQ=255)
           cacacGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGt                                                                                   <  1:3323703/93‑1 (MQ=255)
            acacGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTa                                                                                  >  1:1522529/1‑93 (MQ=255)
            acacGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTa                                                                                  <  2:605053/93‑1 (MQ=255)
            acacGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTa                                                                                  <  2:6789773/93‑1 (MQ=255)
              acGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATc                                                                                <  1:2025826/93‑1 (MQ=255)
              acGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATc                                                                                <  2:5522465/93‑1 (MQ=255)
              acGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATc                                                                                <  2:100090/93‑1 (MQ=255)
              acGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATc                                                                                <  1:6897513/93‑1 (MQ=255)
                gATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGc                                                                              <  1:2805894/93‑1 (MQ=255)
                gATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGc                                                                              <  1:2970156/93‑1 (MQ=255)
                gATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGc                                                                              <  2:5913273/93‑1 (MQ=255)
                 aTGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCt                                                                             >  2:3983157/1‑93 (MQ=255)
                   ggAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTta                                                                           <  2:4133075/93‑1 (MQ=255)
                   ggAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTta                                                                           <  2:1786192/93‑1 (MQ=255)
                    gaagaaGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTtat                                                                          >  2:5118729/1‑93 (MQ=255)
                      agaagCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATaa                                                                        >  2:7164506/1‑93 (MQ=255)
                      agaagCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATaa                                                                        <  2:3664772/93‑1 (MQ=255)
                       gaagCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAAc                                                                       <  1:4993674/93‑1 (MQ=255)
                       gaagCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAAc                                                                       <  2:5117379/93‑1 (MQ=255)
                        aagCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAAcc                                                                      <  1:6336141/93‑1 (MQ=255)
                            gcTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGa                                                                  <  1:2558024/93‑1 (MQ=255)
                            gcTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGa                                                                  >  1:3950238/1‑93 (MQ=255)
                            gcTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGa                                                                  >  2:2551312/1‑93 (MQ=255)
                             cTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGaa                                                                 >  1:4090657/1‑93 (MQ=255)
                              tGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAag                                                                >  1:1338477/1‑93 (MQ=255)
                                  ggTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGtt                                                            >  1:467117/1‑93 (MQ=255)
                                  ggTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGtt                                                            >  1:7465/1‑93 (MQ=255)
                                        cGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTgaaga                                                      >  1:1641867/1‑93 (MQ=255)
                                        cGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTgaaga                                                      <  1:6717144/93‑1 (MQ=255)
                                         gCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTgaagaa                                                     >  2:2854422/1‑93 (MQ=255)
                                          cTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGaaa                                                    <  2:1502879/93‑1 (MQ=255)
                                          cTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGaaa                                                    >  2:7240632/1‑93 (MQ=255)
                                           tGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAAt                                                   >  2:7501180/1‑93 (MQ=255)
                                           tGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAAt                                                   >  2:1356126/1‑93 (MQ=255)
                                           tGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAAt                                                   <  2:2865821/93‑1 (MQ=255)
                                            gACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAAtt                                                  >  1:5548683/1‑93 (MQ=255)
                                             aCGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAAttt                                                 <  1:6026344/93‑1 (MQ=255)
                                                tGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTgcg                                              >  2:6687298/1‑93 (MQ=255)
                                                 gggCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTgcgc                                             >  1:4422379/1‑93 (MQ=255)
                                                 gggCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTgcgc                                             >  2:3850277/1‑93 (MQ=255)
                                                  ggCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGcc                                            >  2:1259341/1‑93 (MQ=255)
                                                     ttCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCcgc                                         <  2:1672080/93‑1 (MQ=255)
                                                      tCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCcgcg                                        >  1:5781994/1‑93 (MQ=255)
                                                         cGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAATTTTGCGCCCGCGGTc                                     <  2:7467001/93‑1 (MQ=255)
                                                         cGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTc                                     <  1:1761737/93‑1 (MQ=255)
                                                         cGTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTc                                     <  1:1157017/93‑1 (MQ=255)
                                                          gTGCATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCt                                    >  2:861910/1‑93 (MQ=255)
                                                             cATTATTCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGa                                 <  2:2025340/93‑1 (MQ=255)
                                                              attattCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAt                                >  2:7025170/1‑93 (MQ=255)
                                                               ttattCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAtc                               <  1:3720803/93‑1 (MQ=255)
                                                                tattCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAtct                              >  2:5295807/1‑93 (MQ=255)
                                                                tattCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAtct                              >  1:5151634/1‑93 (MQ=255)
                                                                 attCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAtctc                             <  1:1023888/93‑1 (MQ=255)
                                                                 attCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAtctc                             <  2:5721796/93‑1 (MQ=255)
                                                                  ttCGCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAtctct                            >  1:3663442/1‑93 (MQ=255)
                                                                     gCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAtctctctc                         <  1:6771957/93‑1 (MQ=255)
                                                                     gCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAtctctctc                         <  2:3267008/93‑1 (MQ=255)
                                                                     gCCCATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAtctctctc                         <  1:6619308/93‑1 (MQ=255)
                                                                     gCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGAtctctctc                         <  1:7006206/93‑1 (MQ=255)
                                                                      cccATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTcc                        >  1:742941/1‑93 (MQ=255)
                                                                      cccATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTcc                        <  1:2438672/93‑1 (MQ=255)
                                                                      cccATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTcc                        <  2:1912958/93‑1 (MQ=255)
                                                                      cccATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTcc                        <  2:6203570/93‑1 (MQ=255)
                                                                        cATCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGa                      >  2:639023/1‑93 (MQ=255)
                                                                        cATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGa                      >  1:7363372/1‑93 (MQ=255)
                                                                         aTCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGAc                     >  2:3469948/1‑93 (MQ=255)
                                                                         aTCCTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGAc                     >  1:3965668/1‑93 (MQ=255)
                                                                           ccTTTACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCa                   >  2:4776063/1‑93 (MQ=255)
                                                                            cTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCaa                  <  2:4204920/93‑1 (MQ=255)
                                                                             tttACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCaaa                 >  2:3909179/1‑93 (MQ=255)
                                                                             tttACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCaaa                 >  2:2965914/1‑93 (MQ=255)
                                                                             tttACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCaaa                 >  2:2228093/1‑93 (MQ=255)
                                                                              ttACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAag                >  2:3520796/1‑93 (MQ=255)
                                                                               tACCATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAaga               <  1:5348382/93‑1 (MQ=255)
                                                                               tACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAaga               <  1:6435909/93‑1 (MQ=255)
                                                                               tACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAaga               <  1:7251673/93‑1 (MQ=255)
                                                                                 ccATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGt             >  1:4205495/1‑93 (MQ=255)
                                                                                 ccATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGt             >  2:6477454/1‑93 (MQ=255)
                                                                                 ccATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGt             >  1:3582074/1‑93 (MQ=255)
                                                                                 ccATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGt             >  2:1088075/1‑93 (MQ=255)
                                                                                  cATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGtt            >  2:1176592/1‑93 (MQ=255)
                                                                                  cATGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGtt            <  1:1262167/93‑1 (MQ=255)
                                                                                  cATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGtt            >  2:558211/1‑93 (MQ=255)
                                                                                   aTGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTg           >  1:4255523/1‑93 (MQ=255)
                                                                                    tGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTGc          >  1:5507466/1‑93 (MQ=255)
                                                                                    tGGGCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTGc          >  1:158101/1‑93 (MQ=255)
                                                                                     gggCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTGCt         <  1:2903601/93‑1 (MQ=255)
                                                                                       gCGGTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTGCtga       <  2:4683074/93‑1 (MQ=255)
                                                                                       gCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTGCtga       >  2:544616/1‑93 (MQ=255)
                                                                                         ggTGGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTGCtgatt     <  1:5883098/93‑1 (MQ=255)
                                                                                         ggTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTGCtgatt     <  2:2664496/93‑1 (MQ=255)
                                                                                            aGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTGCtgattgat  <  1:2874055/93‑1 (MQ=255)
                                                                                            |                                                                                            
CCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGTGGGCTTCCCGTGCATTATTCGCCCATCCTTTACCATGGGCGGTAGCGGCGGCGGTATCGCTTATAACCGTGAAGAGTTTGAAGAAATTTGCGCCCGCGGTCTGGATCTCTCTCCGACCAAAGAGTTGCTGATTGAT  >  NC_000913/31253‑31437

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