Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 192,573:1 +G 5.2% 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 5.2% 611.3 / 27.5 135coding (719/726 nt)pyrHuridylate kinase
Reads supporting (aligned to +/- strand):  ref base . (61/67);  new base G (4/3);  total (65/70)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

Only 100 of 139 total aligned reads displayed.

GGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGG  >  NC_000913/192433‑192721
                                                                                                                                             |                                                                                                                                                    
ggAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGaataat                                                                                                                                              >  4:1100140/1‑149 (MQ=255)
 gAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACg                                                                                                                                                     <  4:386706/141‑1 (MQ=255)
 gAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACg                                                                                                                                                     >  3:386706/1‑141 (MQ=255)
  aaaaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTc                                                                                                                                            <  4:609221/149‑1 (MQ=255)
  aaaaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTc                                                                                                                                            >  1:1466686/1‑149 (MQ=255)
   aaaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTc                                                                                                                                            >  2:98798/1‑149 (MQ=255)
    aaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATCACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATAccc                                                                                                                                          <  1:221999/149‑1 (MQ=255)
        agCTGAAAGTCATGGACCTGGCGCCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGa                                                                                                                                      <  3:760383/149‑1 (MQ=255)
             aaaGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTCCGTGAttaaa                                                                                                                                 <  1:98798/149‑5 (MQ=255)
                   aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataaataa                                                                                                                           >  3:1152715/1‑149 (MQ=255)
                   aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataaataa                                                                                                                           <  3:1100140/149‑1 (MQ=255)
                    tGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAg                                                                                                                          <  1:125094/149‑1 (MQ=255)
                     ggACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                         <  2:1344672/149‑1 (MQ=255)
                         cTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataa                                                                                                                               <  1:566399/139‑1 (MQ=255)
                          tGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTc                                                                                                                                            >  2:448939/1‑125 (MQ=255)
                          tGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTc                                                                                                                                            <  1:448939/125‑1 (MQ=255)
                           ggcggcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTccgtga                                                                                                                                       <  3:156810/130‑5 (MQ=255)
                           ggcggcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTccgtga                                                                                                                                       >  4:156810/1‑126 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGat                                                                                                                                 >  2:647935/1‑129 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGat                                                                                                                                 <  1:647935/129‑1 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                             >  2:493060/1‑149 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                             <  4:1194920/149‑1 (MQ=255)
                                  ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGat                                                                                                                                 <  3:695705/128‑1 (MQ=255)
                                  ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                            >  1:1304741/1‑149 (MQ=255)
                                  ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                            >  3:1161497/1‑149 (MQ=255)
                                   tCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGt                                                                                                           <  4:71394/149‑1 (MQ=255)
                                     aCGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                         <  1:807752/149‑1 (MQ=255)
                                          ggCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATc                                                                                                    <  1:1294514/149‑1 (MQ=255)
                                             tCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGa                                                                                                                                      >  2:577232/1‑112 (MQ=255)
                                                 gACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggggg                                                                                             >  4:1091852/1‑149 (MQ=255)
                                                 gACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGCAAGATTCCGCGTAAGTATCGCggggg                                                                                             >  2:1277617/1‑149 (MQ=255)
                                                  aCCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGc                                                                                            <  1:460028/149‑1 (MQ=255)
                                                     aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                         >  1:428773/1‑133 (MQ=255)
                                                     aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGAGTCCGCGTAAGTATCGCGGGGGCGTa                                                                                         >  1:1407425/1‑149 (MQ=255)
                                                      tAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcg                                                                                                              >  3:1308543/1‑127 (MQ=255)
                                                      tAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcg                                                                                                              <  4:1308543/127‑1 (MQ=255)
                                                          ttACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        >  1:1291685/1‑145 (MQ=255)
                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTggg                                                                                 <  1:1277617/149‑2 (MQ=255)
                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 <  2:1458720/149‑1 (MQ=255)
                                                                 ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGt                                                                                                                                        >  3:965484/1‑90 (MQ=255)
                                                                 ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGt                                                                                                                                        <  4:965484/90‑1 (MQ=255)
                                                                 ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTtata                                                                             >  4:4143/1‑149 (MQ=255)
                                                                    gtgtTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAgg                                                                          >  3:941930/1‑149 (MQ=255)
                                                                      gtTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCg                                                                        <  2:1304741/149‑1 (MQ=255)
                                                                        tttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGtt                                                                      <  4:941930/149‑1 (MQ=255)
                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTcc                                                                                                               >  4:1352899/1‑107 (MQ=255)
                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTcc                                                                                                               <  3:1352899/107‑1 (MQ=255)
                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                            <  4:576355/142‑1 (MQ=255)
                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                            >  3:576355/1‑142 (MQ=255)
                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTa                                                                     <  4:1152715/149‑1 (MQ=255)
                                                                          tCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTaa                                                                                                                     <  3:1207943/100‑1 (MQ=255)
                                                                          tCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTaa                                                                                                                     >  4:1207943/1‑100 (MQ=255)
                                                                            aaTATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttg                                                                  >  1:55240/1‑149 (MQ=255)
                                                                             atatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                            <  4:1450627/138‑1 (MQ=255)
                                                                              tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgtt                                                                <  3:149160/149‑1 (MQ=255)
                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcag                                                            >  2:299415/1‑149 (MQ=255)
                                                                                   aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggg                                                                                               >  2:384806/1‑113 (MQ=255)
                                                                                    cAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcgc                                                                                                  <  3:470948/109‑1 (MQ=255)
                                                                                     aaaCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAATAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAt                                                                    <  1:1109736/138‑1 (MQ=255)
                                                                                           ggTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAg                                                                                                                    <  3:619500/84‑1 (MQ=255)
                                                                                           ggTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAg                                                                                                                    >  4:619500/1‑84 (MQ=255)
                                                                                           ggTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTgg                                                   <  2:55240/149‑1 (MQ=255)
                                                                                                 ctgcgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACgg                                             <  1:1401687/149‑1 (MQ=255)
                                                                                                      ccGTGTGGTAATGGGTGAAAAAGAAGGGACTTTGATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 >  3:1378237/1‑108 (MQ=255)
                                                                                                      ccGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 <  4:1378237/108‑1 (MQ=255)
                                                                                                          gtgGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGt                                                                                      >  1:309851/1‑99 (MQ=255)
                                                                                                           tgGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                           >  4:783869/1‑109 (MQ=255)
                                                                                                           tgGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                           <  3:783869/109‑1 (MQ=255)
                                                                                                              tAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        <  3:480353/93‑1 (MQ=255)
                                                                                                              tAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        >  4:480353/1‑93 (MQ=255)
                                                                                                                   ggTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTa                                                                     <  3:1212336/107‑1 (MQ=255)
                                                                                                                    gTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          <  2:1407425/149‑1 (MQ=255)
                                                                                                                     tGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcaggca                                                         <  2:1394624/117‑1 (MQ=255)
                                                                                                                     tGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAg                         >  3:252673/1‑149 (MQ=255)
                                                                                                                          aagaagGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTg                                                    <  2:1289937/117‑1 (MQ=255)
                                                                                                                          aagaagGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTAc                    >  3:846708/1‑149 (MQ=255)
                                                                                                                           agaagGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACa                   >  2:328773/1‑149 (MQ=255)
                                                                                                                              agGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaaa                <  4:846708/149‑1 (MQ=255)
                                                                                                                               gggACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGATTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAg               >  4:1270773/1‑149 (MQ=255)
                                                                                                                                ggACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                              >  4:834203/1‑101 (MQ=255)
                                                                                                                                ggACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                              <  4:827808/101‑1 (MQ=255)
                                                                                                                                ggACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                              <  3:834203/101‑1 (MQ=255)
                                                                                                                                ggACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                              >  3:827808/1‑101 (MQ=255)
                                                                                                                                  aCTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgtt                                                                <  1:296638/97‑1 (MQ=255)
                                                                                                                                  aCTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgtt                                                                >  2:296638/1‑97 (MQ=255)
                                                                                                                                   cTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCt                                                                                    >  4:1405814/1‑76 (MQ=255)
                                                                                                                                   cTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCt                                                                                    <  3:1405814/76‑1 (MQ=255)
                                                                                                                                        aTCACGGGAATTATTCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGtt                                                                      >  1:783238/1‑85 (MQ=255)
                                                                                                                                        aTCACGGGAATTATTCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGtt                                                                      <  2:783238/85‑1 (MQ=255)
                                                                                                                                        aTCACGGGAATTATTCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg      >  2:1065660/1‑149 (MQ=255)
                                                                                                                                         tCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGt                                                  >  1:3091/1‑104 (MQ=255)
                                                                                                                                         tCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGt                                                  <  2:3091/104‑1 (MQ=255)
                                                                                                                                         tCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          <  4:150583/128‑1 (MQ=255)
                                                                                                                                         tCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          >  3:150583/1‑128 (MQ=255)
                                                                                                                                         tCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa     >  2:942355/1‑149 (MQ=255)
                                                                                                                                          cACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGAt    <  4:252673/149‑1 (MQ=255)
                                                                                                                                           aCGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          <  4:647137/126‑1 (MQ=255)
                                                                                                                                           aCGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          >  3:647137/1‑126 (MQ=255)
                                                                                                                                            cGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATgg  >  3:1327468/1‑149 (MQ=255)
                                                                                                                                            cGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATgg  <  1:942355/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGG  >  NC_000913/192433‑192721

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