Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 192,579 A→T 5.8% *242L (TAA→TTA)  pyrH → uridylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913192,5790AT 5.8% 497.6 / 15.4 137*242L (TAA→TTA) pyrHuridylate kinase
Reads supporting (aligned to +/- strand):  ref base A (61/68);  new base T (4/4);  total (65/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 140 total aligned reads displayed.

GGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGG  >  NC_000913/192433‑192721
                                                                                                                                                   |                                                                                                                                              
ggAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGaataat                                                                                                                                              >  4:1100140/1‑149 (MQ=255)
  aaaaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTc                                                                                                                                            <  4:609221/149‑1 (MQ=255)
  aaaaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTc                                                                                                                                            >  1:1466686/1‑149 (MQ=255)
        agCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGa                                                                                                                                      <  2:1466686/149‑1 (MQ=255)
        agCTGAAAGTCATGGACCTGGCGCCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGa                                                                                                                                      <  3:760383/149‑1 (MQ=255)
         gCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGAt                                                                                                                                     <  2:883187/149‑1 (MQ=255)
         gCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGAt                                                                                                                                     >  3:71394/1‑149 (MQ=255)
            gAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGa                                                                                                                                  >  1:1458720/1‑149 (MQ=255)
             aaaGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATT‑CCGTGAttaaa                                                                                                                                 <  1:98798/149‑5 (MQ=255)
                   aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataaataa                                                                                                                           <  3:1100140/149‑1 (MQ=255)
                   aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataaataa                                                                                                                           >  3:1152715/1‑149 (MQ=255)
                    tGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAg                                                                                                                          <  1:125094/149‑1 (MQ=255)
                     ggACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                         <  2:1344672/149‑1 (MQ=255)
                         cTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataa                                                                                                                               >  2:566399/1‑139 (MQ=255)
                         cTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataa                                                                                                                               <  1:566399/139‑1 (MQ=255)
                              ggcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTc                                                                                                                >  3:1194920/1‑149 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGat                                                                                                                                 >  2:647935/1‑129 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                             <  1:493060/149‑1 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                             >  2:493060/1‑149 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                             <  4:1194920/149‑1 (MQ=255)
                                  ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                            >  1:1304741/1‑149 (MQ=255)
                                   tCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGt                                                                                                           <  4:71394/149‑1 (MQ=255)
                                     aCGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                         <  1:807752/149‑1 (MQ=255)
                                     aCGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                         <  4:1161497/149‑1 (MQ=255)
                                             tCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGa                                                                                                                                      <  1:577232/112‑1 (MQ=255)
                                             tCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGa                                                                                                                                      >  2:577232/1‑112 (MQ=255)
                                             tCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcgcg                                                                                                 >  4:149160/1‑149 (MQ=255)
                                                 gACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggggg                                                                                             >  4:1091852/1‑149 (MQ=255)
                                                 gACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGCAAGATTCCGCGTAAGTATCGCggggg                                                                                             >  2:1277617/1‑149 (MQ=255)
                                                  aCCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGc                                                                                            <  1:460028/149‑1 (MQ=255)
                                                    cATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                          >  4:233457/1‑149 (MQ=255)
                                                     aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                         <  2:428773/133‑1 (MQ=255)
                                                     aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                         >  1:428773/1‑133 (MQ=255)
                                                      tAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcg                                                                                                              <  4:1308543/127‑1 (MQ=255)
                                                      tAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcg                                                                                                              >  3:1308543/1‑127 (MQ=255)
                                                          ttACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        <  2:1291685/145‑1 (MQ=255)
                                                          ttACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        >  1:1291685/1‑145 (MQ=255)
                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 <  3:1091852/149‑1 (MQ=255)
                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 <  2:1458720/149‑1 (MQ=255)
                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 <  3:233457/149‑1 (MQ=255)
                                                                 ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGt                                                                                                                                        <  4:965484/90‑1 (MQ=255)
                                                                 ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGt                                                                                                                                        >  3:965484/1‑90 (MQ=255)
                                                                 ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTtata                                                                             >  4:4143/1‑149 (MQ=255)
                                                                    gtgtTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAgg                                                                          >  3:941930/1‑149 (MQ=255)
                                                                        tttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGtt                                                                      <  4:941930/149‑1 (MQ=255)
                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                            <  4:576355/142‑1 (MQ=255)
                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                            >  3:576355/1‑142 (MQ=255)
                                                                          tCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTaa                                                                                                                     >  4:1207943/1‑100 (MQ=255)
                                                                          tCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTaa                                                                                                                     <  3:1207943/100‑1 (MQ=255)
                                                                            aaTATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                         >  3:9317/1‑94 (MQ=255)
                                                                             atatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                            >  3:1450627/1‑138 (MQ=255)
                                                                              tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgtt                                                                <  3:149160/149‑1 (MQ=255)
                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcag                                                            >  2:299415/1‑149 (MQ=255)
                                                                                   aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggg                                                                                               >  2:384806/1‑113 (MQ=255)
                                                                                   aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggg                                                                                               <  1:384806/113‑1 (MQ=255)
                                                                                    cAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcgc                                                                                                  <  3:470948/109‑1 (MQ=255)
                                                                                     aaaCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAATAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAt                                                                    <  1:1109736/138‑1 (MQ=255)
                                                                                           ggTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAg                                                                                                                    >  4:619500/1‑84 (MQ=255)
                                                                                           ggTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAg                                                                                                                    <  3:619500/84‑1 (MQ=255)
                                                                                               cgctgcgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCAc                                               <  1:299415/149‑1 (MQ=255)
                                                                                                 ctgcgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACgg                                             <  1:1401687/149‑1 (MQ=255)
                                                                                                      ccGTGTGGTAATGGGTGAAAAAGAAGGGACTTTGATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 >  3:1378237/1‑108 (MQ=255)
                                                                                                      ccGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 <  4:1378237/108‑1 (MQ=255)
                                                                                                           tgGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                           <  3:783869/109‑1 (MQ=255)
                                                                                                           tgGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                           >  4:783869/1‑109 (MQ=255)
                                                                                                              tAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        <  3:480353/93‑1 (MQ=255)
                                                                                                              tAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        >  4:480353/1‑93 (MQ=255)
                                                                                                               aaTGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGCCCCGGGCAgtgctcgccat                               >  1:923398/1‑142 (MQ=255)
                                                                                                                    gTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          <  2:1407425/149‑1 (MQ=255)
                                                                                                                     tGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcaggca                                                         >  1:1394624/1‑117 (MQ=255)
                                                                                                                     tGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcaggca                                                         <  2:1394624/117‑1 (MQ=255)
                                                                                                                     tGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAg                         >  3:252673/1‑149 (MQ=255)
                                                                                                                          aagaagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTg                                                    >  1:1289937/1‑117 (MQ=255)
                                                                                                                          aagaagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTg                                                    <  2:1289937/117‑1 (MQ=255)
                                                                                                                           agaagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACa                   >  2:328773/1‑149 (MQ=255)
                                                                                                                            gaagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaa                  >  4:137244/1‑149 (MQ=255)
                                                                                                                              agGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaaa                <  4:846708/149‑1 (MQ=255)
                                                                                                                               gggACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGATTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAg               >  4:1270773/1‑149 (MQ=255)
                                                                                                                                ggACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                              >  4:834203/1‑101 (MQ=255)
                                                                                                                                ggACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                              <  4:827808/101‑1 (MQ=255)
                                                                                                                                ggACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                              >  3:827808/1‑101 (MQ=255)
                                                                                                                                ggACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                              <  3:834203/101‑1 (MQ=255)
                                                                                                                                  aCTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgtt                                                                <  1:296638/97‑1 (MQ=255)
                                                                                                                                  aCTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTAc            <  3:137244/149‑1 (MQ=255)
                                                                                                                                   cTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCt                                                                                    >  4:1405814/1‑76 (MQ=255)
                                                                                                                                        aTCACGGGAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGtt                                                                      >  1:783238/1‑85 (MQ=255)
                                                                                                                                        aTCACGGGAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg      >  2:1065660/1‑149 (MQ=255)
                                                                                                                                         tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGt                                                  >  1:3091/1‑104 (MQ=255)
                                                                                                                                         tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGt                                                  <  2:3091/104‑1 (MQ=255)
                                                                                                                                         tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          <  4:150583/128‑1 (MQ=255)
                                                                                                                                         tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          >  3:150583/1‑128 (MQ=255)
                                                                                                                                         tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa     >  2:942355/1‑149 (MQ=255)
                                                                                                                                          cAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGAt    <  4:252673/149‑1 (MQ=255)
                                                                                                                                           acgggAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATg   <  1:1065660/146‑1 (MQ=255)
                                                                                                                                           aC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          <  4:647137/126‑1 (MQ=255)
                                                                                                                                           aC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          >  3:647137/1‑126 (MQ=255)
                                                                                                                                           aC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATg   <  3:1270773/149‑1 (MQ=255)
                                                                                                                                            c‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATgg  >  3:1327468/1‑149 (MQ=255)
                                                                                                                                              ggAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCgggg                                                                                              >  2:527796/1‑56 (MQ=255)
                                                                                                                                              ggAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCgggg                                                                                              <  1:527796/56‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGG  >  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 < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑