Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,516,689 T→C T382A (ACG→GCG)  yfeA ← putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,516,6890TC97.0% 338.7 / ‑0.3 102T382A (ACG→GCG) yfeAputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (50/48);  total (52/50)
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.26e-01

Only 100 of 102 total aligned reads displayed.

GTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516593‑2516784
                                                                                                |                                                                                               
gTCGCCCTTCCGTTTCCGGCCCGCCCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACg                                                                                              <  2:2320087/100‑1 (MQ=255)
  tgctcTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            <  1:2188824/96‑1 (MQ=255)
  cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            >  2:1817654/1‑100 (MQ=255)
  cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            >  1:1578589/1‑100 (MQ=255)
    cccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaa                                                                                          >  2:2895712/1‑100 (MQ=255)
     ccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACAGCATTATCGCACGGCaaa                                                                                         >  2:2001321/1‑100 (MQ=255)
      cTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAt                                                                                        >  1:205681/1‑100 (MQ=255)
        tCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTg                                                                                      >  1:237098/1‑100 (MQ=255)
        tCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTg                                                                                      >  2:2365387/1‑100 (MQ=255)
          cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAg                                                                                    <  2:1500711/100‑1 (MQ=255)
           gTTTCCGGCCCGCTCAGCGCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                   <  2:1996569/100‑1 (MQ=255)
           gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                   >  2:2832916/1‑100 (MQ=255)
            tttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCg                                                                                  <  2:1100547/100‑1 (MQ=255)
            tttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCg                                                                                  <  1:3312035/100‑1 (MQ=255)
              tCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGa                                                                                <  1:2991972/100‑1 (MQ=255)
               ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGa                                                                                <  2:667084/99‑1 (MQ=255)
               ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                               <  1:3270159/100‑1 (MQ=255)
                cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              >  1:1364579/1‑100 (MQ=255)
                cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              >  1:2118672/1‑100 (MQ=255)
                 ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                             <  1:1086782/100‑1 (MQ=255)
                 ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                             >  1:587234/1‑100 (MQ=255)
                  gCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                             <  2:1997137/99‑1 (MQ=255)
                   cccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAg                                                                           >  1:1090555/1‑100 (MQ=255)
                    ccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGt                                                                          <  1:871078/100‑1 (MQ=255)
                     cGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTg                                                                         <  1:1261127/100‑1 (MQ=255)
                       cTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGaa                                                                       >  1:2365240/1‑100 (MQ=255)
                           gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   <  2:1300096/100‑1 (MQ=255)
                           gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   >  2:1112157/1‑100 (MQ=255)
                             aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 <  2:2210445/100‑1 (MQ=255)
                                aaCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAActgctgtta                                                              >  2:1930486/1‑92 (MQ=255)
                                aaCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatca                                                              >  2:3107841/1‑100 (MQ=255)
                                 aCAGCCGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatcat                                                             <  1:2832929/100‑1 (MQ=255)
                                  cagcagTTCACTCCCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGCTACAGTGAACGCGCATCAtt                                                            <  2:2493432/100‑1 (MQ=255)
                                  cagcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCAtt                                                            >  1:1012599/1‑100 (MQ=255)
                                   agcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTa                                                           <  1:2217385/100‑1 (MQ=255)
                                    gcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTaa                                                          >  2:1189994/1‑100 (MQ=255)
                                      agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                        <  2:688919/100‑1 (MQ=255)
                                      agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                        >  1:2652614/1‑100 (MQ=255)
                                      agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                        <  1:3256328/100‑1 (MQ=255)
                                       gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                       >  1:238282/1‑100 (MQ=255)
                                        ttCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg                                                      <  2:846274/100‑1 (MQ=255)
                                         tCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGt                                                     >  1:851041/1‑100 (MQ=255)
                                               cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               >  2:7363/1‑100 (MQ=255)
                                                ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACCGTGAACGCGCATCATTAAGCCGTAATgacg                                              <  1:2365398/100‑1 (MQ=255)
                                                ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTTAGCCGTAATgacg                                              <  1:2279023/100‑1 (MQ=255)
                                                ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacg                                              >  2:3228275/1‑100 (MQ=255)
                                                ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacg                                              >  1:2819287/1‑100 (MQ=255)
                                                  ggCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacgac                                            <  1:2514200/100‑1 (MQ=255)
                                                  ggCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacgac                                            >  1:167116/1‑100 (MQ=255)
                                                  ggCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacgac                                            >  2:109738/1‑100 (MQ=255)
                                                   gCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACt                                           <  2:238281/100‑1 (MQ=255)
                                                   gCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACt                                           >  2:2093260/1‑100 (MQ=255)
                                                   gCAATTGATACAACTTTTCGTTTTCCCGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACt                                           <  2:1648390/100‑1 (MQ=255)
                                                    cAATTTATACAACTTTTCGTTTTCCTGCCACAACGGCATCGGCGAACGGCAAAGTGAGCGGATACATTTAACGCTCTTCATTAAGCCGTAATGACGACTc                                          <  2:1626023/100‑1 (MQ=255)
                                                     aaTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCa                                         <  1:3322597/100‑1 (MQ=255)
                                                     aaTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCa                                         >  1:2893366/1‑100 (MQ=255)
                                                      aTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCAt                                        <  2:1336385/100‑1 (MQ=255)
                                                           tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   >  1:562451/1‑100 (MQ=255)
                                                           tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   >  1:2136400/1‑100 (MQ=255)
                                                           tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   >  1:3110811/1‑100 (MQ=255)
                                                             cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                 >  2:129576/1‑100 (MQ=255)
                                                               aCTTTTCGTTTTCCTGCATCAACGTCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                               <  2:465685/100‑1 (MQ=255)
                                                               aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                               >  1:2773961/1‑100 (MQ=255)
                                                               aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                               >  1:293393/1‑100 (MQ=255)
                                                                cTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGa                              >  1:1698456/1‑100 (MQ=255)
                                                                cTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGa                              >  2:2328906/1‑100 (MQ=255)
                                                                   ttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTa                           >  1:3142401/1‑100 (MQ=255)
                                                                   tccGTTTTCCGCCATAACCGCCGGAACCGCACGCCAAATTAAGCGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTAt                          <  1:45427/98‑1 (MQ=255)
                                                                    tCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTAt                          >  2:506713/1‑100 (MQ=255)
                                                                     ccttttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATc                         <  2:2028476/98‑1 (MQ=255)
                                                                     cGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATc                         <  1:1498300/100‑1 (MQ=255)
                                                                     cGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATc                         <  1:942722/100‑1 (MQ=255)
                                                                     cGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATc                         <  1:506714/100‑1 (MQ=255)
                                                                      gTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCa                        <  2:3236431/100‑1 (MQ=255)
                                                                      gTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCa                        >  2:2810585/1‑100 (MQ=255)
                                                                      gTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCa                        >  2:2792604/1‑100 (MQ=255)
                                                                        tttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAAt                      <  1:122660/100‑1 (MQ=255)
                                                                         ttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATg                     <  1:2903519/100‑1 (MQ=255)
                                                                         ttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATg                     <  1:1142164/100‑1 (MQ=255)
                                                                           ccTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAAt                      >  2:541411/1‑97 (MQ=255)
                                                                            cggCATAAACGGCAGCATCGCACGGCAAAGTGAGCGGATACAGTGAACGCCCATCATTAATCCGTAATGACGCAGCATATTCTCAACATTATCa                        <  2:758006/92‑1 (MQ=255)
                                                                            cTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATgcgc                  >  2:1549430/1‑100 (MQ=255)
                                                                             tGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgca                 >  1:2516734/1‑100 (MQ=255)
                                                                              gCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcag                >  2:2650636/1‑100 (MQ=255)
                                                                              gCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcag                <  1:1112158/100‑1 (MQ=255)
                                                                              gCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcag                >  1:3172878/1‑100 (MQ=255)
                                                                               cATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcagg               <  2:1090554/100‑1 (MQ=255)
                                                                                aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc              >  2:2291098/1‑100 (MQ=255)
                                                                                aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc              >  1:2141260/1‑100 (MQ=255)
                                                                                aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc              <  1:3096518/100‑1 (MQ=255)
                                                                                  gatctgcagcaGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggcag            <  1:12842/95‑1 (MQ=255)
                                                                                  cAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggcag            <  1:1940594/100‑1 (MQ=255)
                                                                                  cAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggcag            <  2:2118662/100‑1 (MQ=255)
                                                                                   aaCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagc           >  2:1042/1‑100 (MQ=255)
                                                                                      ggCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaa        <  1:3269922/100‑1 (MQ=255)
                                                                                       gcagcaGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggcagagatc       >  2:12842/1‑95 (MQ=255)
                                                                                        cagcagCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaaa      >  1:216947/1‑100 (MQ=255)
                                                                                          ncagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAct    >  2:31129/2‑100 (MQ=255)
                                                                                          gcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAct    <  1:1189998/100‑1 (MQ=255)
                                                                                            agcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActct  <  2:478680/100‑1 (MQ=255)
                                                                                                |                                                                                               
GTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516593‑2516784

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