Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 732,506 T→C 17.4% L975P (CTG→CCG)  rhsC → Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913732,5060TC17.4% 185.6 / 35.3 97L975P (CTG→CCG) rhsCRhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (42/38);  new base C (11/6);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

Only 100 of 103 total aligned reads displayed.

ACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCG  >  NC_000913/732407‑732605
                                                                                                   |                                                                                                   
aCGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCt                                                                                                     <  1:1365461/100‑1 (MQ=33)
 cgtcAGTCTGACATAGAAATCCTACCTCATCCCTTAATATGTTATCCGCACGTATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCt                                                                                                     <  1:3330740/96‑1 (MQ=11)
 cGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTg                                                                                                    <  2:2784167/100‑1 (MQ=32)
 cGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTg                                                                                                    >  2:1782385/1‑100 (MQ=32)
  ggCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTgg                                                                                                   <  1:2207051/100‑1 (MQ=32)
   gCAGGCTGACAGAGAAACCCGCCCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCGCCCACCGGTACCATTACGACAGTCAGCACCAGCCGGt                                                                                                  <  2:1351606/100‑1 (MQ=11)
     aGGCTGACAGAGAAAACCGACCTCATCCCGGGAGGGGTTATCCGAACGGATGATGAGCGCACCCACCGGTACCATTACTACTGTCAGCACCGGCTGGTGc                                                                                                <  2:1239351/100‑1 (MQ=11)
     aGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGc                                                                                                <  1:2363776/100‑1 (MQ=33)
      ggCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCa                                                                                               >  1:851790/1‑100 (MQ=33)
       gCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATTATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCAc                                                                                              <  2:3321556/100‑1 (MQ=18)
       gCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCAc                                                                                              <  1:3389733/100‑1 (MQ=33)
       gCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCAc                                                                                              >  2:1599915/1‑100 (MQ=32)
        cTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACt                                                                                             >  2:1770877/1‑100 (MQ=33)
        cTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACt                                                                                             <  1:85476/100‑1 (MQ=33)
        cTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACt                                                                                             >  1:598766/1‑100 (MQ=33)
        cTGACAGAGAAAACCGACCTCATCCCCGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACt                                                                                             <  2:3784481/100‑1 (MQ=17)
         tGACAGAGAAAACTGACCTCATCCCGGAAGGGGTTATCCGCACGGATGAGGAGCGCACCCACCGGTACGTTTACGACAGTCAGCACCGGCTGGTGCACTa                                                                                            <  2:3423605/100‑1 (MQ=11)
         tGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTa                                                                                            >  1:3091283/1‑100 (MQ=32)
         tGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCAGCCGGTGctgac                                                                                            <  2:1769247/100‑5 (MQ=11)
         tGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCAGCCGGTGctgac                                                                                            <  1:891012/100‑5 (MQ=11)
          gACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTac                                                                                           >  1:3928669/1‑100 (MQ=32)
          gACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCAGCCGGTGctgact                                                                                           <  2:3609670/100‑6 (MQ=11)
           aCAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCTCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCAGCCGGTGctgactg                                                                                          <  1:3157424/100‑7 (MQ=11)
           aCAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCTCCGGTACCATTACGACAGTCAGCACCAGCCGGTGctgactg                                                                                          <  2:1249446/100‑7 (MQ=11)
           aCAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTaca                                                                                          <  1:2944829/100‑1 (MQ=32)
            cAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTacac                                                                                         <  1:1222657/100‑1 (MQ=31)
            cAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTacac                                                                                         <  2:2251303/100‑1 (MQ=31)
            cAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTacac                                                                                         <  1:1665732/100‑1 (MQ=31)
            cAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTacac                                                                                         <  2:3662955/100‑1 (MQ=31)
             agagAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATAATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACAcg                                                                                        >  1:85686/1‑100 (MQ=11)
                                               cGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGt                                                       <  1:3742046/99‑1 (MQ=11)
                                               cGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTc                                                      >  1:2080201/1‑100 (MQ=31)
                                               cGCACGGATGATAAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACAAAATATGCAGAGCCGCTGGTCGAAAGTc                                                      >  2:1266373/1‑100 (MQ=9)
                                                 gaCGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCATAGCCGCTGGTCGAAAGTCg                                                     <  2:1983791/98‑1 (MQ=17)
                                                  aCGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGc                                                    <  2:2428527/99‑1 (MQ=32)
                                                      atgatgAGCGCACCCACCGGTACCATTACTACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCt                                               <  1:3425670/100‑1 (MQ=18)
                                                      atgatgAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCt                                               >  2:3220094/1‑100 (MQ=32)
                                                       tgatgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCtt                                              <  2:1170066/100‑1 (MQ=33)
                                                        gatgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCttt                                             >  1:1498890/1‑100 (MQ=32)
                                                        gatgaGCGCACACACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCtct                                             >  2:43158/1‑98 (MQ=18)
                                                         atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGCAACACTATGCAGAGCCGCTGGTCGAAGGTCGCTATGttt                                             >  2:2453726/1‑99 (MQ=11)
                                                         atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGTTCGAAAGTCGCTATCtt                                              <  2:1408807/98‑1 (MQ=17)
                                                         atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTa                                            <  1:3193648/100‑1 (MQ=33)
                                                         atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTa                                            <  1:2051309/100‑1 (MQ=33)
                                                          tgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCCCGGCGAAAGTCGCTATCTTTAc                                           >  2:3425533/1‑100 (MQ=255)
                                                           gaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACg                                          <  2:393253/100‑1 (MQ=33)
                                                            aGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGa                                         >  1:61408/1‑100 (MQ=33)
                                                            aGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGa                                         >  1:343514/1‑100 (MQ=32)
                                                            aGCGCACCCACCGGTACCATTACGACAGTCAGCACAGGCTGGTGCACTACACGCGGCCACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGa                                         >  2:1385736/1‑100 (MQ=12)
                                                            aGCGCACCCACCGGTAACATTACGACAGTCAGCACCGGCTGGTGCACTACACCCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGa                                         >  2:14617/1‑100 (MQ=12)
                                                             gcgcACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATACAAAACCACTGGTCGAAAGTCGCTATCTTTACGAc                                        >  1:1962829/1‑100 (MQ=11)
                                                              cgcACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGCCACAATATGCAGAGCCGCTGGCCGAAAGTCGCTATCTTTACGAcc                                       >  2:3468219/1‑100 (MQ=12)
                                                                cacccaccGGTACCATTACGACAGTCAGCACCGCCTGGTGCACTACACCCGGACACAATATGCACAGCTCCTACGCGAAAGTCGCTATCTTTACGAcc                                       >  2:1457147/1‑98 (MQ=11)
                                                                 acccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGc                                    >  2:3857307/1‑100 (MQ=32)
                                                                  cccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCg                                     >  2:2910973/1‑98 (MQ=32)
                                                                  cccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCt                                   <  2:1197891/100‑1 (MQ=33)
                                                                    ccccGGTACCATTACGACAGTCAGCACAGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTg                                  >  2:2788864/3‑99 (MQ=17)
                                                                     accGGTACCATTACGACAGTCAGCACCGGCTGGTGCATTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTggg                                <  1:2067729/100‑1 (MQ=17)
                                                                      ccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGc                               <  1:86270/100‑1 (MQ=32)
                                                                      ccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGc                               <  1:2864314/100‑1 (MQ=33)
                                                                          tACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACAAAATATGCAGAGCCGTTGGTCGAAAGTCGCTATCTTTATGACTCGCTGGGCCGCa                           <  2:1601488/100‑1 (MQ=11)
                                                                            ccATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAgg                         >  2:2229181/1‑100 (MQ=33)
                                                                             cATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAggg                        >  2:3008912/1‑100 (MQ=32)
                                                                             cATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAggg                        <  1:2053596/100‑1 (MQ=32)
                                                                               ttACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTg                      <  2:409531/100‑1 (MQ=33)
                                                                                tACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTg                      >  2:1738316/1‑99 (MQ=32)
                                                                                 aCGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGc                    >  2:3616223/1‑100 (MQ=32)
                                                                                  cGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGc                    >  2:2856794/1‑99 (MQ=33)
                                                                                  cGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCa                   <  1:3439274/100‑1 (MQ=33)
                                                                                  cGACAGTCAGCACCGGCTGGTGAACTACACGCGGACACTATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGc                    <  2:3666535/99‑1 (MQ=12)
                                                                                  cGACACTCAGCACCGGCTGGTGCACTACAGTGGGACATAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGc                    <  2:3915079/99‑1 (MQ=11)
                                                                                   cgtgtagtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCaa                  >  1:1789544/9‑100 (MQ=11)
                                                                                    gtgtagtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACACGCTGGGCCGCAGGGTGGCaaa                 >  1:830515/8‑100 (MQ=11)
                                                                                    aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTAGCTATCTTTACGACAAGCTAGGCCACAGGGTGGCaaa                 >  1:1894112/1‑100 (MQ=11)
                                                                                     tgtagtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCaaaa                >  2:1012741/7‑100 (MQ=11)
                                                                                     cAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCaaaa                >  1:2924075/1‑100 (MQ=33)
                                                                                      gtagtgCACCAGCCGGTGCACTACACGCGGAAACAATATGCAGAGACGCTGGTCGAAAGTCGCGATCTTAACGACACGCTGGGCCGCAGGGGGGCAAAAc               >  2:786942/6‑100 (MQ=11)
                                                                                      gtagtgCACCAGCCGGTGCACTAAAAGAGGACACAATATGCAGAGACGATGGTAGAAAGTCGCTATCTTTACGAAACGCTGGGCCGGAGGGTGGCaa                  >  2:728097/6‑97 (MQ=2)
                                                                                      aGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAAc               <  1:292898/100‑1 (MQ=32)
                                                                                       tagtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACg              >  1:1963115/5‑100 (MQ=11)
                                                                                       tagtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAAGGTGGCAAAAc               >  1:2681838/5‑99 (MQ=11)
                                                                                       gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACg              >  2:2905514/1‑100 (MQ=32)
                                                                                        agtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg             >  2:1020391/4‑100 (MQ=255)
                                                                                        agtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg             >  2:509615/4‑100 (MQ=255)
                                                                                        agtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg             >  1:2335715/4‑100 (MQ=255)
                                                                                        agtgCACCAGCCGGTGCACTACACGCAGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGAACCGCTGGGCCGCAGGGTGGCAAAACg              >  1:2375987/4‑99 (MQ=11)
                                                                                         cAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACggg            >  2:3573219/1‑100 (MQ=32)
                                                                                          aGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGt           <  2:2071268/100‑1 (MQ=33)
                                                                                           gCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACggg            >  2:3923891/1‑98 (MQ=32)
                                                                                           gCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTa          >  2:2087615/1‑100 (MQ=35)
                                                                                            cACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTAt         >  2:3138257/1‑100 (MQ=35)
                                                                                             aCCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATg        <  1:2802687/100‑1 (MQ=35)
                                                                                              ccGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATgg       >  1:1092849/1‑100 (MQ=35)
                                                                                              ccGGCTGGTGCACTACACACGTACACAATATACAGAGCCGCTAGTCGAAAGTCGCTATCTTTAGGCCCAGCCGGGCCGCAGGGTGGCCAAAGGGGTATg        >  2:1825461/1‑99 (MQ=11)
                                                                                                 gCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATggcgg    <  1:3864232/100‑1 (MQ=38)
                                                                                                 gCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGAACCGCTGGGCCGCAGGGTGGCAAAACGGGTATggcgg    >  1:1439960/1‑100 (MQ=25)
                                                                                                   tGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGgcggcg  <  2:3018422/100‑1 (MQ=35)
                                                                                                   tGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGgcggcg  >  1:3640021/1‑100 (MQ=38)
                                                                                                   tGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGgcggcg  <  2:143769/100‑1 (MQ=38)
                                                                                                   tGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTAGGCCGCAGGGTGGCAAAACGGGTATgacgg    >  2:3491810/1‑98 (MQ=255)
                                                                                                   |                                                                                                   
ACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCG  >  NC_000913/732407‑732605

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