Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913730,6440AG50.7% 22.0 / 104.5 73R354R (CGA→CGGrhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base A (15/21);  new base G (18/19);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTG  >  NC_000913/730511‑730779
                                                                                                                                     |                                                                                                                                       
cgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGaa                                                                                                                                   <  1:338181/140‑1 (MQ=255)
  cTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGTGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAAt                                                                                                                                 >  1:195756/1‑140 (MQ=255)
             tGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCACACCGTCACACGGGCCGACCgg                                                                                                                                     <  2:189265/125‑1 (MQ=255)
             tGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCACACCGTCACACGGGCCGACCgg                                                                                                                                     >  1:189265/1‑125 (MQ=255)
                     cGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGc                                                                                                              >  2:111760/1‑140 (MQ=255)
                             ggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcgggcg                                                                                                      <  1:394920/140‑1 (MQ=255)
                                  aCTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGACGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGa                                                                                                 <  1:129727/140‑1 (MQ=255)
                                      gcggcgGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGa                                                                                             <  1:452301/140‑1 (MQ=255)
                                                    ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCgg                                                                               <  1:111760/140‑1 (MQ=255)
                                                    ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCgg                                                                               <  2:277365/140‑1 (MQ=255)
                                                       cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg                                                                            >  1:516770/1‑140 (MQ=255)
                                                       cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg                                                                            >  1:20734/1‑140 (MQ=255)
                                                         gcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                          <  2:99208/140‑1 (MQ=40)
                                                         gcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                          <  1:441005/140‑1 (MQ=40)
                                                           aaTACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCtt                                                                        <  2:440446/140‑1 (MQ=39)
                                                                  aGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTac                                                                 <  2:516770/140‑1 (MQ=37)
                                                                   ggTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTaca                                                                <  2:274283/140‑1 (MQ=18)
                                                                   ggTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTaca                                                                <  1:88421/140‑1 (MQ=38)
                                                                   ggTGAGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTaca                                                                <  2:500269/140‑1 (MQ=17)
                                                                    gTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTacac                                                               >  1:51339/1‑140 (MQ=17)
                                                                      gcgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGt                                                             >  1:439314/1‑140 (MQ=17)
                                                                       cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTa                                                            >  1:237232/1‑140 (MQ=17)
                                                                       cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTa                                                            <  2:437520/140‑1 (MQ=18)
                                                                       cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTa                                                            >  2:338543/1‑140 (MQ=37)
                                                                       cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTa                                                            <  1:504729/140‑1 (MQ=38)
                                                                       cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTa                                                            >  1:373405/1‑140 (MQ=37)
                                                                       cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTa                                                            >  1:369768/1‑140 (MQ=37)
                                                                        gcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  2:9392/140‑1 (MQ=18)
                                                                        gcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  1:197583/140‑1 (MQ=18)
                                                                        gcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  2:81784/140‑1 (MQ=18)
                                                                        gcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  2:195756/140‑1 (MQ=38)
                                                                        gcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  1:104867/140‑1 (MQ=38)
                                                                        gcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  2:537060/140‑1 (MQ=38)
                                                                        gcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  1:338543/140‑1 (MQ=38)
                                                                        gcAGCTTTACTTACGATGATAAATACCGGGGCCGAATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  1:257071/140‑1 (MQ=18)
                                                                          aGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCa                                                         <  1:275722/140‑1 (MQ=17)
                                                                            cTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGt                                                             <  1:199955/134‑1 (MQ=17)
                                                                            cTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGt                                                             >  2:199955/1‑134 (MQ=17)
                                                                            cTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGt                                                       <  2:468083/140‑1 (MQ=18)
                                                                            cTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGt                                                       <  2:247518/140‑1 (MQ=38)
                                                                                    acGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAg                                               <  1:42478/140‑1 (MQ=18)
                                                                                    acGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAg                                               >  1:260088/1‑140 (MQ=18)
                                                                                    acGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAg                                               >  1:119208/1‑140 (MQ=18)
                                                                                    acGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAg                                               <  2:171643/140‑1 (MQ=18)
                                                                                    acGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAg                                               <  1:455071/140‑1 (MQ=39)
                                                                                      gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGt                                                       >  1:174849/1‑130 (MQ=17)
                                                                                      gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGt                                                       <  2:174849/130‑1 (MQ=18)
                                                                                      gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                             >  1:410867/1‑140 (MQ=21)
                                                                                      gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                             <  2:349755/140‑1 (MQ=25)
                                                                                      gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                             >  1:347828/1‑140 (MQ=21)
                                                                                      gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                             >  1:407085/1‑140 (MQ=40)
                                                                                      gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                             >  1:258335/1‑140 (MQ=40)
                                                                                          atAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGc                                         <  2:463069/140‑1 (MQ=255)
                                                                                           tAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGca                                        >  1:216859/1‑140 (MQ=21)
                                                                                           tAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGca                                        >  1:15652/1‑140 (MQ=21)
                                                                                            aaaTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcat                                       <  1:22799/140‑1 (MQ=25)
                                                                                              aTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatca                                     <  2:222301/140‑1 (MQ=37)
                                                                                                aCCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcaacacc                                   >  1:467230/1‑140 (MQ=21)
                                                                                                    gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                               <  2:216859/140‑1 (MQ=37)
                                                                                                    gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                               >  1:353405/1‑140 (MQ=21)
                                                                                                    gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGACAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                               >  1:460076/1‑140 (MQ=21)
                                                                                                    gggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                               >  1:518549/1‑140 (MQ=40)
                                                                                                      gCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCatcaccatcacc                             >  1:78189/1‑140 (MQ=21)
                                                                                                       ccGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCg                            <  2:22186/140‑1 (MQ=37)
                                                                                                       ccGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCg                            <  2:20734/140‑1 (MQ=255)
                                                                                                        cGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGa                           <  2:347828/140‑1 (MQ=37)
                                                                                                          gATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACa                         >  1:479787/1‑140 (MQ=40)
                                                                                                                 gcgcACCGCCACAGGGGCCGCCGGGAAATTTTTTACCTTTACGACAGCGGGGGGGGGGTGAAAGAACAGCAAACCCCGGAGGGCTAAATTTACACGTATCAGTATGAGAAGGACCGCACCCCCACCCCCGACAGCCTGaa                  <  2:479787/140‑1 (MQ=9)
                                                                                                                  cgcCCCGTCACAGGGGCGGACGGAAAATTTTTTACCGTAACGACGGGGACGGGGGGGTGAAAAAACAGCAAACCCGGGGAGGCTTAAGTTACACGTATGAGTATGAGAGGGCCCGCACCACCACCACCGACGGCCTGAAc                 <  2:353268/140‑1 (MQ=255)
                                                                                                                  cgcACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAc                 >  1:292493/1‑140 (MQ=37)
                                                                                                                    cACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccg               >  1:117737/1‑140 (MQ=255)
                                                                                                                     aCCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccgc              >  1:187144/1‑140 (MQ=37)
                                                                                                                         tCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTg          >  1:329990/1‑140 (MQ=255)
                                                                                                                           acacGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGaa        >  1:534428/1‑140 (MQ=37)
                                                                                                                                 gccggccgGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTg  >  1:236751/1‑140 (MQ=255)
                                                                                                                                 gCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTg  >  1:23084/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTG  >  NC_000913/730511‑730779

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