Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,197,756 T→A 5.6% L14Q (CTG→CAG)  hemE → uroporphyrinogen decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,197,7560TA 5.6% 328.6 / 8.8 120L14Q (CTG→CAG) hemEuroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (57/53);  new base A (2/5);  total (59/61)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

Only 100 of 126 total aligned reads displayed.

GCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTG  >  NC_000913/4197615‑4197898
                                                                                                                                             |                                                                                                                                              
gCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGctgctg                                                                                                                                               >  2:176596/1‑143 (MQ=255)
    ggtcgcgcgcCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGcc                                                                                                                                           <  1:300476/139‑1 (MQ=255)
      tAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAg                                                                                                                                         >  2:46877/1‑143 (MQ=255)
        gcgcgcCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGcc                                                                                                                                       <  2:259370/143‑1 (MQ=255)
        gcgcgcCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCAGCAGCGCCCGca                                                                                                                                       <  1:606742/143‑2 (MQ=255)
        gcgcgcCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCAGCAGCGCCCGca                                                                                                                                       <  2:1150762/143‑2 (MQ=255)
          gcgcCGTCTGTTAGGAGATACGGTGGCTATCTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAGCAGCCAAGATGACCGAACTTAAAAACGATCGTTATCTGCGTGCGCAGCAGCGcccgcaga                                                                                                                                     <  2:250404/143‑7 (MQ=255)
          gcgcCGTCTGATAGGAGATACGGTGGCGAGGTGGCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAGGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCAGCAGCGCCCGcaga                                                                                                                                     <  2:610603/143‑4 (MQ=255)
           cgcCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTACGGAACAGCCAAAATGACCGCACTTAAAAACGATCGTTATCTGCGGGCGCTCCTGCGCCAGCCCgc                                                                                                                                    >  1:1261882/1‑142 (MQ=255)
           cgcCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGTTGCGCCAGCCCGt                                                                                                                                    <  2:1050797/143‑1 (MQ=255)
            gcCGTCTGACAGAAGAGACGGTGGCGAGGGGTCGGGCAGGGTATGAGTGATGATACACGGACCGCCGGACGCACGAAGGAACAGCCCAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGtt                                                                                                                                   <  1:424577/143‑1 (MQ=255)
            gcCGTCGGAGAGAAGATACCGTCGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGtt                                                                                                                                   <  2:244982/143‑1 (MQ=255)
             ccGGCTGATAGAAGATCCGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTg                                                                                                                                  <  1:450867/143‑1 (MQ=255)
                tCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGCACTTAAAAACGATCGTTATCTGCGCGCGCTGCAGCGCCAGCCCGTCGATg                                                                                                                               >  2:1169150/1‑143 (MQ=255)
                  tGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGGTgca                                                                                                                             >  2:662699/1‑141 (MQ=255)
                  tGATAGAAGATACGGTGGCGAAGTGTCGGGCAGTGTATGAGGGATGAGAGACTGACAGGCTGACGCACTAAGGAACAGCCAAAAGCACCGAACTTAAAAACGAGCGTTATCGGCGGCCGCCGCTGCGCCAGCGCGtgcacctc                                                                                                                             >  2:547403/1‑136 (MQ=255)
                   gATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTAATGTCa                                                                                                                            >  2:86307/1‑143 (MQ=255)
                    aTAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCAc                                                                                                                           >  1:185809/1‑143 (MQ=255)
                      agaagaTACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTTACTc                                                                                                                         >  1:185404/1‑143 (MQ=255)
                       gaagaTACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACAGCCTGACGCACGAAGGGACAGCCAAAATGGCCGAACGTAAAAACGATCGTTGTCTGCGAGCGCTGCTGCGCCCGCCCGGCGATGTCCCTcc                                                                                                                        >  2:817297/1‑143 (MQ=255)
                        aagaTACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGCTGTCACTCCa                                                                                                                       >  2:1229675/1‑143 (MQ=255)
                        aagaTACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCa                                                                                                                       >  1:619738/1‑143 (MQ=255)
                        aagaTACGGTGCGATGTGTCGGGCAGAGTATGAGGGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAg                                                                                                                      >  2:76807/1‑143 (MQ=255)
                         gggtACGGGGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAg                                                                                                                      <  2:422614/140‑1 (MQ=255)
                          gctACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGt                                                                                                                     <  1:752500/141‑1 (MQ=255)
                           atccGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTa                                                                                                                    <  2:169382/140‑1 (MQ=255)
                           aTACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGCCCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCCGCCCGTTGATGTCACTCCAGTa                                                                                                                    >  1:658804/1‑143 (MQ=255)
                            tACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGATCAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTAt                                                                                                                   <  2:1030787/143‑1 (MQ=255)
                            tACGGTGGCGAGGTGGCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAGATGACCGAGCTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTAATGTCACTCCAGTAt                                                                                                                   <  2:1080140/143‑1 (MQ=255)
                             aCGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATg                                                                                                                  <  1:114796/143‑1 (MQ=255)
                              gggTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATgg                                                                                                                 <  2:374419/142‑1 (MQ=255)
                               ggtggtGATGTCGCGGCCAGAGTAGTCGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGCCCGACCTTAAAATCGATCGTTATCTGCGTGCGCTGCCGCGCCATTCCGTTGATGTCACTCCAGTATGga                                                                                                                <  1:380307/143‑1 (MQ=255)
                                   gCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTGATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGgatgat                                                                                                            <  2:482337/143‑1 (MQ=255)
                                    cGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCAAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGTTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGatgatg                                                                                                           <  2:838893/143‑1 (MQ=255)
                                       tgtgTCGGGAAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATTACCGAACTTAAATATGATCGTTATCTGCGGGCGATGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATgcgc                                                                                                        <  1:128650/143‑1 (MQ=255)
                                       cctgtCGGGCAGCGTATGAGTGATGATACACTGACCGCCTGACGCAGTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCAGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATgcgc                                                                                                        <  2:463655/141‑1 (MQ=255)
                                        gtgtCTGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGCAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGcc                                                                                                       <  2:404206/143‑1 (MQ=255)
                                         tgtCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCa                                                                                                      <  1:771514/143‑1 (MQ=255)
                                           tCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGCATGGATGATGCGCCAgg                                                                                                    >  2:494180/1‑143 (MQ=255)
                                              ggCAGAGTATGAGTGATGATACACTGACGGCCTGACGCACTAAGGAACAGCCACAATGACCGAACTTAAAACCGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGGTGATCTCACTCCAGTCTGTCTGATGCGCCAGCCgg                                                                                                 >  2:921242/1‑143 (MQ=255)
                                              ggCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCgg                                                                                                 <  1:46877/143‑1 (MQ=255)
                                              ggCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCgg                                                                                                 >  1:348961/1‑143 (MQ=255)
                                               gCAGAGTATGAGTGATGATACACTGACCGCCTCTCGCACTGACGAACAGCCCAGATGAACGATCTGAAAAACGATCGTTATCTGCGGGCGCTGGTGCGGGAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCggg                                                                                                <  1:1140021/143‑1 (MQ=255)
                                                  ggttATGAGTCATGATACGCGTACCGCCTGACGCACTTAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCg                                                                                             <  2:486686/140‑1 (MQ=255)
                                                             atgatTCACTGACCGCCTGACGGACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCATTCCAGTATGGATGATGCGCCTGGCGGGTCGCTATCTACCgg                                                                                  <  1:1176588/143‑1 (MQ=255)
                                                               gatACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGaa                                                                                <  1:816108/143‑1 (MQ=255)
                                                                atACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAat                                                                               >  1:747364/1‑143 (MQ=255)
                                                                acacacTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAat                                                                               <  1:254951/141‑1 (MQ=255)
                                                                 tACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAata                                                                              <  2:727979/143‑1 (MQ=255)
                                                                   cacTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCCCCAGGCGGGTCGCTACCTACCGGACtata                                                                            >  2:559421/1‑143 (MQ=255)
                                                                    acTGACCGACTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGGCCCAGGCGGGTCGCCATCTACCGGAATATaa                                                                           >  2:784657/1‑143 (MQ=255)
                                                                     cagACCGCCGGACGCACTCAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATaaa                                                                          <  1:466754/141‑1 (MQ=255)
                                                                     cTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATaaa                                                                          <  2:228411/143‑1 (MQ=255)
                                                                      tGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGACGCGCCAGGCGGGTCGCTATCTACCGGGAGATAAAg                                                                         >  1:623371/1‑143 (MQ=255)
                                                                       gACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGc                                                                        >  2:45321/1‑143 (MQ=255)
                                                                          cGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGGCACTCCAGTATGGATGATGCCCCAGGCGGGTCGCTATCTACCGAAATATAAAGCCAc                                                                     >  2:42301/1‑143 (MQ=255)
                                                                          cGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGGTGATGTCACTCCAGTATGGATGGTGCGCCAGGCGGGCGGCTATCTACTGGAACATGAAGCCAc                                                                     >  2:1005463/1‑143 (MQ=255)
                                                                            ccTGACGCACTCAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCAcgc                                                                   <  1:178122/143‑1 (MQ=255)
                                                                            ccTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTGGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCAAcga                                                                   >  2:12881/1‑142 (MQ=255)
                                                                              ggACGCACAAAGGAACAGCCAGAATGGCAGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCAcgcgc                                                                 <  1:41523/142‑1 (MQ=255)
                                                                               gACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATATTAATGAAGAGAAAGCCCcgcgcg                                                                >  1:1199533/1‑143 (MQ=255)
                                                                                aCGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTGTAGTAATGACGGGCCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACgcgcgc                                                               >  2:269295/1‑143 (MQ=255)
                                                                                  gCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGCATATAAAGCCACGCGCGccc                                                             >  1:1015960/1‑143 (MQ=255)
                                                                                  gCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGATATCTGAGGGCGCTGCTGCGCCAGCCCGTTGATATCACTCCAGTATGTATAATACACCAGGCAGGCCGCTATCAACCATAAAATAAAACCAcgcaacaca                                                             >  2:970959/1‑137 (MQ=255)
                                                                                   cACTAAGGAACAGACTAGATGGCGGAATTTATAATTGATCGGTATCTGCGGGCGCGGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCca                                                            <  1:1065730/143‑1 (MQ=255)
                                                                                    aCTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCGGGCGGGGCGCTATCTACCGGAATATAAAGCCACGCGCGcccgc                                                           >  1:707039/1‑141 (MQ=255)
                                                                                    aCTAAGGAACAGCCAAAATGACCGAACTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGAGGCGCCAGGCGGGTCGCGATTTACCGAGATTTAAGGCCCCGGGCGCCccgg                                                          >  1:231308/1‑143 (MQ=255)
                                                                                      tAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCcaggc                                                         >  1:576299/1‑143 (MQ=255)
                                                                                      tAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAAAATAGATGCGCGGGCGACcaggt                                                         >  1:582773/1‑142 (MQ=255)
                                                                                          gAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGACTATAAAGCCACGCGCGCCCcggcaggc                                                     >  1:817956/1‑143 (MQ=255)
                                                                                              agtcAAAATGACGTGACTGAAATACGGTCGTTATCTGCGGGCGCTGCTGCGTCAGTCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCCGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGAtt                                                 <  2:963076/140‑1 (MQ=255)
                                                                                                ccAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGCTGCGCCCGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAAGAGAtgat                                               >  1:27579/1‑140 (MQ=255)
                                                                                                ccAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGTCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGAtttt                                               >  2:146684/1‑143 (MQ=255)
                                                                                                      tGACCGAACTTAAAGACGATCGTTCTCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCg                                         <  2:190295/143‑1 (MQ=255)
                                                                                                      tGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGcc                                                              <  1:834195/122‑1 (MQ=255)
                                                                                                           gAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATACAGCCACGCGCGCCCATGCAGGCGATCTTATGTCGCtgtg                                    >  1:1167689/1‑143 (MQ=255)
                                                                                                           caaCTTAAAAACGGTCGGTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCtgtg                                    <  1:165818/142‑1 (MQ=255)
                                                                                                              cTTAAAAACGGTCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGGCACGCCAGTATGGATGATGCGTCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCaa                                 <  2:106325/143‑1 (MQ=255)
                                                                                                               ttAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGGTGCGCCAGGCGGGTCGCTATCTACCGGAATATCAAGCCACGCGCGCCCAGGCAGGCGACAGTATGTCGCTGTGCagg                                >  1:137026/1‑141 (MQ=255)
                                                                                                               gtgaagaCGAGCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTGCAGTATGGATGATGCGCCGGGCGGGTCGCTAGCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCaaa                                <  2:3651/137‑1 (MQ=255)
                                                                                                                tAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGCTTTTATGTCGCTGTGCaaaa                               >  2:1230839/1‑143 (MQ=255)
                                                                                                                  caatcGATCGTTATCTGCGGGTGCTTCTGCGCCAGCCGGTTGATGTCAATCCAGTATGGATGATGCGCCAGGCGCGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAAc                             <  2:561805/139‑1 (MQ=255)
                                                                                                                  aaaaCGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATGTTATGTCGCTGTGCGAAAAc                             >  2:845913/1‑143 (MQ=255)
                                                                                                                   aaaCGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGGCGCTATCTACCGGAATATAAGGCCACGCGCGCCCAGGCAGGCGCTGTTATGTCGCTGTGCGAAAACg                            >  1:730765/1‑143 (MQ=255)
                                                                                                                    aaCGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGc                           >  2:1005650/1‑143 (MQ=255)
                                                                                                                     aCGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCGCTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGcc                          <  2:986069/143‑1 (MQ=255)
                                                                                                                     aCGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGcc                          <  1:86307/143‑1 (MQ=255)
                                                                                                                            ttATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCCAAAACGCCGAGCTgg                   >  1:894485/1‑143 (MQ=255)
                                                                                                                              aTCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATCTTATGTCGCTGTGCAAAAACGCCGAGCTGGcc                 >  2:724468/1‑142 (MQ=255)
                                                                                                                               tCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAACAACGCCGAGCTGGCGt                >  2:1056149/1‑143 (MQ=255)
                                                                                                                                cTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTCTGTCGCTGTGCCAAAACGCGGAGCTGGCGTg               >  1:935967/1‑143 (MQ=255)
                                                                                                                                  gctggCGCAGCAGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCg             >  1:1203295/4‑143 (MQ=255)
                                                                                                                                  gCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGTGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCCAAACCGCCGAGCTGGcgcgcg             >  2:571347/1‑143 (MQ=255)
                                                                                                                                   ggggCGCGGCCGCGCCAGCCCGTCGATGTCACCCCAGTATGGATGATGCGCCCGGCGGGTCGCTATCTACCGGAAGATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGa            <  2:510833/142‑1 (MQ=255)
                                                                                                                                   cGGGCGCTGCTGCGCCGGCCCGTTGATGTCACTCCAGTATGGTTGATGCGCCAGGCGGGCCGCTCACTAACGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTCATGTCGCTGTGCAAAAACGCCGAGCTGCCGGGcca            >  2:405562/1‑141 (MQ=255)
                                                                                                                                   cGGGCGCGCCGGCGCGAGCCCGTGGATGTCCCTCCAGGACGGATGATGCGCCAGGCGGGTCGCTAGCTACCGGAATATCAAGCCACGCGCGCCCAGGTAGGCGATTTTATGGCGCTGTGCAAAAACGCCGAGCTGGCGTGCGa            <  2:775164/143‑1 (MQ=255)
                                                                                                                                    gggCGCGGCTGCGCCGGCCCGTTGATGTCACTCCAGTATGGATGGTGCGCCCGCCACCTCGCTATCTACCGGAAGATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGaa           <  1:588864/143‑1 (MQ=255)
                                                                                                                                      gcgcTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGt         <  1:260244/143‑1 (MQ=255)
                                                                                                                                       cgcTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTg        >  1:576601/1‑143 (MQ=255)
                                                                                                                                            ctgcTCCAGCCCGTTGATGTCACCCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACttt   <  1:985295/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGTTATCTGCGGGCGCTGCTGCGCCAGCCCGTTGATGTCACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTG  >  NC_000913/4197615‑4197898

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