Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,202 T→A 7.9% intergenic (+34/‑193) yagF → / → yagG CP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913285,2020TA 7.9% 265.3 / 7.5 87intergenic (+34/‑193)yagF/yagGCP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base T (41/39);  new base A (2/5);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCAT  >  NC_000913/285054‑285348
                                                                                                                                                    |                                                                                                                                                  
ccGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATAtt                                                                                                                                                    >  7:94609/1‑149 (MQ=255)
          gcggcACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGAc                                                                                                                                          <  1:369562/149‑1 (MQ=255)
            ggcACTGCAGTCGGTGAGCGGCGGCACCTGGACAGGCTGTATTTATGACACCGATAAACTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGa                                                                                                                                           <  1:25671/146‑1 (MQ=255)
            ggcACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAACTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGa                                                                                                                                           >  2:25671/1‑146 (MQ=255)
            ggcACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACgg                                                                                                                                        >  6:285443/1‑149 (MQ=255)
            ggcACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACgg                                                                                                                                        <  4:87509/149‑1 (MQ=255)
                   cAGTCGTTGAGCGGAGTCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACGGttttatc                                                                                                                                 <  2:66460/149‑4 (MQ=38)
                   cAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTACCGCAGACTGCtttttt                                                                                                                                 >  8:285833/1‑149 (MQ=255)
                   cAGTCGGTGAGAGGCGGCACCTCGAAAGGCTGTATTTATGACACAGATAAAATTATAGAGGTAATTAAAGCAGGTAAAAAAGCGCTCGTAATTTAATTATTTTAAGAGATAAAACAGTCTGCGGAATATTTCCCGCAGACTGCAttgtt                                                                                                                                 >  6:189562/1‑149 (MQ=255)
                         gTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGttttt                                                                                                                           >  7:289064/1‑149 (MQ=255)
                         gTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGttttt                                                                                                                           <  2:344381/149‑1 (MQ=255)
                          tGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGAc                                                                                                                                          >  1:131817/1‑133 (MQ=255)
                          tGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGAc                                                                                                                                          <  2:131817/133‑1 (MQ=255)
                          tGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTg                                                                                                                          <  7:377180/149‑1 (MQ=255)
                           gAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTg                                                                                                                          >  4:237376/1‑148 (MQ=255)
                                gcggcACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAttt                                                                                                                    <  8:84736/149‑1 (MQ=255)
                                gcggcACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAtt                                                                                                                     >  2:170670/1‑148 (MQ=255)
                                        tGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAAtt                                                                                                             >  3:220076/1‑148 (MQ=255)
                                        tGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAAtt                                                                                                             >  5:121023/1‑148 (MQ=255)
                                           aaaGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGTCTTTGTTGTTTTTGAAATTTATTAATTTaaa                                                                                                         <  5:396513/149‑1 (MQ=255)
                                           aaaGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTaaa                                                                                                         <  4:127858/149‑1 (MQ=255)
                                            aaGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTaaa                                                                                                         <  5:266169/148‑1 (MQ=255)
                                               gCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAAt                                                                                                     >  1:60728/1‑149 (MQ=255)
                                                   tATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTa                                                                                                   <  5:285443/147‑1 (MQ=255)
                                                   tATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGt                                                                                                 <  3:237376/149‑1 (MQ=255)
                                                   tATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGt                                                                                                 <  5:255720/149‑1 (MQ=255)
                                                     tttATGACACCGATAAAATTTTCGAGGTTATTTACGCCGGTATAATAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATAATTCCCGCAGTCGGCTTTGTTGTTTTTGATATTTATTAATTTAAACAATTAGtt                                                                                                <  7:358082/148‑1 (MQ=255)
                                                       tATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTga                                                                                              <  3:6913/148‑1 (MQ=255)
                                                        aTGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGACTATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTgag                                                                                             <  7:285833/148‑1 (MQ=255)
                                                         tGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTgaga                                                                                            >  7:396204/1‑148 (MQ=255)
                                                         tGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTgaga                                                                                            <  1:170670/148‑1 (MQ=255)
                                                         tGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTgaga                                                                                            <  8:396204/148‑1 (MQ=255)
                                                            cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatata                                                                                        >  6:83788/1‑149 (MQ=255)
                                                             acCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatatat                                                                                       >  5:400301/1‑149 (MQ=255)
                                                             acCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTACCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatatat                                                                                       >  4:63760/1‑149 (MQ=255)
                                                              ccGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATc                                                                                      >  6:27082/1‑149 (MQ=255)
                                                              ccGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATc                                                                                      >  8:394981/1‑149 (MQ=255)
                                                               cGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCg                                                                                     >  3:333011/1‑149 (MQ=255)
                                                               cGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCGTTGTTGTTTTTGATATTTATTATGTTAAACGAGTAGTTGCGatatatgg                                                                                     >  2:264984/1‑147 (MQ=255)
                                                                  tAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTg                                                                                  <  8:216881/149‑1 (MQ=255)
                                                                  tAAAATTATAGCCGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGGCCGCGGAATAATTCCCGCAGTAGGCTTTGTTGGTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTg                                                                                  <  8:300658/149‑1 (MQ=255)
                                                                     aaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTa                                                                                                   >  2:325912/1‑129 (MQ=255)
                                                                     aaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTa                                                                                                   <  1:325912/129‑1 (MQ=255)
                                                                          tCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTcac                                                                           <  7:210614/148‑1 (MQ=255)
                                                                          tCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTcac                                                                           <  6:122364/148‑1 (MQ=255)
                                                                            gAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaaa                                                                        <  5:83788/149‑1 (MQ=255)
                                                                            gAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaa                                                                         >  7:386089/1‑148 (MQ=255)
                                                                             aGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaaaa                                                                       <  6:168391/149‑1 (MQ=255)
                                                                             aGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaaaa                                                                       <  2:60728/149‑1 (MQ=255)
                                                                             aGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaaaa                                                                       <  4:220076/149‑1 (MQ=255)
                                                                             aGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaa                                                                         >  7:125418/1‑147 (MQ=255)
                                                                                  aTTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCaa                                                                   >  7:164853/1‑148 (MQ=255)
                                                                                  aTTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCaa                                                                   >  3:252724/1‑148 (MQ=255)
                                                                                      aCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTaaa                                                                                                         <  8:406688/106‑1 (MQ=255)
                                                                                      aCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTaaa                                                                                                         >  7:406688/1‑106 (MQ=255)
                                                                                      aCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATaa                                                              <  8:164853/149‑1 (MQ=255)
                                                                                           ggTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTTGCGTCACAAAAGCAAAATAACGTaa                                                         >  7:255454/1‑149 (MQ=255)
                                                                                           ggTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTaa                                                         >  5:321558/1‑149 (MQ=255)
                                                                                            gTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTaa                                                         <  2:162146/148‑1 (MQ=255)
                                                                                            gTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTaa                                                         >  1:162146/1‑148 (MQ=255)
                                                                                             tAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGTAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAAtt                                                       <  8:125418/149‑1 (MQ=255)
                                                                                              aaaaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTc                                                      >  8:309111/1‑149 (MQ=255)
                                                                                                aaaaGCGCTCGGAATTTTATTATTTTAAGAGATAAAACCGTCTGCGGATTATTTCCCGCAGACGCCTTTGTTGTTTTTGAAATTAATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGAAAAATAACGTAATTCg                                                     >  8:182305/1‑148 (MQ=255)
                                                                                                aaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAAGTAGTTGAGATATATCGTTGGCGTCACAAAAGCAACATAACGTAATTCg                                                     >  6:225787/1‑148 (MQ=255)
                                                                                                  aaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCAGAAAAGCAAAATAACTTAATTCGGa                                                   >  1:58571/1‑148 (MQ=255)
                                                                                                  aaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACTTAATTCGGa                                                   >  1:108683/1‑148 (MQ=255)
                                                                                                         cGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGat                                           <  6:400301/149‑1 (MQ=255)
                                                                                                             aTTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGa                                       >  6:387545/1‑149 (MQ=255)
                                                                                                               ttAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAATAGAAAAATAACGTAATTCGGAAATAAGATATGAcc                                     >  7:263422/1‑149 (MQ=255)
                                                                                                               ttAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGAc                                      >  5:233603/1‑148 (MQ=255)
                                                                                                                 aaTTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCAt                                   >  2:337970/1‑149 (MQ=255)
                                                                                                                        ttAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGATTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTTATATATATAGTTTGCGTCACAAAAGCAAAATAACTTAATTCGGAAATAATATATGACCATTGCTGGt                            >  2:281634/1‑149 (MQ=255)
                                                                                                                            gagaTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGtt                                                                                   <  6:90759/90‑1 (MQ=255)
                                                                                                                            gagaTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGtt                                                                                   >  5:90759/1‑90 (MQ=255)
                                                                                                                             agaTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAAtt                       <  5:387545/149‑1 (MQ=255)
                                                                                                                                aaaaaCCGTATGCGTAATATATCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAATATCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTTGTTAATTGaa                    <  5:189562/148‑1 (MQ=255)
                                                                                                                                  aaaCCTTATTCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGTCGTCACAAAAGCAAAATAACGTTATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAAt                   <  2:139353/148‑1 (MQ=255)
                                                                                                                                  aaaCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATa                  <  4:333011/149‑1 (MQ=255)
                                                                                                                                  aaaCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATa                  <  8:386089/149‑1 (MQ=255)
                                                                                                                                   aaCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAg                 <  1:264984/149‑1 (MQ=255)
                                                                                                                                    aCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGc                <  7:309111/149‑1 (MQ=255)
                                                                                                                                      cGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCt               >  4:216675/1‑148 (MQ=255)
                                                                                                                                       gTCTGCTGAATATTTCCCTCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCa             <  7:411901/149‑1 (MQ=255)
                                                                                                                                       gTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCa             >  6:400192/1‑149 (MQ=255)
                                                                                                                                             ggAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAATAGCAAAATAACGTAATTCGGAAATAAGATATGAGCATTGCTGGTTAATTGAATAGCTCATTacac       >  5:399104/1‑149 (MQ=255)
                                                                                                                                              gAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCATAATAAAGAAATTCCGAAATAAGATATGACCATTCCTGGTTAATTGAATAGATCATTACACt      >  6:23609/1‑149 (MQ=255)
                                                                                                                                                 tatTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTTCa   <  8:289064/149‑1 (MQ=255)
                                                                                                                                                  atTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGGTAATTGAATAGCTCCTTACACTCCAt  <  8:255454/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCAT  >  NC_000913/285054‑285348

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