Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 31,762 G→A G316S (GGT→AGT)  carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091331,7620GA99.5% 658.9 / ‑5.8 184G316S (GGT→AGT) carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (92/90);  total (93/91)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 187 total aligned reads displayed.

TTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCACTGGCGGACGTACTCC  >  NC_000913/31672‑31850
                                                                                          |                                                                                        
tttGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGCCCAGt                                                                                        >  2:1408730/1‑93 (MQ=255)
 ttGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGtt                                                                                       <  1:2150705/93‑1 (MQ=255)
   gCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTc                                                                                     <  1:1707482/93‑1 (MQ=255)
     ggTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTccc                                                                                   >  1:503543/1‑93 (MQ=255)
     ggTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTccc                                                                                   <  2:2190461/93‑1 (MQ=255)
       tGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGa                                                                                 <  2:767838/93‑1 (MQ=255)
       tGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGa                                                                                 <  2:140633/93‑1 (MQ=255)
       tGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGa                                                                                 <  1:630647/93‑1 (MQ=255)
       tGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGa                                                                                 <  2:277340/93‑1 (MQ=255)
        gAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGAt                                                                                >  1:2542265/1‑93 (MQ=255)
        gAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGAt                                                                                <  1:4627704/93‑1 (MQ=255)
         aaCCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGAtt                                                                               <  1:1061502/93‑1 (MQ=255)
            ccGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCt                                                                            <  1:2110563/93‑1 (MQ=255)
            ccGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCt                                                                            >  1:3018331/1‑93 (MQ=255)
             cGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTa                                                                           <  1:3095875/93‑1 (MQ=255)
               aaaaaCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTaaa                                                                         <  1:4579232/93‑1 (MQ=255)
               aaaaaCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTaaa                                                                         <  2:1282722/93‑1 (MQ=255)
                 aaaCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGt                                                                       <  2:2426300/93‑1 (MQ=255)
                 aaaCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGt                                                                       <  2:1590462/93‑1 (MQ=255)
                 aaaCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGt                                                                       <  1:2696021/93‑1 (MQ=255)
                 aaaCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAACGACCAGTTTCCCGATTCCTAAAgg                                                                       >  1:3247549/1‑92 (MQ=255)
                  aaCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTg                                                                      <  1:1591546/93‑1 (MQ=255)
                  aaCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTg                                                                      <  1:4422347/93‑1 (MQ=255)
                  aaCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTg                                                                      >  2:496038/1‑93 (MQ=255)
                   aCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTgg                                                                     <  2:1838926/93‑1 (MQ=255)
                   aCGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTgg                                                                     <  2:1628714/93‑1 (MQ=255)
                    cGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTggg                                                                    >  1:3143307/1‑92 (MQ=255)
                    cGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTggc                                                                    <  2:2107920/93‑1 (MQ=255)
                     ggTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTggcg                                                                   <  2:1957564/93‑1 (MQ=255)
                     ggTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTggcg                                                                   >  2:2883968/1‑93 (MQ=255)
                     ggTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTggcg                                                                   <  2:4135619/93‑1 (MQ=255)
                     ggTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTggcg                                                                   >  1:191654/1‑93 (MQ=255)
                      gtcgtcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTggcgg                                                                  >  2:1990050/1‑93 (MQ=255)
                      gtcgtcTCATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTggcgg                                                                  >  2:4662169/1‑93 (MQ=255)
                        cgtcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGgcggcg                                                                <  2:2078451/93‑1 (MQ=255)
                         gtcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGa                                                               >  2:1009791/1‑93 (MQ=255)
                          tcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGgcggcgca                                                              >  2:3778622/1‑91 (MQ=255)
                          tcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGgcggcgca                                                              >  1:1740476/1‑91 (MQ=255)
                          tcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGaa                                                              >  2:399328/1‑93 (MQ=255)
                           cTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGaaa                                                             >  2:842757/1‑93 (MQ=255)
                            tGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAAc                                                            <  2:1776377/93‑1 (MQ=255)
                            tGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAAc                                                            <  2:4563815/93‑1 (MQ=255)
                            tGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAAc                                                            <  2:3923077/93‑1 (MQ=255)
                             gATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACt                                                           >  1:3733446/1‑93 (MQ=255)
                             gATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACt                                                           <  2:3487668/93‑1 (MQ=255)
                              aTTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTg                                                          >  1:1829532/1‑93 (MQ=255)
                               ttGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTgg                                                         >  1:3690444/1‑93 (MQ=255)
                                tGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGc                                                        <  2:3848523/93‑1 (MQ=255)
                                    aTCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTg                                                    >  2:2350454/1‑93 (MQ=255)
                                     tCGAAATGACCCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTgg                                                   <  1:4133681/93‑1 (MQ=255)
                                     tCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTgg                                                   <  1:579452/93‑1 (MQ=255)
                                       gAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAAATGGCGGTGGGt                                                 >  2:1388102/1‑93 (MQ=255)
                                           tGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTaca                                             <  2:2255095/93‑1 (MQ=255)
                                           tGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTaca                                             <  1:2748064/93‑1 (MQ=255)
                                            gAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTacac                                            <  1:4179793/93‑1 (MQ=255)
                                            gAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGGTacac                                            >  2:1128277/1‑93 (MQ=255)
                                                ccACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTc                                        >  2:4293446/1‑93 (MQ=255)
                                                ccACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTc                                        >  1:3405351/1‑93 (MQ=255)
                                                 cACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTcg                                       >  1:3149744/1‑93 (MQ=255)
                                                 cACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGGTACACCCTcg                                       >  1:810837/1‑93 (MQ=255)
                                                  aCGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTcga                                      <  2:1117876/93‑1 (MQ=255)
                                                   cgcgTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTcgac                                     >  2:3636711/1‑93 (MQ=255)
                                                   cgcgTGTCCCGTTCTTCCGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTcgac                                     <  2:874693/93‑1 (MQ=255)
                                                    gcgTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTcgacg                                    >  1:712716/1‑93 (MQ=255)
                                                    gcgTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGGGGGTTACACCCTcgacg                                    >  1:730529/1‑93 (MQ=255)
                                                    gcgTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGGGGGCTAGACCCTccacc                                    >  1:344257/1‑92 (MQ=255)
                                                       tgtCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAAc                                 <  1:1350567/93‑1 (MQ=255)
                                                        gtCCCGTTCTTCGGCGCTGGCGTCGAAAGAGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGGTACACCCTCGACGCACt                                >  2:1534751/1‑93 (MQ=255)
                                                          cccGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACtga                              <  1:451328/93‑1 (MQ=255)
                                                          cccGTTCGTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACtga                              <  1:351554/93‑1 (MQ=255)
                                                           ccGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACtgat                             >  1:3890822/1‑93 (MQ=255)
                                                           ccGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACtgat                             >  1:1599373/1‑93 (MQ=255)
                                                           ccGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGGTACACCCTCGACGAACtgat                             >  1:3204259/1‑93 (MQ=255)
                                                            cGTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACtgatg                            <  1:2630155/93‑1 (MQ=255)
                                                             gTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACtgatga                           <  1:1302823/93‑1 (MQ=255)
                                                             gTTCTTCGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACtgatga                           >  1:4160758/1‑93 (MQ=255)
                                                               tcttcGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAAc                         <  1:124801/93‑1 (MQ=255)
                                                               tcttcGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAAc                         <  2:3771248/93‑1 (MQ=255)
                                                               tcttcGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAAc                         <  1:953889/93‑1 (MQ=255)
                                                               tcttcGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGGTACACCCTCGACGAACTGATGAAc                         >  2:253075/1‑93 (MQ=255)
                                                                cttcGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACg                        <  1:2378449/93‑1 (MQ=255)
                                                                 ttcGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGa                       <  1:3337060/93‑1 (MQ=255)
                                                                   cGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACa                     <  2:3094437/93‑1 (MQ=255)
                                                                   cGGCGCTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACa                     <  2:1702612/93‑1 (MQ=255)
                                                                       gcTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCAc                 <  2:4215319/93‑1 (MQ=255)
                                                                       gcTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATCAACCACATcgc                 >  1:640020/1‑91 (MQ=255)
                                                                       gcTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACCAACTGATGAACGACATCAc                 >  1:1462809/1‑93 (MQ=255)
                                                                       gcTGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGGTACACCCTCGACGAACTGATGAACGACATCAc                 >  2:2405763/1‑93 (MQ=255)
                                                                         tGGCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCACTg               <  1:1047040/93‑1 (MQ=255)
                                                                         cggCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAATTGGGGGGGGGTTAAACCCCCGCCGAACTGATGAACAAAATCCCTg               >  2:2168716/2‑93 (MQ=255)
                                                                           gCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGACCTGATGAACGACATCACTGGc             <  1:2046464/93‑1 (MQ=255)
                                                                           gCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCACTGGc             >  1:3540354/1‑93 (MQ=255)
                                                                           gCGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAAATGATGAACGACATACCTGGc             >  2:2339931/1‑93 (MQ=255)
                                                                            cGTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCACTGGCg            <  1:379354/93‑1 (MQ=255)
                                                                             gTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCACTGgggg           >  1:1144781/1‑93 (MQ=255)
                                                                             gTCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCACTGGCgg           >  1:3463160/1‑93 (MQ=255)
                                                                              tCGAAAGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGATGAACTGATGAACGACATCACTGGCggg          >  2:705086/1‑92 (MQ=255)
                                                                                 aaaGCGACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCACTGGCGGACgg       >  2:1577017/1‑92 (MQ=255)
                                                                                  aaGCGACCCGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCACTGGCGGGCGTa      >  2:4561050/1‑93 (MQ=255)
                                                                                      gACCAGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCTCTGGCGGACGTACTcc  <  2:1761397/93‑1 (MQ=255)
                                                                                          |                                                                                        
TTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCGATTGCTAAAGTGGCGGCGAAACTGGCGGTGGGTTACACCCTCGACGAACTGATGAACGACATCACTGGCGGACGTACTCC  >  NC_000913/31672‑31850

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