Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,847,491 T→C L91L (TTA→TTG)  b3668 ← sensor histidine protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,847,4910TC99.8% 1376.3 / ‑5.2 417L91L (TTA→TTG) b3668sensor histidine protein kinase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (206/210);  total (206/211)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 430 total aligned reads displayed.

CGCCGTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACG  >  NC_000913/3847401‑3847581
                                                                                          |                                                                                          
 gccgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCaaa                                                                                         >  2:5841306/1‑93 (MQ=255)
 gccgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCaaa                                                                                         <  1:2155883/93‑1 (MQ=255)
  ccgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAAc                                                                                        <  1:7562224/93‑1 (MQ=255)
  ccgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAAc                                                                                        <  1:1729085/93‑1 (MQ=255)
  ccgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAAc                                                                                        >  2:8318343/1‑93 (MQ=255)
   cgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACt                                                                                       <  1:8340417/93‑1 (MQ=255)
    gTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTa                                                                                      <  2:7445066/93‑1 (MQ=255)
     ttAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTAc                                                                                     >  2:9817289/1‑93 (MQ=255)
     ttAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTAc                                                                                     <  2:10090/93‑1 (MQ=255)
     ttAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTAc                                                                                     >  2:6896072/1‑93 (MQ=255)
     ttAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTAc                                                                                     >  1:4911698/1‑93 (MQ=255)
     ttAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTAc                                                                                     >  1:3903055/1‑93 (MQ=255)
     ttAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTAc                                                                                     >  1:3770591/1‑93 (MQ=255)
           ccgccCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGatca                                                                               <  2:2097864/93‑1 (MQ=255)
             gccCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGatcatc                                                                             >  2:9402682/1‑93 (MQ=255)
               ccccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCaa                                                                           >  2:2996504/1‑93 (MQ=255)
                cccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAAt                                                                          >  2:471349/1‑93 (MQ=255)
                cccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAAt                                                                          >  1:7740162/1‑93 (MQ=255)
                cccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAAt                                                                          >  1:588914/1‑93 (MQ=255)
                cccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAAt                                                                          >  1:4336865/1‑93 (MQ=255)
                 ccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATa                                                                         <  1:1386635/93‑1 (MQ=255)
                 ccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATa                                                                         >  2:2472777/1‑93 (MQ=255)
                 catGTAACAGCAAGGGGCGCCAGTCACGCTGATGGCGATAGCCCGAGACCAGCGTTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATa                                                                         <  2:7688878/91‑1 (MQ=255)
                  cTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCCGTAACGTCAGAAAACTACCGATCATCAATac                                                                        >  1:2753313/1‑92 (MQ=255)
                  cTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATaa                                                                        >  1:194045/1‑93 (MQ=255)
                  cTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCAATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATaa                                                                        <  2:4177590/93‑1 (MQ=255)
                   tGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAAc                                                                       >  1:1607451/1‑93 (MQ=255)
                   tGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAAc                                                                       <  1:1115016/93‑1 (MQ=255)
                    gTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACg                                                                      <  1:2356670/93‑1 (MQ=255)
                    gTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACg                                                                      >  2:6761898/1‑93 (MQ=255)
                      aaCAGCAAGGTGCGCCAGTCACGCTGAGGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGa                                                                    <  2:1873681/93‑1 (MQ=255)
                       aCAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGaa                                                                   <  2:2332444/93‑1 (MQ=255)
                       aCAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGaa                                                                   >  2:1500720/1‑93 (MQ=255)
                       aCAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGaa                                                                   >  1:3678314/1‑93 (MQ=255)
                        cAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGaaa                                                                  <  2:6544024/93‑1 (MQ=255)
                        cAGCAAGGTGCGCCAGACACGCTGATGGCGATAGCGCGAGAGCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGaaa                                                                  <  2:4431620/93‑1 (MQ=255)
                         aGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGaaaa                                                                 >  2:3728001/1‑93 (MQ=255)
                          gCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAAt                                                                <  2:779766/93‑1 (MQ=255)
                          gCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAAt                                                                <  2:4895517/93‑1 (MQ=255)
                           cAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATg                                                               <  1:308988/93‑1 (MQ=255)
                            aaGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATgg                                                              <  1:8188851/93‑1 (MQ=255)
                             aGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATggg                                                             >  1:9956396/1‑93 (MQ=255)
                              ggTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGt                                                            >  1:2559675/1‑93 (MQ=255)
                              ggTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGt                                                            <  1:8988393/93‑1 (MQ=255)
                               gTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGtt                                                           >  2:3551351/1‑93 (MQ=255)
                               gTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGta                                                           >  1:1842670/1‑92 (MQ=255)
                                tGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGGTa                                                          >  1:6245149/1‑93 (MQ=255)
                                 tctccAGGCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTaa                                                         <  2:2510442/90‑1 (MQ=255)
                                  cgcCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTaaa                                                        <  1:8640219/93‑1 (MQ=255)
                                  cgcCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTaaa                                                        <  2:7293739/93‑1 (MQ=255)
                                   gcCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAac                                                       <  1:4760877/93‑1 (MQ=255)
                                    ccAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAacg                                                      >  1:6662289/1‑92 (MQ=255)
                                    ccAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAacc                                                      >  1:5515466/1‑93 (MQ=255)
                                     cAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccg                                                     >  2:3439407/1‑93 (MQ=255)
                                      aGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccga                                                    <  2:612394/93‑1 (MQ=255)
                                      aGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccga                                                    <  2:9117023/93‑1 (MQ=255)
                                      aGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccga                                                    <  2:8219240/93‑1 (MQ=255)
                                       gTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccgac                                                   <  2:8466947/93‑1 (MQ=255)
                                        tCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccgacc                                                  >  1:9605428/1‑93 (MQ=255)
                                        tCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccgacc                                                  <  2:5596120/93‑1 (MQ=255)
                                        tCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccgacc                                                  <  1:9684115/93‑1 (MQ=255)
                                        tCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccgacc                                                  <  1:6376068/93‑1 (MQ=255)
                                        tCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccgacc                                                  >  1:339124/1‑93 (MQ=255)
                                          aCGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGc                                                >  1:532339/1‑93 (MQ=255)
                                           cGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGcc                                               <  2:3877919/93‑1 (MQ=255)
                                             cTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTg                                             >  2:1354936/1‑93 (MQ=255)
                                              tGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGc                                            >  1:8039732/1‑93 (MQ=255)
                                                aTGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGt                                          <  2:5977617/93‑1 (MQ=255)
                                                aTGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGt                                          >  1:8792959/1‑93 (MQ=255)
                                                 tGGCGATAGCGCGTGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGtt                                         <  1:9350240/93‑1 (MQ=255)
                                                 tGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGtt                                         >  2:1683427/1‑93 (MQ=255)
                                                      aTAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCa                                    >  2:5222880/1‑93 (MQ=255)
                                                      aTAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCa                                    >  1:3250455/1‑93 (MQ=255)
                                                         gcgcGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTa                                 >  1:4431965/1‑93 (MQ=255)
                                                           gcgAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTaaa                               >  2:2001055/1‑93 (MQ=255)
                                                            cgAGATCAGGGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAAt                              <  2:6496919/93‑1 (MQ=255)
                                                             gagaTCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATc                             <  1:9288451/93‑1 (MQ=255)
                                                             gagaTCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATc                             <  2:8057924/93‑1 (MQ=255)
                                                               gaTCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATcag                           >  1:1276245/1‑93 (MQ=255)
                                                                   aGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCagcagc                       >  1:2848928/1‑93 (MQ=255)
                                                                    gcgcTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGcc                      <  2:5791695/93‑1 (MQ=255)
                                                                     cgcTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCa                     <  1:7900982/93‑1 (MQ=255)
                                                                     cgcTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCa                     <  1:7399278/93‑1 (MQ=255)
                                                                      gcTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCAc                    <  1:10032087/93‑1 (MQ=255)
                                                                        tACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTc                  <  1:8216904/93‑1 (MQ=255)
                                                                        tACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTc                  <  2:8724599/93‑1 (MQ=255)
                                                                         aCGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTcc                 >  1:1031139/1‑93 (MQ=255)
                                                                         aCGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTcc                 <  2:3173341/93‑1 (MQ=255)
                                                                           gggCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCcgc               <  1:9455322/93‑1 (MQ=255)
                                                                            ggCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCcgcg              >  2:8971389/1‑93 (MQ=255)
                                                                            ggCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCcgcg              <  2:2746464/93‑1 (MQ=255)
                                                                            ggCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCcgcg              >  2:1799035/1‑93 (MQ=255)
                                                                               aGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCcccc           >  1:3589042/1‑93 (MQ=255)
                                                                               aGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGATAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGccc           >  1:7238468/1‑93 (MQ=255)
                                                                                gTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCa          <  1:6134522/93‑1 (MQ=255)
                                                                                 tAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAg         >  1:692795/1‑93 (MQ=255)
                                                                                   aCGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCa       >  1:4053187/1‑93 (MQ=255)
                                                                                    cGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCaa      >  1:4528641/1‑93 (MQ=255)
                                                                                     gTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAAt     <  2:5009982/93‑1 (MQ=255)
                                                                                      tCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATa    >  2:9679175/1‑93 (MQ=255)
                                                                                          |                                                                                          
CGCCGTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACG  >  NC_000913/3847401‑3847581

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