Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,694,196 T→A 5.7% E185V (GAA→GTA)  U712_13885 ← Porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,694,1960TA 5.7% 217.9 / 3.8 70E185V (GAA→GTA) U712_13885Porphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base T (36/30);  new base A (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

GTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCATTTAATTGTAAGGTCAGGG  >  CP006881/2694070‑2694333
                                                                                                                              |                                                                                                                                         
gTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGa                                                                                                                                  <  6:124037/136‑1 (MQ=255)
gTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGa                                                                                                                                  >  5:124037/1‑136 (MQ=255)
gTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACAt                                                                                                                               >  8:186645/1‑139 (MQ=255)
gTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACAt                                                                                                                               >  4:99951/1‑139 (MQ=255)
gTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACAt                                                                                                                               >  3:314470/1‑139 (MQ=255)
gTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACAt                                                                                                                               >  2:342877/1‑139 (MQ=255)
gTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATACGGTTACGACAt                                                                                                                               >  5:177366/1‑139 (MQ=255)
gTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATACGGTTACGACAt                                                                                                                               >  5:171639/1‑139 (MQ=255)
gTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAAAAGCTCTTCATCCGATTATCGGCACTCAATCGCCAGGGCTCCCTGCCCCAAAGCAGTCAAACAGCGCACAGGCTCAAGCAATTCGGTTACGACa                                                                                                                                >  8:242694/1‑138 (MQ=255)
 tCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATc                                                                                                                              >  1:156882/1‑139 (MQ=255)
 tCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATc                                                                                                                              >  2:97710/1‑139 (MQ=255)
           atatTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCa                                                                                                                     <  7:186645/138‑1 (MQ=255)
            tatTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAAc                                                                                                                   <  1:342877/139‑1 (MQ=255)
            tatTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAAc                                                                                                                   <  6:130448/139‑1 (MQ=255)
            tatTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAAc                                                                                                                   <  2:10663/139‑1 (MQ=255)
            tatTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAAc                                                                                                                   <  4:314470/139‑1 (MQ=255)
            tatTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATACGGTTACGACATCTTGCTTCCAAc                                                                                                                   <  6:171639/139‑1 (MQ=255)
            tatTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATACGGTTACGACATCTTGCTTCCAAc                                                                                                                   <  6:177366/139‑1 (MQ=255)
                cATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCa                                                                                                                <  4:335573/138‑1 (MQ=255)
                cATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCa                                                                                                                <  8:132201/138‑1 (MQ=255)
                 aTCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACg                                                                                                                                   >  8:92161/1‑118 (MQ=255)
                 aTCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACg                                                                                                                                   <  7:92161/118‑1 (MQ=255)
                    tGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCACCCCATTCTg                                                                                                           >  3:184700/1‑139 (MQ=255)
                         aCTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAg                                                                                                      >  5:280818/1‑139 (MQ=255)
                             agaAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCCTTCTGGAAAGGCCg                                                                                                  >  6:181636/1‑139 (MQ=255)
                             agaAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCg                                                                                                  >  8:88172/1‑139 (MQ=255)
                                     acGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGc                                                                                          >  5:193376/1‑139 (MQ=255)
                                           aCAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTaaaa                                                                                     >  6:201669/1‑138 (MQ=255)
                                                    cATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCg                                                                             <  3:99951/137‑1 (MQ=255)
                                                            tctcGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACTTCTAGCTTCCAACCCATTCTGGAAAGGCCGGCCGCCGCTAAAATAATTGCGTCATAATCct                                                                   >  6:17884/1‑139 (MQ=255)
                                                             ctcGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCct                                                                   <  1:366890/138‑1 (MQ=255)
                                                              tcGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCctct                                                                 >  8:260925/1‑139 (MQ=255)
                                                                       aaTCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGc                                                        >  4:25311/1‑139 (MQ=255)
                                                                            ccAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCtttt                                                    >  4:216431/1‑138 (MQ=255)
                                                                            ccAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCtttt                                                    >  2:145512/1‑138 (MQ=255)
                                                                             cAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGa                                                  <  1:145512/139‑1 (MQ=255)
                                                                                 gCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTc                                              <  2:184557/139‑1 (MQ=255)
                                                                                     ccTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATa                                                                        >  5:147825/1‑109 (MQ=255)
                                                                                     ccTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATa                                                                        <  6:147825/109‑1 (MQ=255)
                                                                                     ccTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCtt                                            <  7:88172/137‑1 (MQ=255)
                                                                                     ccTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCtt                                            <  2:156882/137‑1 (MQ=255)
                                                                                       tGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCCGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTAt                                        >  6:216128/1‑139 (MQ=255)
                                                                                         ccccACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTAt                                        <  6:280818/137‑1 (MQ=255)
                                                                                           ccACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCtat                                    >  2:244133/1‑139 (MQ=255)
                                                                                            cacaGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCtata                                   >  2:78951/1‑139 (MQ=255)
                                                                                             acaGCAGGCAAACAGCGCTCATGCTCAATGAATTCGGTTACTACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCATCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  <  5:181636/139‑1 (MQ=255)
                                                                                             acaGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGATAGGGCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  <  7:242694/139‑1 (MQ=255)
                                                                                             acaGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  <  1:312914/139‑1 (MQ=255)
                                                                                               agcagGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  >  3:216250/1‑137 (MQ=255)
                                                                                               agcagGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  <  4:216250/137‑1 (MQ=255)
                                                                                               agcagGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATAtt                                 <  4:341532/138‑1 (MQ=255)
                                                                                               agcagGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATAtt                                 >  3:341532/1‑138 (MQ=255)
                                                                                                   ggCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACctc                            >  6:126125/1‑139 (MQ=255)
                                                                                                    gCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACctc                            >  4:330252/1‑138 (MQ=255)
                                                                                                     cAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTa                          <  6:193376/139‑1 (MQ=255)
                                                                                                       aaCAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAAt                        <  1:97710/139‑1 (MQ=255)
                                                                                                       aaCAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAAt                        <  3:216431/139‑1 (MQ=255)
                                                                                                           gcgcTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGc                                                        >  3:193979/1‑103 (MQ=255)
                                                                                                           gcgcTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGc                                                        <  4:193979/103‑1 (MQ=255)
                                                                                                            cgcTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACACATTCTGGAAAGGCCGGCAGCCGCTAAAATATTTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCAtt                   >  3:309075/1‑139 (MQ=255)
                                                                                                              cTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCCGCTTTTGAAGTCTTGTATCTATATTACCGCTAATCCAttt                  >  7:44005/1‑138 (MQ=255)
                                                                                                              cTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  <  7:121428/122‑1 (MQ=255)
                                                                                                              cTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  >  8:121428/1‑122 (MQ=255)
                                                                                                              cTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCAttt                  >  7:14933/1‑138 (MQ=255)
                                                                                                              cTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCAttt                  >  5:46127/1‑138 (MQ=255)
                                                                                                                  ggCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCATTTAATTg             <  6:46127/139‑1 (MQ=255)
                                                                                                                   gCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCATTTAATTGt            <  3:25311/139‑1 (MQ=255)
                                                                                                                   gCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCATTTAATTGt            <  2:250703/139‑1 (MQ=255)
                                                                                                                        aGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCGTCTGTTTCCAGCTTTTGAAGTCTTGTTTCTATATTACCTCTATTCCATTTAATTGTAAGGt       >  8:267672/1‑139 (MQ=255)
                                                                                                                        aGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCATTTAATTGTAAgg        >  5:39373/1‑138 (MQ=255)
                                                                                                                             ttCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  <  6:364023/107‑1 (MQ=255)
                                                                                                                             ttCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  >  3:9001/1‑107 (MQ=255)
                                                                                                                             ttCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  <  4:9001/107‑1 (MQ=255)
                                                                                                                             ttCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTatat                                  >  5:364023/1‑107 (MQ=255)
                                                                                                                             ttCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCATTTAATTGTa           <  6:46505/130‑1 (MQ=255)
                                                                                                                             ttCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCATTTAATTGTa           >  5:46505/1‑130 (MQ=255)
                                                                                                                             ttCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATATATATTACCTCTAATCCATTTAATTGTAAGGTCAggg  >  6:20845/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                                         
GTCCGTTTTGTATATTCATCTGTAAACTGAGAAAACAACGCCAACAGCTCTTCATCCGATTCTCGGCACTCAATCGCCAGGGCTCCCTGCCCCACAGCAGGCAAACAGCGCTCAGGCTCAAGGAATTCGGTTACGACATCTTGCTTCCAACCCATTCTGGAAAGGCCGGCAGCCGCTAAAATAATTGCGTCATAATCCTCTGTTTCCAGCTTTTGAAGTCTTGTATCTATATTACCTCTAATCCATTTAATTGTAAGGTCAGGG  >  CP006881/2694070‑2694333

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: