Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,791,864 C→A 7.7% intergenic (‑147/+48) hemE ← / ← gltD uroporphyrinogen decarboxylase/glutamate synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,791,8640CA 7.7% 199.2 / 2.4 64intergenic (‑147/+48)hemE/gltDuroporphyrinogen decarboxylase/glutamate synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base A (3/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCGAACCGCGCACCATGTCGCCACCGGCA  >  NC_002947/5791722‑5792011
                                                                                                                                                |                                                                                                                                                   
gCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCt                                                                                                                                               >  2:256232/1‑149 (MQ=255)
  aTCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCtt                                                                                                                                              >  4:41002/1‑148 (MQ=255)
  aTCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCtt                                                                                                                                              >  2:344960/1‑148 (MQ=255)
   tCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTaa                                                                                                                                            <  7:206438/149‑1 (MQ=255)
         aaaaTAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGc                                                                                                                                      <  7:797/149‑1 (MQ=255)
          aaaTAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTc                                                                                                                                                     >  1:333333/1‑133 (MQ=255)
          aaaTAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTc                                                                                                                                                     <  2:333333/133‑1 (MQ=255)
                   gggCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCggg                                                                                                                             >  5:364748/1‑148 (MQ=255)
                   gggCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCggg                                                                                                                             >  2:77535/1‑148 (MQ=255)
                   gggCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCggg                                                                                                                             >  2:58933/1‑148 (MQ=255)
                    ggCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCTATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGc                                                                                                                            <  1:77535/149‑1 (MQ=255)
                     gCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCAt                                                                                                                          <  4:310516/149‑1 (MQ=255)
                       aTTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCCGCTCCCGCAGGGTTAAGCGCAGAACCCGTGGGGGCAGGCATGc                                                                                                                        <  3:206562/149‑1 (MQ=21)
                       aTTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCCGCTCCCGCAGGGTTAAGCGCAGAACCCGTGGGGGCAGGCATGc                                                                                                                        <  3:211488/149‑1 (MQ=21)
                                gACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGaa                                                                                                                 >  5:194517/1‑147 (MQ=255)
                                gACACCAGCGCAAAAGGCACGGCCATGGCCATGTGGTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACa                                                                                                               <  3:302174/149‑1 (MQ=255)
                                       gcgcAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGaa                                                                                                        <  3:41002/149‑1 (MQ=255)
                                        cgcAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATTCCCGCGAACAGTTCGAAg                                                                                                       >  5:215592/1‑149 (MQ=255)
                                        cgcAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAg                                                                                                       <  6:282821/149‑1 (MQ=255)
                                           aaaaGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGaa                                                                                                        >  5:161722/1‑145 (MQ=255)
                                           aaaaGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGaa                                                                                                        <  6:161722/145‑1 (MQ=255)
                                           aaaaGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCt                                                                                                     >  3:322037/1‑148 (MQ=255)
                                            aaaGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCa                                                                                                   <  1:344960/149‑1 (MQ=255)
                                            aaaGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCa                                                                                                   <  6:364748/149‑1 (MQ=255)
                                                  aCGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACACGTCGAAGCTCAGACTTc                                                                                             >  1:174444/1‑149 (MQ=255)
                                                    ggccatggccatGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGa                                                                                           <  1:256232/149‑1 (MQ=255)
                                                    ggccatggccatGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGa                                                                                           <  3:198977/149‑1 (MQ=255)
                                                       catggccatgTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGt                                                                                        <  1:58933/149‑1 (MQ=255)
                                                            ccatgTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTc                                                                                             >  7:180183/1‑139 (MQ=255)
                                                            ccatgTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTc                                                                                             <  8:180183/139‑1 (MQ=255)
                                                             catgTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCa                                                                                  <  2:308684/149‑1 (MQ=255)
                                                              atgTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGTAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAg                                                                                 <  3:198546/149‑1 (MQ=255)
                                                                   ttttGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGAATGCCTGCCCCCACGGGTTATGAGCTTAACCCTGCGGGACCGGGCCTGCCCGCCAACCGGTCCAAGCTCAGACTTCGAGGTATTCAAGGAAAc                                                                            >  5:41600/1‑149 (MQ=255)
                                                                      tGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATAccc                                                                          <  4:65950/148‑1 (MQ=255)
                                                                         cTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGccc                                                                                                                      >  3:362234/1‑101 (MQ=255)
                                                                         cTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGccc                                                                                                                      <  4:362234/101‑1 (MQ=255)
                                                                          tATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCgg                                                                      >  3:44487/1‑148 (MQ=255)
                                                                            ttGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGa                                                                                           <  7:288682/125‑1 (MQ=255)
                                                                            ttGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGa                                                                                           >  8:288682/1‑125 (MQ=255)
                                                                            ttGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCt                                                                         <  6:12810/143‑1 (MQ=255)
                                                                            ttGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCt                                                                         >  5:12810/1‑143 (MQ=255)
                                                                                   ccACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACg                                                            <  5:243217/149‑1 (MQ=255)
                                                                                    cacaGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGAATGCCTGCCCCCACGGGTTATGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACACGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCATCCGCAGCCGTACGc                                                           >  8:292479/1‑149 (MQ=255)
                                                                                          tCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCTCGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAg                                                                                 >  1:260712/1‑121 (MQ=255)
                                                                                          tCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAg                                                                                 <  2:260712/121‑1 (MQ=255)
                                                                                            aaGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGaa                                                   <  6:325060/149‑1 (MQ=255)
                                                                                                gTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATc                                               <  3:260334/149‑1 (MQ=255)
                                                                                                gTTTTTGTATAAAAAAGCTTCGCGGGCATGCCTGCCCCCACGTGTTCTGCGCTTAACCCTGCGTGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATc                                               <  1:52268/149‑1 (MQ=255)
                                                                                                 tttttGTATACAAAAGCTTCGGGGGCATGCCTGCCCCCACGGGTGCTGCGCTTAACCCTGCGGGAGCGGGAATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCg                                              <  6:41600/149‑1 (MQ=255)
                                                                                                         taCAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTcac                                       >  7:215633/1‑148 (MQ=255)
                                                                                                                tgttCGCGGGCATGCCCGCTCCCGCAGGGTTAAGCGCAGAACCCGTGGGGGCAGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGt                               >  7:305101/3‑149 (MQ=18)
                                                                                                                 cTTCGCGTGCATGCATGCCCCCACTGGTTCTGCCCTTAACCCTGCGGGAGCGGGCAGGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTc                              <  6:215592/149‑1 (MQ=255)
                                                                                                                 cTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAgg                                                                                         >  2:123771/1‑90 (MQ=255)
                                                                                                                 cTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTc                              >  8:213678/1‑149 (MQ=255)
                                                                                                                  ttCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTcc                             >  8:169601/1‑149 (MQ=255)
                                                                                                                  ttCGCGGGCATGCCCGCTCCCGCAGGGTTAAGCGCAGAACCCGTGGGGGCAGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTc                              <  8:305101/148‑1 (MQ=18)
                                                                                                                         gCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCGAACcgc                      >  7:148312/1‑149 (MQ=255)
                                                                                                                         gCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCGAACcgc                      <  4:44487/149‑1 (MQ=255)
                                                                                                                          cATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCCCCACCACGTCCGAACcgcg                     >  6:297795/1‑149 (MQ=255)
                                                                                                                            aTGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAgg                                                                                         <  2:355443/80‑1 (MQ=255)
                                                                                                                             tgcctgccCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCGAACCGCGCa                   >  8:22994/1‑149 (MQ=255)
                                                                                                                                 tgccCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCGAACCGCGCACCAt               >  5:88362/1‑149 (MQ=255)
                                                                                                                                     cccACGGGTTCTGCGCTTAACCCTGCGGCAGCTCGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACAACCAGGTCCGAACCGCGCACCATGTCg           >  7:154412/1‑149 (MQ=255)
                                                                                                                                      ccACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCGAACCGCGCACCATGTCGc          >  7:178562/1‑149 (MQ=255)
                                                                                                                                            ggTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCGAACCGCGCACCATGTCGCCACCg     <  6:88362/149‑1 (MQ=255)
                                                                                                                                               tCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTAGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCTCAGTCACCACCAGGTCCGAACCGCGCACCATGTCGCCACCGGc   >  4:200278/1‑148 (MQ=255)
                                                                                                                                                cTGCTCTTAACCATGCGGGAGCGGGCATGCACGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGCTACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGGCACCACCAGGGCCGAACCGCGCACCCTGTCGCCACCGGCa  <  5:297795/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCGAACCGCGCACCATGTCGCCACCGGCA  >  NC_002947/5791722‑5792011

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: