Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 389,169 T→C T187A (ACG→GCG)  hemB ← 5‑aminolevulinate dehydratase (porphobilinogen synthase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913389,1690TC99.2% 482.7 / ‑4.5 126T187A (ACG→GCG) hemB5‑aminolevulinate dehydratase (porphobilinogen synthase)
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (64/61);  total (64/62)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 127 total aligned reads displayed.

CCAGCAGTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGTATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGTTGTC  >  NC_000913/389022‑389313
                                                                                                                                                   |                                                                                                                                                
ccAGCAGTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCAt                                                                                                                                                <  2:1676809/150‑1 (MQ=255)
ccAGCAGTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATc                                                                                                                                               >  1:361928/1‑151 (MQ=255)
 cagcagTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCt                                                                                                                                              >  1:592530/1‑151 (MQ=255)
    cagTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTa                                                                                                                                           <  2:202690/151‑1 (MQ=255)
     agTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTaa                                                                                                                                          <  2:757559/151‑1 (MQ=255)
     agTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTaa                                                                                                                                          <  2:6858/151‑1 (MQ=255)
     agTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTaa                                                                                                                                          <  2:461596/151‑1 (MQ=255)
      gTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTaaa                                                                                                                                         <  2:978428/151‑1 (MQ=255)
        gATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAAt                                                                                                                                        <  1:1029947/150‑1 (MQ=255)
            cACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGc                                                                                                                                   >  1:379716/1‑151 (MQ=255)
              cGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAg                                                                                                                                 <  1:1603610/151‑1 (MQ=255)
               gAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGTCATAATCGCCGCATCTTTAAATCCCGCAGc                                                                                                                                <  2:297030/151‑1 (MQ=255)
               gAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGc                                                                                                                                >  2:343473/1‑151 (MQ=255)
                 aTCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGt                                                                                                                              <  1:951114/151‑1 (MQ=255)
                 aTCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGt                                                                                                                              <  1:41661/151‑1 (MQ=255)
                   cGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTcc                                                                                                                            >  1:1136298/1‑151 (MQ=255)
                    gCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTcc                                                                                                                            >  2:1360147/1‑150 (MQ=255)
                    gCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCa                                                                                                                           <  2:1088075/151‑1 (MQ=255)
                    gCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCa                                                                                                                           <  1:1000392/151‑1 (MQ=255)
                     ccTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTcc                                                                                                                            >  2:1434179/1‑149 (MQ=255)
                     ccTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAg                                                                                                                          <  2:1071613/151‑1 (MQ=255)
                      cTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAgc                                                                                                                         >  1:247765/1‑151 (MQ=255)
                      cTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAgc                                                                                                                         >  2:500185/1‑151 (MQ=255)
                        cACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAgc                                                                                                                         <  1:1217078/149‑1 (MQ=255)
                         acgacgGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGcc                                                                                                                      >  2:1647046/1‑151 (MQ=255)
                          cgacgGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCt                                                                                                                     <  1:802829/151‑1 (MQ=255)
                            acgGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTg                                                                                                                    >  2:121923/1‑150 (MQ=255)
                            acgGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTg                                                                                                                    >  1:262222/1‑150 (MQ=255)
                            acgGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTg                                                                                                                    <  2:151805/150‑1 (MQ=255)
                               gTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGAc                                                                                                                  <  2:432067/149‑1 (MQ=255)
                                      gggTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCt                                                                                                         <  1:1024547/151‑1 (MQ=255)
                                        gTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGt                                                                                                       <  1:1185459/151‑1 (MQ=255)
                                        gTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGt                                                                                                       >  2:1624212/1‑151 (MQ=255)
                                         ttCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGt                                                                                                       >  1:1693244/1‑150 (MQ=255)
                                          tCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTa                                                                                                      <  2:424075/150‑1 (MQ=255)
                                            aTCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCt                                                                                                   >  1:799387/1‑151 (MQ=255)
                                                gATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCCCGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGTATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTgg                                                                                                 <  2:1041514/149‑1 (MQ=255)
                                                gATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTgg                                                                                                 <  2:850283/149‑1 (MQ=255)
                                                 aTAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCg                                                                                              >  1:1468997/1‑151 (MQ=255)
                                                  tAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGt                                                                                             >  2:1121594/1‑151 (MQ=255)
                                                  tAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGt                                                                                             >  2:1209979/1‑151 (MQ=255)
                                                   aGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGt                                                                                             <  1:1646939/150‑1 (MQ=255)
                                                   aGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTc                                                                                            <  2:379743/151‑1 (MQ=255)
                                                   aGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTc                                                                                            >  2:477934/1‑151 (MQ=255)
                                                          gCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATcgcg                                                                                     <  1:282419/151‑1 (MQ=255)
                                                           cGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCgcgg                                                                                    >  2:854702/1‑151 (MQ=255)
                                                            ggTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCgcgg                                                                                    <  1:307248/150‑1 (MQ=255)
                                                                gCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGa                                                                                >  2:190561/1‑150 (MQ=255)
                                                                    tttAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAgggg                                                                           >  1:58313/1‑151 (MQ=255)
                                                                      tAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGc                                                                          <  2:1646590/150‑1 (MQ=255)
                                                                      tAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCg                                                                         >  1:1713442/1‑151 (MQ=255)
                                                                        aTGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGAt                                                                       <  1:234383/151‑1 (MQ=255)
                                                                         tGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATg                                                                      >  2:1387056/1‑151 (MQ=255)
                                                                          gcgcTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATg                                                                      >  1:619814/1‑150 (MQ=255)
                                                                                 cGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGc                                                              >  1:1598496/1‑151 (MQ=255)
                                                                                   gCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCGTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGc                                                              <  1:963472/149‑1 (MQ=255)
                                                                                     aGCTTCACGGAACGGGCCATAAAAGGAGGGGGCGCACTTGGTCGAATTCGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCt                                                          <  1:319654/151‑1 (MQ=255)
                                                                                     aGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCt                                                          >  2:1486023/1‑151 (MQ=255)
                                                                                       cTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTgc                                                        >  2:1341517/1‑151 (MQ=255)
                                                                                        ttCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTgc                                                        >  2:290447/1‑150 (MQ=255)
                                                                                         tCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTgca                                                       <  1:897119/150‑1 (MQ=255)
                                                                                         tCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTgc                                                        <  2:1666374/149‑1 (MQ=255)
                                                                                          cACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTgcagc                                                     <  2:849734/151‑1 (MQ=255)
                                                                                             ggAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCagcagc                                                  >  2:1343904/1‑151 (MQ=255)
                                                                                             ggAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCagcagc                                                  >  2:897838/1‑151 (MQ=255)
                                                                                              gAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCt                                                 <  2:361956/151‑1 (MQ=255)
                                                                                                 cGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGc                                               >  1:487885/1‑150 (MQ=255)
                                                                                                  gggCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCaa                                             <  1:1049023/151‑1 (MQ=255)
                                                                                                     ccATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAAcc                                           >  2:814860/1‑150 (MQ=255)
                                                                                                       aTAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACg                                        <  1:1637257/151‑1 (MQ=255)
                                                                                                        tAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACgg                                       >  1:1670805/1‑151 (MQ=255)
                                                                                                         aaaaGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACgg                                       <  1:622209/150‑1 (MQ=255)
                                                                                                         aaaaGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGgc                                      <  2:1494692/151‑1 (MQ=255)
                                                                                                           aaGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGgctt                                    <  2:1132800/151‑1 (MQ=255)
                                                                                                            aggaggaggCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGgcttg                                   >  1:1584735/1‑151 (MQ=255)
                                                                                                                ggaggCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGcttgcttg                               >  1:240477/1‑151 (MQ=255)
                                                                                                                 gaggCGAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCttgcttgc                              >  1:1572861/1‑151 (MQ=255)
                                                                                                                      gAACTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTaa                          >  2:1236144/1‑150 (MQ=255)
                                                                                                                        aCTTGGTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAAtt                       >  1:1510953/1‑151 (MQ=255)
                                                                                                                            ggTCGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTc                    <  2:771460/150‑1 (MQ=255)
                                                                                                                               cGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTcc                   >  2:761722/1‑148 (MQ=255)
                                                                                                                               cGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTcc                   >  1:385181/1‑148 (MQ=255)
                                                                                                                               cGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTcc                   >  1:179425/1‑148 (MQ=255)
                                                                                                                               cGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCag                 <  2:616672/150‑1 (MQ=255)
                                                                                                                               cGAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCag                 <  2:1584831/150‑1 (MQ=255)
                                                                                                                                gAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATTGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCagag               <  2:419302/151‑1 (MQ=255)
                                                                                                                                gAATACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCagag               <  1:986086/151‑1 (MQ=255)
                                                                                                                                   tACGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTcg            >  1:840129/1‑151 (MQ=255)
                                                                                                                                    aCGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTcgc           <  2:936116/151‑1 (MQ=255)
                                                                                                                                     cGACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTcgc           >  1:552196/1‑150 (MQ=255)
                                                                                                                                      gACATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGcgt         >  2:495802/1‑151 (MQ=255)
                                                                                                                                       aCATAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGcgtc        >  1:1456156/1‑151 (MQ=255)
                                                                                                                                         aTAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGcgtcgt      >  2:582243/1‑151 (MQ=255)
                                                                                                                                         aTAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGcgtcgt      >  1:489302/1‑151 (MQ=255)
                                                                                                                                          tAATCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGtt     >  1:1547013/1‑151 (MQ=255)
                                                                                                                                           aaTCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGtt     <  1:1360055/150‑1 (MQ=255)
                                                                                                                                           aaTCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGTTg    >  2:146424/1‑151 (MQ=255)
                                                                                                                                            aTCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGTTg    >  2:585406/1‑150 (MQ=255)
                                                                                                                                             tCGCCGCATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGTTGTc  >  1:1220886/1‑151 (MQ=255)
                                                                                                                                              cgccgcATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGTTGTc  <  1:763751/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCAGCAGTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATGCGCTTCCGGCAGCTTCACGGAACGGGCCATAAAAGGAGGAGGCGAACTTGGTCGAATACGACATAATCGCCGTATCTTTAAATCCCGCAGCGTCCAGCGCCTGACGAATCGCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGTTGTC  >  NC_000913/389022‑389313

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