Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 192,573:1 +G 7.8% coding (719/726 nt) pyrH → uridylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913192,5731.G 7.8% 373.2 / 28.4 90coding (719/726 nt)pyrHuridylate kinase
Reads supporting (aligned to +/- strand):  ref base . (37/46);  new base G (4/3);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

GAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGA  >  NC_000913/192425‑192718
                                                                                                                                                     |                                                                                                                                                 
gAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAc                                                                                                                                                    <  4:664684/149‑1 (MQ=255)
gAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAc                                                                                                                                                    <  3:178139/149‑1 (MQ=255)
gAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAc                                                                                                                                                    <  2:295675/149‑1 (MQ=255)
gAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGAGTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGGGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAc                                                                                                                                                    <  1:55990/149‑1 (MQ=255)
    tGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGaa                                                                                                                                               <  3:26823/149‑1 (MQ=255)
      cTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGaata                                                                                                                                             >  4:603263/1‑149 (MQ=255)
        ggAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACgggaatta                                                                                                                                           >  2:90053/1‑143 (MQ=255)
         gAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAAtt                                                                                                                                          <  1:358538/149‑1 (MQ=255)
           aaaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTcc                                                                                                                                        <  1:177449/149‑1 (MQ=255)
              agagCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGt                                                                                                                                     <  3:721699/149‑1 (MQ=255)
              agagCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGt                                                                                                                                     <  4:306099/149‑1 (MQ=255)
              agagCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGt                                                                                                                                     <  1:298221/149‑1 (MQ=255)
                agCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTccgtga                                                                                                                                    <  1:90053/149‑5 (MQ=255)
                agCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGa                                                                                                                                   >  4:602889/1‑149 (MQ=255)
                 gCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGAt                                                                                                                                  <  3:215490/149‑1 (MQ=255)
                 gCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGAt                                                                                                                                  <  4:560341/149‑1 (MQ=255)
                 gCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGAt                                                                                                                                  >  4:595108/1‑149 (MQ=255)
                  cTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATg                                                                                                                                 <  3:32818/149‑1 (MQ=255)
                   tGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTccgtga                                                                                                                                    >  2:239223/1‑142 (MQ=255)
                   tGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTccgtga                                                                                                                                    <  1:239223/146‑5 (MQ=255)
                    gAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGa                                                                                                                               <  1:27154/149‑1 (MQ=255)
                    gAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGa                                                                                                                               >  2:137707/1‑149 (MQ=255)
                     aaaGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGat                                                                                                                              <  3:53662/149‑1 (MQ=255)
                           aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataaataa                                                                                                                        <  3:153998/149‑1 (MQ=255)
                           aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataaataa                                                                                                                        <  2:31535/149‑1 (MQ=255)
                           aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataaataa                                                                                                                        <  3:603263/149‑1 (MQ=255)
                            tGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAg                                                                                                                       <  3:481725/149‑1 (MQ=255)
                            tGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAg                                                                                                                       <  4:212431/149‑1 (MQ=255)
                             ggACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                      >  3:506512/1‑149 (MQ=255)
                             ggACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                      <  4:691622/149‑1 (MQ=255)
                              gACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTCTAATCACGGAATAATTCCCGTGATGGATAAATAAggg                                                                                                                     <  1:464012/149‑1 (MQ=255)
                                    gcggcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                      >  4:424583/1‑142 (MQ=255)
                                    gcggcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                      <  3:424583/142‑1 (MQ=255)
                                        ccTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcg                                                                                                           >  1:108849/1‑149 (MQ=255)
                                         cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                          >  2:229759/1‑149 (MQ=255)
                                         cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTAAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                          >  3:466368/1‑149 (MQ=255)
                                          ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                         >  2:34101/1‑149 (MQ=255)
                                          ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                         >  4:417657/1‑149 (MQ=255)
                                                      cGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCgg                                                                                             <  3:145195/149‑1 (MQ=255)
                                                       gTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggg                                                                                            <  1:94360/149‑1 (MQ=255)
                                                       gTGACCATAAATCACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggg                                                                                            <  3:701414/149‑1 (MQ=255)
                                                           ccATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCg                                                                                        >  3:76247/1‑149 (MQ=255)
                                                            cATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                       <  1:719248/149‑1 (MQ=255)
                                                            cATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                       >  3:615008/1‑149 (MQ=255)
                                                             aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTa                                                                                      >  1:581634/1‑149 (MQ=255)
                                                              tAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                     <  4:299153/149‑1 (MQ=255)
                                                                   tACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTg                                                                                >  2:511125/1‑149 (MQ=255)
                                                                   tACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTg                                                                                <  2:547657/149‑1 (MQ=255)
                                                                    aCCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               <  3:261054/149‑1 (MQ=255)
                                                                    aCCGATTCGTGTTTTCAATATGAACAAACCCGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               <  4:112670/149‑1 (MQ=255)
                                                                      cGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGtt                                                                             <  2:10178/149‑1 (MQ=255)
                                                                         ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTtata                                                                          <  4:603469/149‑1 (MQ=255)
                                                                           cGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAATGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                        <  2:453912/149‑1 (MQ=255)
                                                                           cGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                        <  3:67920/149‑1 (MQ=255)
                                                                              gtTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCg                                                                     >  1:431305/1‑149 (MQ=255)
                                                                                tttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGtt                                                                   <  1:455486/149‑1 (MQ=255)
                                                                                   cAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAtt                                                                >  3:185874/1‑149 (MQ=255)
                                                                                    aaTATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttg                                                               >  4:602070/1‑149 (MQ=255)
                                                                                     atatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgt                                                              <  3:602889/149‑1 (MQ=255)
                                                                                      tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               <  2:388115/131‑1 (MQ=255)
                                                                                      tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               >  1:388115/1‑131 (MQ=255)
                                                                                         gAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgca                                                          <  3:595108/149‑1 (MQ=255)
                                                                                          aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAtt                                                                >  2:248887/1‑142 (MQ=255)
                                                                                          aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAtt                                                                <  1:248887/142‑1 (MQ=255)
                                                                                           aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               >  2:362491/1‑126 (MQ=255)
                                                                                           aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               <  1:362491/126‑1 (MQ=255)
                                                                                           aCAAACCGGGTGCGCTGCGCCGTGAGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcagg                                                        >  2:438251/1‑149 (MQ=255)
                                                                                             aaaCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcaggca                                                      <  4:447453/149‑1 (MQ=255)
                                                                                               aCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataaataa                                                                                                                        >  3:78842/1‑81 (MQ=255)
                                                                                               aCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataaataa                                                                                                                        <  4:78842/81‑1 (MQ=255)
                                                                                                 cGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATGCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGttt                                                  <  2:374536/149‑1 (MQ=255)
                                                                                                        gctgcgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACg                                           >  3:376674/1‑149 (MQ=255)
                                                                                                        gctgcgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACg                                           >  4:752819/1‑149 (MQ=255)
                                                                                                              ccGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAgc                                     >  1:655747/1‑149 (MQ=255)
                                                                                                                       aaTGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCaa                            >  4:557385/1‑149 (MQ=255)
                                                                                                                            gTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                       >  4:193541/1‑149 (MQ=255)
                                                                                                                            gTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                       <  3:250912/149‑1 (MQ=255)
                                                                                                                            gTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                       >  4:634195/1‑149 (MQ=255)
                                                                                                                            gTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                       >  2:605502/1‑149 (MQ=255)
                                                                                                                              gAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGc                     <  2:119615/149‑1 (MQ=255)
                                                                                                                                     aagGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaa              >  1:681667/1‑149 (MQ=255)
                                                                                                                                     aagGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaa              <  2:681667/149‑1 (MQ=255)
                                                                                                                                      agGGACTTTAATCACGGTATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaaa             <  2:674203/149‑1 (MQ=255)
                                                                                                                                      agGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaaa             <  2:431305/149‑1 (MQ=255)
                                                                                                                                        ggACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGt           >  4:288210/1‑149 (MQ=255)
                                                                                                                                            tttAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGt       >  3:371206/1‑149 (MQ=255)
                                                                                                                                                aTCACGGGAATTATTCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTa                  <  2:293961/134‑1 (MQ=255)
                                                                                                                                                aTCACGGGAATTATTCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTa                  >  1:293961/1‑134 (MQ=255)
                                                                                                                                                aTCACGGGAATTATTCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg   >  2:13757/1‑149 (MQ=255)
                                                                                                                                                aTCACGGGAATTATTCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg   >  4:424874/1‑149 (MQ=255)
                                                                                                                                                aTCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg   <  1:591076/149‑1 (MQ=255)
                                                                                                                                                 tCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa  >  1:724220/1‑149 (MQ=255)
                                                                                                                                                 tCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa  >  3:270009/1‑149 (MQ=255)
                                                                                                                                                 tCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCCGAGCGTACAAAAGTACGTGAGGa  >  3:262453/1‑149 (MQ=255)
                                                                                                                                                 tCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATGCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa  >  4:439553/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGA  >  NC_000913/192425‑192718

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