Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 732,503 G→A 9.0% R974Q (CGG→CAG)  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,5030GA 9.0% 222.1 / 11.8 78R974Q (CGG→CAG) rhsCRhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (33/38);  new base A (3/4);  total (36/42)
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.09e-01

AACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAA  >  NC_000913/732361‑732651
                                                                                                                                              |                                                                                                                                                    
aaCCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGt                                                                                                                                                >  1:150044/1‑149 (MQ=32)
 aCCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGGCAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCGCCGGTACCATTACGACAGTCAGCACCGGCTGGTg                                                                                                                                               <  1:498947/149‑1 (MQ=17)
 aCCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTg                                                                                                                                               <  1:313081/149‑1 (MQ=32)
     tATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACt                                                                                                                                           >  1:61767/1‑149 (MQ=32)
      aTCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCAGCCGGTGctgac                                                                                                                                          <  1:531095/149‑5 (MQ=255)
        cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTaca                                                                                                                                        <  2:118575/149‑1 (MQ=32)
          cccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACAcg                                                                                                                                      >  4:181631/1‑149 (MQ=32)
           ccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACAcgc                                                                                                                                     <  3:350759/149‑1 (MQ=32)
             gTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCgg                                                                                                                                   >  1:98051/1‑149 (MQ=32)
                   cgcACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACaa                                                                                                                             <  2:314078/149‑1 (MQ=32)
                    gcACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGGTGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAat                                                                                                                            >  1:180886/1‑149 (MQ=17)
                     cACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAata                                                                                                                           <  2:150044/149‑1 (MQ=32)
                     cACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAata                                                                                                                           >  2:610781/1‑149 (MQ=32)
                     cACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAata                                                                                                                           <  1:104639/149‑1 (MQ=32)
                      aCTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAGGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAatat                                                                                                                          >  2:82160/1‑149 (MQ=17)
                      aCTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAatat                                                                                                                          >  4:492992/1‑149 (MQ=32)
                       cTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCgg                                                                                                                                   >  4:174711/1‑139 (MQ=32)
                       cTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCgg                                                                                                                                   <  3:174711/139‑1 (MQ=32)
                        tATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGc                                                                                                                        >  4:273778/1‑149 (MQ=32)
                        tATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGc                                                                                                                        <  3:396437/149‑1 (MQ=32)
                         aTCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCa                                                                                                                       <  4:491149/149‑1 (MQ=32)
                          tCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCag                                                                                                                      <  4:85142/149‑1 (MQ=32)
                                   gTATGACCGTCACGGCAGGCTGACAGAGATAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTgg                                                                                                             >  3:505448/1‑149 (MQ=17)
                                           gTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGCTGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATCTGCAGAGCCGCTGGTCGAAAGt                                                                                                     >  2:198851/1‑149 (MQ=17)
                                           gTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGt                                                                                                     <  4:588777/149‑1 (MQ=32)
                                           gTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGt                                                                                                     >  3:361052/1‑149 (MQ=33)
                                            tCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTc                                                                                                    <  4:368572/149‑1 (MQ=34)
                                               cGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCt                                                                                                 <  3:621536/149‑1 (MQ=34)
                                                   aGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCt                                                                                             <  3:76820/149‑1 (MQ=35)
                                                      cTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTa                                                                                          <  3:181631/149‑1 (MQ=34)
                                                         aCAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGa                                                                                       <  2:64923/149‑1 (MQ=34)
                                                           agagAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGAcc                                                                                     >  2:472770/1‑149 (MQ=32)
                                                           agagAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGAcc                                                                                     >  3:480395/1‑149 (MQ=32)
                                                                 aaCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTgg                                                                               >  3:279635/1‑149 (MQ=32)
                                                                    cGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTa                                                                                          >  3:42016/1‑135 (MQ=32)
                                                                    cGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTa                                                                                          <  4:42016/135‑1 (MQ=32)
                                                                        cTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAg                                                                        <  2:222125/149‑1 (MQ=32)
                                                                         tCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAgg                                                                       >  2:73336/1‑149 (MQ=32)
                                                                           aTCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGt                                                                     >  4:20591/1‑149 (MQ=32)
                                                                              ccGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGc                                                                  >  2:436033/1‑149 (MQ=32)
                                                                                      ggTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATTCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACggg                                                          <  4:424467/149‑1 (MQ=17)
                                                                                       gTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGt                                                         >  4:9317/1‑149 (MQ=32)
                                                                                       gTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGt                                                         <  1:621454/149‑1 (MQ=32)
                                                                                          aTCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATg                                                      <  1:172707/149‑1 (MQ=35)
                                                                                           tCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATgg                                                     <  1:131404/149‑1 (MQ=34)
                                                                                              gCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTAc                                                                                         <  2:184148/110‑1 (MQ=33)
                                                                                              gCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTAc                                                                                         >  1:184148/1‑110 (MQ=32)
                                                                                              gCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATggcgg                                                  >  4:301601/1‑149 (MQ=35)
                                                                                               cACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATggcggc                                                 <  1:610781/149‑1 (MQ=35)
                                                                                                  ggATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTg                                              >  1:413465/1‑149 (MQ=37)
                                                                                                   gatgatGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGa                                             >  2:126644/1‑149 (MQ=37)
                                                                                                      gatgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACg                                          >  4:108587/1‑149 (MQ=37)
                                                                                                       atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACgg                                         <  1:346370/149‑1 (MQ=37)
                                                                                                       atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACgg                                         <  4:480395/149‑1 (MQ=37)
                                                                                                          aGCGCACCCACCGGTACCATTACGACAGTCAGCATCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGAc                                      <  1:389479/149‑1 (MQ=38)
                                                                                                          aGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGAc                                      <  2:180886/149‑1 (MQ=38)
                                                                                                          aGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGAc                                      <  3:563113/149‑1 (MQ=38)
                                                                                                           gcgcACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGAcc                                     >  4:90725/1‑149 (MQ=39)
                                                                                                              cacccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGa                                  <  2:590334/149‑1 (MQ=38)
                                                                                                                 ccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACgg                               >  4:87070/1‑149 (MQ=35)
                                                                                                                    ccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCt                            <  1:126644/149‑1 (MQ=35)
                                                                                                                        tACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGAt                        >  3:239136/1‑149 (MQ=35)
                                                                                                                        tACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGAt                        <  2:577585/149‑1 (MQ=35)
                                                                                                                            aTTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCg                    <  3:273778/149‑1 (MQ=34)
                                                                                                                              tACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCt                  <  3:87070/149‑1 (MQ=35)
                                                                                                                              tACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCt                  >  2:248827/1‑149 (MQ=35)
                                                                                                                                 gACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGaa                                            <  1:494122/120‑1 (MQ=38)
                                                                                                                                 gACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGaa                                            >  2:494122/1‑120 (MQ=38)
                                                                                                                                  aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCa              <  1:472770/149‑1 (MQ=34)
                                                                                                                                     tagtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATg                       <  1:429596/133‑1 (MQ=255)
                                                                                                                                     tagtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATg                       >  2:429596/5‑137 (MQ=255)
                                                                                                                                      tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGa          >  3:222535/1‑149 (MQ=35)
                                                                                                                                      agtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTa                                                                                          <  3:476189/66‑1 (MQ=11)
                                                                                                                                      agtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTa                                                                                          >  4:476189/4‑69 (MQ=11)
                                                                                                                                       gtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTAc                                                                                         >  4:571005/3‑69 (MQ=11)
                                                                                                                                       gtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTAc                                                                                         <  3:571005/67‑1 (MQ=11)
                                                                                                                                         gCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAAc       <  2:169425/149‑1 (MQ=35)
                                                                                                                                         gCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAAc       >  3:91381/1‑149 (MQ=34)
                                                                                                                                              ggCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCaa  <  4:279635/149‑1 (MQ=35)
                                                                                                                                              |                                                                                                                                                    
AACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAA  >  NC_000913/732361‑732651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: