Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,005,312 G→A 5.2% R182H (CGC→CAC)  pyrD → dihydro‑orotate oxidase, FMN‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,005,3120GA 5.2% 259.7 / 3.4 77R182H (CGC→CAC) pyrDdihydro‑orotate oxidase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base A (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

CGTCCTGGGTATTAACATCGGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGC  >  NC_000913/1005169‑1005459
                                                                                                                                               |                                                                                                                                                   
cGTCCTGGGTATTAACATCGGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTacgcacgc                                                                                                                                                <  1:587585/149‑1 (MQ=255)
  tCCTGGGTATTAACATCGGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTg                                                                                                                                              >  2:791696/1‑149 (MQ=255)
    cTGGGTATTAACATCGGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCa                                                                                                                                            >  3:575901/1‑149 (MQ=255)
      gggTATTAACATCGGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAat                                                                                                                                          >  3:178178/1‑149 (MQ=255)
      gggTATTAACATCGGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAat                                                                                                                                          >  3:48503/1‑149 (MQ=255)
       ggTATTAACATCGGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAata                                                                                                                                         <  2:481458/149‑1 (MQ=255)
         tATTAACATCGGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATg                                                                                                                                       >  4:797715/1‑149 (MQ=255)
          aTTAACATCAGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATgg                                                                                                                                      >  3:254615/1‑149 (MQ=255)
            tAACATCGGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTg                                                                                                                                    >  1:448495/1‑149 (MQ=255)
                  ctgCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAgcgc                                                                                                                              <  2:736101/147‑1 (MQ=255)
                   ggCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCt                                                                                                                             <  3:187246/149‑1 (MQ=255)
                    gCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGt                                                                                                                                     <  4:424665/140‑1 (MQ=255)
                    gCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGt                                                                                                                                     >  3:424665/1‑140 (MQ=255)
                       aaaaTAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGgat                                                                                                                         <  1:150787/149‑1 (MQ=255)
                           tAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCGGGATTACGCACGCTGCAATATTGTGAAGCGCTGGATGAtc                                                                                                                     <  1:235949/149‑1 (MQ=255)
                           tAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGAtc                                                                                                                     <  2:164284/149‑1 (MQ=255)
                            aaaGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGAtct                                                                                                                    <  3:690017/149‑1 (MQ=255)
                               gATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCtt                                                                                                                 >  1:376956/1‑149 (MQ=255)
                               gATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCtt                                                                                                                 <  4:254615/149‑1 (MQ=255)
                                    gCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACcg                                                                                                            >  4:128190/1‑149 (MQ=255)
                                     ccAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCt                                                                                                                                               >  3:604123/1‑113 (MQ=255)
                                     ccAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCt                                                                                                                                               <  4:604123/113‑1 (MQ=255)
                                       aGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGa                                                                                                         <  4:624681/149‑1 (MQ=255)
                                             gCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTaaaa                                                                                                   >  1:727957/1‑149 (MQ=255)
                                               aGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTacgcac                                                                                                                                                  <  1:663244/100‑1 (MQ=255)
                                               aGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTacgcac                                                                                                                                                  >  2:663244/1‑100 (MQ=255)
                                                 ggCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATaa                                                                                               <  2:695424/149‑1 (MQ=255)
                                                  gCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAg                                                                                              >  3:240539/1‑149 (MQ=255)
                                                         tGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATg                                                                                       <  4:48503/149‑1 (MQ=255)
                                                         tGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATg                                                                                       <  2:727957/149‑1 (MQ=255)
                                                         tGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATg                                                                                       <  4:178178/149‑1 (MQ=255)
                                                            cTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGAtt                                                                                    <  2:448495/149‑1 (MQ=255)
                                                               tCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGc                                                                                 <  3:285844/149‑1 (MQ=255)
                                                                  gATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAg                                                                              <  4:151328/149‑1 (MQ=255)
                                                                            gAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCa                                                                    <  4:575901/149‑1 (MQ=255)
                                                                            gAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACACACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGa                                                                                      >  2:410600/1‑131 (MQ=255)
                                                                            gAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACACACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGa                                                                                      >  2:65520/1‑131 (MQ=255)
                                                                            gAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACACACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGa                                                                                      <  1:410600/131‑1 (MQ=255)
                                                                            gAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACACACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGa                                                                                      <  1:65520/131‑1 (MQ=255)
                                                                             aaaaaaTCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCAt                                                                   >  1:550713/1‑149 (MQ=255)
                                                                                      atgcctatgcGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATAtgtg                                                          <  1:123155/149‑1 (MQ=255)
                                                                                           tatgcGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGAt                                                     <  2:642009/149‑1 (MQ=255)
                                                                                           tatgcGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGAt                                                     >  1:642009/1‑149 (MQ=255)
                                                                                               cGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCa                                                 <  1:416812/149‑1 (MQ=255)
                                                                                               cGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCa                                                 <  3:178675/149‑1 (MQ=255)
                                                                                                    tatatCGCCATGAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGaa                                            >  3:516719/1‑149 (MQ=255)
                                                                                                    tatatCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGaa                                            <  4:276263/149‑1 (MQ=255)
                                                                                                       atCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGAt                                         >  4:773109/1‑149 (MQ=255)
                                                                                                       atCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGAt                                         >  2:13969/1‑149 (MQ=255)
                                                                                                            cATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCtt                                                                                                                 <  3:346958/72‑1 (MQ=255)
                                                                                                            cATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCtt                                                                                                                 >  4:346958/1‑72 (MQ=255)
                                                                                                             aTCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGcc                                   >  3:371951/1‑149 (MQ=255)
                                                                                                              tCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCg                                  >  3:338070/1‑149 (MQ=255)
                                                                                                               cAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCgg                                 <  4:371951/149‑1 (MQ=255)
                                                                                                                aaTATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGa                                >  4:245787/1‑149 (MQ=255)
                                                                                                                 atatTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGAt                               >  3:239424/1‑149 (MQ=255)
                                                                                                                   atTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCt                             <  1:584977/149‑1 (MQ=255)
                                                                                                                       cATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCt                         >  4:191498/1‑149 (MQ=255)
                                                                                                                        aTCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTg                        >  1:538421/1‑149 (MQ=255)
                                                                                                                          cGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaa                      >  4:653434/1‑149 (MQ=255)
                                                                                                                          cGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaa                      >  4:604414/1‑149 (MQ=255)
                                                                                                                          cGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaa                      >  4:410328/1‑149 (MQ=255)
                                                                                                                            ccGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaaga                    >  2:182437/1‑149 (MQ=255)
                                                                                                                                    cccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaa                      <  3:633315/139‑1 (MQ=255)
                                                                                                                                    cccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaa                      >  4:633315/1‑139 (MQ=255)
                                                                                                                                    cccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGa            >  1:721798/1‑149 (MQ=255)
                                                                                                                                    cccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGa            <  3:245787/149‑1 (MQ=255)
                                                                                                                                    cccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGa            >  1:269796/1‑149 (MQ=255)
                                                                                                                                    cccAGGATTACGCACGCTGCAATAAGGTTAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGa            >  4:595750/1‑149 (MQ=255)
                                                                                                                                     ccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGcc                                   <  4:426506/125‑1 (MQ=255)
                                                                                                                                     ccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGcc                                   >  3:426506/1‑125 (MQ=255)
                                                                                                                                     ccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaagaa                >  2:103715/1‑144 (MQ=255)
                                                                                                                                     ccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaagaa                <  1:103715/144‑1 (MQ=255)
                                                                                                                                     ccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGAt           >  1:328342/1‑149 (MQ=255)
                                                                                                                                     ccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGAt           >  3:476833/1‑149 (MQ=255)
                                                                                                                                     ccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGAt           >  2:548311/1‑149 (MQ=255)
                                                                                                                                     ccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGAt           <  4:239424/149‑1 (MQ=255)
                                                                                                                                         gATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAata                                                              >  4:391833/1‑94 (MQ=255)
                                                                                                                                         gATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAata                                                              <  3:391833/94‑1 (MQ=255)
                                                                                                                                          aTTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAgg      <  4:476833/149‑1 (MQ=255)
                                                                                                                                              cgcacgcTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGc  <  1:182437/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CGTCCTGGGTATTAACATCGGCAAAAATAAAGATACGCCAGTGGAGCAGGGCAAAGATGACTATCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGC  >  NC_000913/1005169‑1005459

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