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.

GGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGG  >  NC_000913/192433‑192721
                                                                                                                                                   |                                                                                                                                              
ggAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGaataat                                                                                                                                              >  4:1100140/1‑149 (MQ=255)
  aaaaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTc                                                                                                                                            <  4:609221/149‑1 (MQ=255)
  aaaaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTc                                                                                                                                            >  1:1466686/1‑149 (MQ=255)
        agCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGa                                                                                                                                      <  2:1466686/149‑1 (MQ=255)
        agCTGAAAGTCATGGACCTGGCGCCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGa                                                                                                                                      <  3:760383/149‑1 (MQ=255)
         gCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGAt                                                                                                                                     <  2:883187/149‑1 (MQ=255)
         gCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGAt                                                                                                                                     >  3:71394/1‑149 (MQ=255)
            gAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGa                                                                                                                                  >  1:1458720/1‑149 (MQ=255)
             aaaGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTCCGTGAttaaa                                                                                                                                 <  1:98798/149‑5 (MQ=255)
                   aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataaataa                                                                                                                           <  3:1100140/149‑1 (MQ=255)
                   aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataaataa                                                                                                                           >  3:1152715/1‑149 (MQ=255)
                    tGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAg                                                                                                                          <  1:125094/149‑1 (MQ=255)
                     ggACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                         <  2:1344672/149‑1 (MQ=255)
                         cTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataa                                                                                                                               >  2:566399/1‑139 (MQ=255)
                         cTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGataa                                                                                                                               <  1:566399/139‑1 (MQ=255)
                              ggcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTc                                                                                                                >  3:1194920/1‑149 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGat                                                                                                                                 >  2:647935/1‑129 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                             <  1:493060/149‑1 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                             >  2:493060/1‑149 (MQ=255)
                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                             <  4:1194920/149‑1 (MQ=255)
                                  ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                            >  1:1304741/1‑149 (MQ=255)
                                   tCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGt                                                                                                           <  4:71394/149‑1 (MQ=255)
                                     aCGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                         <  1:807752/149‑1 (MQ=255)
                                     aCGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                         <  4:1161497/149‑1 (MQ=255)
                                             tCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGa                                                                                                                                      <  1:577232/112‑1 (MQ=255)
                                             tCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGa                                                                                                                                      >  2:577232/1‑112 (MQ=255)
                                             tCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcgcg                                                                                                 >  4:149160/1‑149 (MQ=255)
                                                 gACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggggg                                                                                             >  4:1091852/1‑149 (MQ=255)
                                                 gACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGCAAGATTCCGCGTAAGTATCGCggggg                                                                                             >  2:1277617/1‑149 (MQ=255)
                                                  aCCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGc                                                                                            <  1:460028/149‑1 (MQ=255)
                                                    cATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                          >  4:233457/1‑149 (MQ=255)
                                                     aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                         <  2:428773/133‑1 (MQ=255)
                                                     aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                         >  1:428773/1‑133 (MQ=255)
                                                      tAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcg                                                                                                              <  4:1308543/127‑1 (MQ=255)
                                                      tAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcg                                                                                                              >  3:1308543/1‑127 (MQ=255)
                                                          ttACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        <  2:1291685/145‑1 (MQ=255)
                                                          ttACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        >  1:1291685/1‑145 (MQ=255)
                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 <  3:1091852/149‑1 (MQ=255)
                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 <  2:1458720/149‑1 (MQ=255)
                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                 <  3:233457/149‑1 (MQ=255)
                                                                 ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGt                                                                                                                                        <  4:965484/90‑1 (MQ=255)
                                                                 ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGt                                                                                                                                        >  3:965484/1‑90 (MQ=255)
                                                                 ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTtata                                                                             >  4:4143/1‑149 (MQ=255)
                                                                    gtgtTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAgg                                                                          >  3:941930/1‑149 (MQ=255)
                                                                        tttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGtt                                                                      <  4:941930/149‑1 (MQ=255)
                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                            <  4:576355/142‑1 (MQ=255)
                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                            >  3:576355/1‑142 (MQ=255)
                                                                          tCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTaa                                                                                                                     >  4:1207943/1‑100 (MQ=255)
                                                                          tCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTaa                                                                                                                     <  3:1207943/100‑1 (MQ=255)
                                                                            aaTATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                         >  3:9317/1‑94 (MQ=255)
                                                                             atatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                            >  3:1450627/1‑138 (MQ=255)
                                                                              tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgtt                                                                <  3:149160/149‑1 (MQ=255)
                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcag                                                            >  2:299415/1‑149 (MQ=255)
                                                                                   aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggg                                                                                               >  2:384806/1‑113 (MQ=255)
                                                                                   aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggg                                                                                               <  1:384806/113‑1 (MQ=255)
                                                                                    cAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcgc                                                                                                  <  3:470948/109‑1 (MQ=255)
                                                                                     aaaCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAATAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAt                                                                    <  1:1109736/138‑1 (MQ=255)
                                                                                           ggTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAg                                                                                                                    >  4:619500/1‑84 (MQ=255)
                                                                                           ggTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAg                                                                                                                    <  3:619500/84‑1 (MQ=255)
                                                                                               cgctgcgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCAc                                               <  1:299415/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)
                                                                                                           tgGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                           <  3:783869/109‑1 (MQ=255)
                                                                                                           tgGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                           >  4:783869/1‑109 (MQ=255)
                                                                                                              tAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        <  3:480353/93‑1 (MQ=255)
                                                                                                              tAATGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                        >  4:480353/1‑93 (MQ=255)
                                                                                                               aaTGGGTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGCCCCGGGCAgtgctcgccat                               >  1:923398/1‑142 (MQ=255)
                                                                                                                    gTGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          <  2:1407425/149‑1 (MQ=255)
                                                                                                                     tGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcaggca                                                         >  1:1394624/1‑117 (MQ=255)
                                                                                                                     tGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcaggca                                                         <  2:1394624/117‑1 (MQ=255)
                                                                                                                     tGAAAAAGAAGGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAg                         >  3:252673/1‑149 (MQ=255)
                                                                                                                          aagaagGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTg                                                    >  1:1289937/1‑117 (MQ=255)
                                                                                                                          aagaagGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTg                                                    <  2:1289937/117‑1 (MQ=255)
                                                                                                                           agaagGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACa                   >  2:328773/1‑149 (MQ=255)
                                                                                                                            gaagGGACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaa                  >  4:137244/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:827808/1‑101 (MQ=255)
                                                                                                                                ggACTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                              <  3:834203/101‑1 (MQ=255)
                                                                                                                                  aCTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgtt                                                                <  1:296638/97‑1 (MQ=255)
                                                                                                                                  aCTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTAc            <  3:137244/149‑1 (MQ=255)
                                                                                                                                   cTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCt                                                                                    >  4:1405814/1‑76 (MQ=255)
                                                                                                                                        aTCACGGGAATTATTCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGtt                                                                      >  1:783238/1‑85 (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)
                                                                                                                                           acgggAATTATTCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATg   <  1:1065660/146‑1 (MQ=255)
                                                                                                                                           aCGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          <  4:647137/126‑1 (MQ=255)
                                                                                                                                           aCGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                          >  3:647137/1‑126 (MQ=255)
                                                                                                                                           aCGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATg   <  3:1270773/149‑1 (MQ=255)
                                                                                                                                            cGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATgg  >  3:1327468/1‑149 (MQ=255)
                                                                                                                                              ggAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCgggg                                                                                              >  2:527796/1‑56 (MQ=255)
                                                                                                                                              ggAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCgggg                                                                                              <  1:527796/56‑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 < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: