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.5% 1268.1 / ‑3.6 399L91L (TTA→TTG) b3668sensor histidine protein kinase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (209/188);  total (210/189)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

Only 100 of 410 total aligned reads displayed.

CGCCGTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATAC  >  NC_000913/3847401‑3847580
                                                                                          |                                                                                         
cgccgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCaa                                                                                         <  2:1642730/93‑1 (MQ=255)
ccccgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCaa                                                                                         <  2:3589507/91‑1 (MQ=255)
 gccgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCaaa                                                                                        >  2:2255612/1‑93 (MQ=255)
 gccgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCaaa                                                                                        >  2:667877/1‑93 (MQ=255)
  ccgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAAc                                                                                       <  1:2457046/93‑1 (MQ=255)
   cgTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACt                                                                                      >  2:1429184/1‑93 (MQ=255)
    gTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTa                                                                                     >  2:1216269/1‑93 (MQ=255)
    gTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTa                                                                                     >  2:5705991/1‑93 (MQ=255)
     ttAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTAc                                                                                    >  2:3221459/1‑93 (MQ=255)
     ttAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTAc                                                                                    >  1:2275643/1‑93 (MQ=255)
     gtAACGCCTCCCCCTGTAACAGCAAGGTGCCCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTAc                                                                                    <  1:1345612/92‑1 (MQ=255)
      tAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTAcc                                                                                   >  1:3607298/1‑93 (MQ=255)
       aaCGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCg                                                                                  <  2:6064454/93‑1 (MQ=255)
       aaCGCCGACCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCg                                                                                  >  2:7761713/1‑93 (MQ=255)
        aCGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGa                                                                                 <  1:6513979/93‑1 (MQ=255)
        aCGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGa                                                                                 <  1:2537562/93‑1 (MQ=255)
         acccgccCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGat                                                                                <  1:6851246/91‑1 (MQ=255)
          gccgccCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGatc                                                                               <  1:4318002/93‑1 (MQ=255)
          gccgccCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGatc                                                                               <  1:7527013/93‑1 (MQ=255)
           ccgccCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGatca                                                                              <  2:1868681/93‑1 (MQ=255)
             gccCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGatcatc                                                                            <  1:2352932/93‑1 (MQ=255)
               ccccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCaa                                                                          >  2:4692832/1‑93 (MQ=255)
                cccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAAt                                                                         >  1:3253392/1‑93 (MQ=255)
                cccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAAt                                                                         <  1:1318974/93‑1 (MQ=255)
                 ccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATa                                                                        >  1:7128720/1‑93 (MQ=255)
                 ccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATa                                                                        >  2:1054703/1‑93 (MQ=255)
                 ccTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATa                                                                        >  2:3463635/1‑93 (MQ=255)
                  cTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATaa                                                                       <  2:4892987/93‑1 (MQ=255)
                  cTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATaa                                                                       <  2:4486172/93‑1 (MQ=255)
                   tGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAAc                                                                      >  1:7243367/1‑93 (MQ=255)
                   tGTAACAGCAAGGTGCGCCAGCCACGCTGATTGCGATAGCGCGAGATCAGCGCCACTGGCAGTAACGTCAGCAAACTACCGATCATCAATaaa                                                                      >  2:215275/1‑92 (MQ=255)
                    gTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACg                                                                     <  2:6661110/93‑1 (MQ=255)
                    gTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGACCAGCGCTACGGGCAGTAACGTCAGCAAACCACAGATCATCAATAACg                                                                     >  1:2004618/1‑93 (MQ=255)
                     tAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACgg                                                                    >  1:1006402/1‑93 (MQ=255)
                     tAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACgg                                                                    <  2:2372830/93‑1 (MQ=255)
                     tAACAGCAAGGTGCGCCAGACACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACgg                                                                    <  2:5229041/93‑1 (MQ=255)
                      aaCAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGa                                                                   <  2:2858157/93‑1 (MQ=255)
                         aGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGaaaa                                                                >  1:276777/1‑93 (MQ=255)
                          gCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAAt                                                               >  1:4003257/1‑93 (MQ=255)
                          gCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAAt                                                               >  1:6311340/1‑93 (MQ=255)
                            caGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAGCGTCAGCAAACTACCGATCATCAATAACGGAAAATgg                                                             <  1:7803760/92‑1 (MQ=255)
                            aaGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATgg                                                             >  1:3934985/1‑93 (MQ=255)
                            aaGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATgg                                                             >  1:2896599/1‑93 (MQ=255)
                            aaGGTGCGCCAGTCACGCTGATGACGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATgg                                                             <  2:5406850/93‑1 (MQ=255)
                             aGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATggg                                                            >  1:1439891/1‑93 (MQ=255)
                              ggTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGt                                                           >  1:4871353/1‑93 (MQ=255)
                              ggTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGt                                                           >  1:423254/1‑93 (MQ=255)
                               gTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGtt                                                          <  2:2710146/93‑1 (MQ=255)
                               gTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGtt                                                          <  2:3141852/93‑1 (MQ=255)
                                  cgcCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTaaa                                                       <  1:4112397/93‑1 (MQ=255)
                                   gcCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAac                                                      <  1:3073368/93‑1 (MQ=255)
                                   gcCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAac                                                      <  1:5695295/93‑1 (MQ=255)
                                    ccAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAacc                                                     >  1:5596186/1‑93 (MQ=255)
                                     cAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccg                                                    >  2:1167013/1‑93 (MQ=255)
                                      aGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccga                                                   >  1:3676424/1‑93 (MQ=255)
                                       gTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccgac                                                  <  2:7631229/93‑1 (MQ=255)
                                        tCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAaccgacc                                                 >  1:4592386/1‑93 (MQ=255)
                                           cGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGcc                                              <  2:4308296/93‑1 (MQ=255)
                                            gCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCt                                             >  1:4793048/1‑93 (MQ=255)
                                             cTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTg                                            >  1:4485947/1‑93 (MQ=255)
                                             cTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTg                                            >  2:6997665/1‑93 (MQ=255)
                                              tGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGCCCGCCtcc                                           >  1:5287366/1‑91 (MQ=255)
                                               gATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCg                                          <  1:160376/93‑1 (MQ=255)
                                                aTGGCGATAGCGCGAGATGAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAGCCGACCGCCTGCGt                                         <  2:6453868/93‑1 (MQ=255)
                                                aTGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGt                                         >  1:3467303/1‑93 (MQ=255)
                                                 tGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGtt                                        <  2:5526816/93‑1 (MQ=255)
                                                 tGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGtt                                        <  2:2987167/93‑1 (MQ=255)
                                                     gATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGcc                                    >  1:7610370/1‑93 (MQ=255)
                                                      aTAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAACCTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCa                                   <  2:7794603/93‑1 (MQ=255)
                                                      aTAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCa                                   >  1:531181/1‑93 (MQ=255)
                                                       tAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAg                                  <  2:4409248/93‑1 (MQ=255)
                                                        aGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAACCCGACCGCCTGCGTTAGCCAGt                                 <  1:3803077/93‑1 (MQ=255)
                                                           gcgAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTaaa                              <  2:7678641/93‑1 (MQ=255)
                                                           gcgAGATCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTaaa                              >  2:2074593/1‑93 (MQ=255)
                                                             gagaTCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATc                            <  1:59130/93‑1 (MQ=255)
                                                              agaTCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATca                           >  1:7460720/1‑93 (MQ=255)
                                                              agaTCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATca                           <  1:1055669/93‑1 (MQ=255)
                                                              agaTCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATca                           <  1:5228708/93‑1 (MQ=255)
                                                                aTCAGCGCTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATcagc                         <  2:2522797/93‑1 (MQ=255)
                                                                  cAGCGTTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATcagcag                       <  1:6641226/93‑1 (MQ=255)
                                                                    gcgcTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGcc                     >  2:3844676/1‑93 (MQ=255)
                                                                    gcgcTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGcc                     >  1:4997280/1‑93 (MQ=255)
                                                                      gcTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCAc                   <  1:6918568/93‑1 (MQ=255)
                                                                       cTACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACt                  >  2:7386510/1‑93 (MQ=255)
                                                                        tACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTc                 >  1:6767973/1‑93 (MQ=255)
                                                                        tACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTc                 >  1:5093906/1‑93 (MQ=255)
                                                                        tACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTc                 >  1:1130083/1‑93 (MQ=255)
                                                                        tACGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTc                 >  2:680649/1‑93 (MQ=255)
                                                                         aCGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTcc                <  2:7151268/93‑1 (MQ=255)
                                                                         aCGGGCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTcc                <  1:1699733/93‑1 (MQ=255)
                                                                          cGGGCAGTAACGTCAGCACACTACCGCTCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTccc               >  2:5951144/1‑92 (MQ=255)
                                                                             gCAGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCgcgc            <  2:6258596/93‑1 (MQ=255)
                                                                               aGTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGccc          >  2:5368940/1‑93 (MQ=255)
                                                                               aGTAACGTCAGCAAACTACCCATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGccc          >  2:1589358/1‑93 (MQ=255)
                                                                                gTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCa         >  2:6363418/1‑93 (MQ=255)
                                                                                gTAACGTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCa         >  1:302379/1‑93 (MQ=255)
                                                                                  aaCGTCAGCAATCTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAga       >  2:3999594/1‑92 (MQ=255)
                                                                                     gTCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAAt    >  1:3385572/1‑93 (MQ=255)
                                                                                      tCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATa   <  1:1822323/93‑1 (MQ=255)
                                                                                      tCAGCAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATa   >  1:6092636/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGCCGTTAACGCCGCCCCCTGTAACAGCAAGGTGCGCCAGTCACGCTGATGGCGATAGCGCGAGATCAGCGCTACGGGCAGTAACGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATAC  >  NC_000913/3847401‑3847580

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