Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 872,840 A→G 5.4% intergenic (+39/‑139) gsiD → / → yliE glutathione ABC transporter permease/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913872,8400AG 5.4% 238.3 / 3.3 77intergenic (+39/‑139)gsiD/yliEglutathione ABC transporter permease/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (38/33);  new base G (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

TCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTT  >  NC_000913/872701‑872988
                                                                                                                                           |                                                                                                                                                    
tCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAattat                                                                                                                                             >  2:815877/1‑149 (MQ=255)
tCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAattat                                                                                                                                             >  1:670310/1‑149 (MQ=255)
 cGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAattatt                                                                                                                                            <  4:127530/149‑1 (MQ=255)
  gCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTg                                                                                                                                           >  2:655134/1‑149 (MQ=255)
   cTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGc                                                                                                                                          <  1:240757/149‑1 (MQ=255)
         ttCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAg                                                                                                                                    >  2:93991/1‑149 (MQ=255)
            cGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCaa                                                                                                                                 <  4:727191/149‑1 (MQ=255)
                ccTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAgtgt                                                                                                                             >  4:42101/1‑149 (MQ=255)
                ccTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAgtgt                                                                                                                             <  3:855179/149‑1 (MQ=255)
                  tGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACGGGAAattatt                                                                                                                                            <  2:533588/132‑1 (MQ=255)
                  tGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACGGGAAattatt                                                                                                                                            >  3:300192/1‑132 (MQ=255)
                  tGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACGGGAAattatt                                                                                                                                            <  4:300192/132‑1 (MQ=255)
                  tGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACGGGAAattatt                                                                                                                                            >  1:533588/1‑132 (MQ=255)
                    gCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTa                                                                                                                         >  3:858741/1‑149 (MQ=255)
                     cTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTaa                                                                                                                        <  1:655134/149‑1 (MQ=255)
                     cTATTTTTCTGAACGTACTGGCGTTCAATTTTTTGGGCGATGGTTTACGCGATGCGCTGGATCCGGAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTaa                                                                                                                        <  1:258759/149‑1 (MQ=255)
                      tATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAAc                                                                                                                       >  4:118198/1‑149 (MQ=255)
                      tATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAAc                                                                                                                       <  3:853401/149‑1 (MQ=255)
                       aTTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACt                                                                                                                      <  3:570251/149‑1 (MQ=255)
                         ttttCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACttt                                                                                                                    <  1:265918/149‑1 (MQ=255)
                         ttttCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACttt                                                                                                                    <  4:674133/149‑1 (MQ=255)
                                 ccGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGAt                                                                                                            <  3:118198/149‑1 (MQ=255)
                                    tACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCCAAGCAAGTGTTGTAACTTTCTGCTGAtttt                                                                                                         <  4:633676/149‑1 (MQ=255)
                                       tGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTa                                                                                                      <  4:310684/149‑1 (MQ=255)
                                       tGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTa                                                                                                      <  1:815877/149‑1 (MQ=255)
                                           gTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTaa                                                                                                                                     <  3:777357/114‑1 (MQ=255)
                                           gTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTaa                                                                                                                                     >  4:777357/1‑114 (MQ=255)
                                           gTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAAt                                                                                                  <  1:195893/149‑1 (MQ=255)
                                            ttCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAatta                                                                                                                                              <  1:322562/104‑1 (MQ=255)
                                            ttCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAatta                                                                                                                                              >  2:322562/1‑104 (MQ=255)
                                             tCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCg                                                                                                <  4:43526/149‑1 (MQ=255)
                                              cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCgg                                                                                               <  4:488527/149‑1 (MQ=255)
                                               aaTTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCggg                                                                                              >  4:163508/1‑149 (MQ=255)
                                               aaTTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCggg                                                                                              >  4:705931/1‑149 (MQ=255)
                                                    gttgGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAAtt                                                                                         >  1:400503/1‑149 (MQ=255)
                                                        ggCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGt                                                                                     <  2:422595/149‑1 (MQ=255)
                                                        ggCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGt                                                                                     <  2:196843/149‑1 (MQ=255)
                                                            aTGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTaaa                                                                                 >  2:438835/1‑149 (MQ=255)
                                                             tGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTaaaa                                                                                >  3:230854/1‑149 (MQ=255)
                                                                 ttACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGcc                                                                            <  1:438835/149‑1 (MQ=255)
                                                                 ttACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGcc                                                                            >  3:701357/1‑149 (MQ=255)
                                                                  tACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCg                                                                           >  2:683649/1‑149 (MQ=255)
                                                                  tACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCg                                                                           >  2:502194/1‑149 (MQ=255)
                                                                       gATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTaaa                                                                                 <  3:205063/138‑1 (MQ=255)
                                                                       gATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTaaa                                                                                 >  4:205063/1‑138 (MQ=255)
                                                                         tGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACCGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAg                                                                    <  4:693700/149‑1 (MQ=255)
                                                                         tGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAg                                                                    <  4:230854/149‑1 (MQ=255)
                                                                              tGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGccc                                                                                                                                        <  1:497887/76‑1 (MQ=255)
                                                                              tGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGccc                                                                                                                                        >  2:497887/1‑76 (MQ=255)
                                                                               ggATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACaa                                                              <  1:201148/149‑1 (MQ=255)
                                                                                gATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAAt                                                             >  1:395651/1‑149 (MQ=255)
                                                                                     gAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTg                                                        >  3:494450/1‑149 (MQ=255)
                                                                                          ttAAAGGATAGTTACGGTTGAATATTGCTTGAAAGGGTAATCTCCTCAATGGAAATTATTGACCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAATCCTCAGCAAGGGACAATTTTTGTAGCg                                                   >  2:18978/1‑149 (MQ=255)
                                                                                              aGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGcac                                               >  4:386699/1‑149 (MQ=255)
                                                                                               ggATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGcaca                                              >  2:74903/1‑149 (MQ=255)
                                                                                                     ttACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTc                                        >  4:160006/1‑149 (MQ=255)
                                                                                                      tACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCa                                       <  1:103018/149‑1 (MQ=255)
                                                                                                         gTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGAt                                    >  2:443148/1‑149 (MQ=255)
                                                                                                          tttGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCACTAAGCAAGTGGTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCGGGAGGGGTCAGTTTTTGAGGCGGCACAGTGttaaaaaa                                   >  4:756592/1‑143 (MQ=255)
                                                                                                           ttGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAg                                  >  1:777203/1‑149 (MQ=255)
                                                                                                              aaTATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTa                               >  3:27111/1‑149 (MQ=255)
                                                                                                                 atTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTAttt                            >  1:19718/1‑149 (MQ=255)
                                                                                                                 atTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTAttt                            <  2:395651/149‑1 (MQ=255)
                                                                                                                  ttGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTAtttt                           >  3:481683/1‑149 (MQ=255)
                                                                                                                  ttGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTAtttt                           <  2:670310/149‑1 (MQ=255)
                                                                                                                  ttGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTAtttt                           >  3:188426/1‑149 (MQ=255)
                                                                                                                  ttGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTAtttt                           <  4:484046/149‑1 (MQ=255)
                                                                                                                      ttGAAAGGGTAATCACCTTACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTa                       >  4:317000/1‑149 (MQ=255)
                                                                                                                          aaGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGt                                                                                     <  4:282588/83‑1 (MQ=255)
                                                                                                                          aaGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGt                                                                                     >  3:282588/1‑83 (MQ=255)
                                                                                                                            gggTAATTACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTTTAGGATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTa                 <  3:756592/149‑1 (MQ=255)
                                                                                                                                    aCCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTaaaaa                                                                               <  2:274289/79‑1 (MQ=255)
                                                                                                                                    aCCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTaaaaa                                                                               >  1:274289/1‑79 (MQ=255)
                                                                                                                                    aCCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATg         >  3:242460/1‑149 (MQ=255)
                                                                                                                                    aCCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATg         >  4:556610/1‑149 (MQ=255)
                                                                                                                                    aCCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATg         >  3:34372/1‑149 (MQ=255)
                                                                                                                                       tCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTg      >  3:32222/1‑149 (MQ=255)
                                                                                                                                        cacaGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGa     <  2:319822/149‑1 (MQ=255)
                                                                                                                                        cacaGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGa     <  2:134282/149‑1 (MQ=255)
                                                                                                                                         acaGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAg    >  4:491918/1‑149 (MQ=255)
                                                                                                                                         acaGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAg    >  3:778477/1‑149 (MQ=255)
                                                                                                                                         acaGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAg    >  1:617473/1‑149 (MQ=255)
                                                                                                                                           aGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGtt  >  3:163359/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
TCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTT  >  NC_000913/872701‑872988

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