Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,958,082 T→A 8.5% L281Q (CTG→CAG)  PP_2589 → aldehyde dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,958,0820TA 8.5% 168.8 / 5.6 59L281Q (CTG→CAG) PP_2589aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base A (2/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

ACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGTGGCTGCCCGG  >  NC_002947/2957935‑2958229
                                                                                                                                                   |                                                                                                                                                   
aCCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAggacctg                                                                                                                                                    >  5:252681/1‑149 (MQ=255)
           gACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGTCCAGGTCCTGGcag                                                                                                                                          <  2:54044/148‑3 (MQ=255)
           gACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGc                                                                                                                                         >  3:31577/1‑149 (MQ=255)
           gACTGCCTGGACATCAGCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGc                                                                                                                                         <  6:252681/149‑1 (MQ=255)
                 cTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGTCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGaaaa                                                                                                                                   <  2:171928/149‑1 (MQ=255)
                 cTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGaaaa                                                                                                                                   <  7:263204/149‑1 (MQ=255)
                 cTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATCCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGaaaa                                                                                                                                   <  5:314977/149‑1 (MQ=255)
                        tCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGgccgcc                                                                                                                            >  7:32014/1‑149 (MQ=255)
                          aCCGGCTCCACCCAGGTGGGCAAGTAGTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCtt                                                                                                                          >  6:155886/1‑149 (MQ=255)
                                   aCCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGAGGAAAAGGCCGCCTTCGGCAttt                                                                                                                  >  8:338934/1‑148 (MQ=255)
                                      cAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTttctt                                                                                                              >  1:216707/1‑149 (MQ=255)
                                        ggTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCa                                                                                                            <  1:185938/149‑1 (MQ=255)
                                         gTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCaa                                                                                                           <  5:155886/149‑1 (MQ=255)
                                          tGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAAc                                                                                                          <  2:132658/149‑1 (MQ=255)
                                               aaGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAggg                                                                                                     >  5:53417/1‑149 (MQ=255)
                                                   aCTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAg                                                                                                 >  8:25504/1‑149 (MQ=255)
                                                   aCTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTATCTTCATCCAGTGCGAg                                                                                                 >  7:143177/1‑149 (MQ=255)
                                                    cTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAgg                                                                                                >  5:167081/1‑149 (MQ=255)
                                                     ttCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGgt                                                                                               >  1:143449/1‑149 (MQ=255)
                                                         tGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGacctggacctgg                                                                                                                                             <  3:81665/99‑1 (MQ=255)
                                                         tGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGacctggacctgg                                                                                                                                             >  4:81665/1‑99 (MQ=255)
                                                            aGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCg                                                                                        <  2:216707/149‑1 (MQ=255)
                                                            aGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCg                                                                                        <  6:167081/149‑1 (MQ=255)
                                                             gTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCgg                                                                                       <  8:32014/149‑1 (MQ=255)
                                                              tACTCGGCACAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGc                                                                                      <  1:100881/149‑1 (MQ=255)
                                                                 tCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCa                                                                                    >  6:368655/1‑148 (MQ=255)
                                                                 tCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCa                                                                                    >  3:156577/1‑148 (MQ=255)
                                                                           ccAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTg                                                                         >  2:234685/1‑149 (MQ=255)
                                                                                 tCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTg                                                                   >  6:203211/1‑149 (MQ=255)
                                                                                 tCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTg                                                                   >  6:208238/1‑149 (MQ=255)
                                                                                 tCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTg                                                                   >  7:29412/1‑149 (MQ=255)
                                                                                   aaGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTg                                                                   <  8:29412/147‑1 (MQ=255)
                                                                                         gtgtGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCg                                                                                        <  1:13256/120‑1 (MQ=255)
                                                                                         gtgtGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCg                                                                                        >  2:13256/1‑120 (MQ=255)
                                                                                              gCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATc                                                       >  4:102011/1‑148 (MQ=255)
                                                                                              gCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGATGCGGAAAAGGCCGCCTTCGGCATATTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCTGCTGTACGTGCAGCGTGCCATcc                                                      >  4:23280/1‑149 (MQ=255)
                                                                                                tGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTc                                                                                         >  7:133148/1‑112 (MQ=255)
                                                                                                tGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTc                                                                                         <  8:133148/112‑1 (MQ=255)
                                                                                                  gAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGTAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGa                                                  <  5:368655/149‑1 (MQ=255)
                                                                                                       cggcggCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAg                                              <  6:53417/148‑1 (MQ=255)
                                                                                                       cggcggCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAg                                              <  7:338934/148‑1 (MQ=255)
                                                                                                         gcggcAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTg                                                            <  8:240794/132‑1 (MQ=255)
                                                                                                         gcggcAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTg                                                            >  7:240794/1‑132 (MQ=255)
                                                                                                            gcAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATc                                        <  3:1314/149‑1 (MQ=255)
                                                                                                              aaGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAAACAGGCCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGa                                      >  3:352544/1‑149 (MQ=255)
                                                                                                                gagCCCAAACCTGGTGTTGGACAAGTGCCAGGTCCTGGACCTGGCTGCGGAATAGGCCGCCTTCGGCATTTTCTTCAACCAgg                                                                                                      <  6:75971/83‑1 (MQ=255)
                                                                                                                gagCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAgc                                    <  7:25504/149‑1 (MQ=255)
                                                                                                                  gCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAgcgc                                  >  7:80886/1‑149 (MQ=255)
                                                                                                                    ccAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCAt                                         >  4:2067/1‑140 (MQ=255)
                                                                                                                    ccAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCAt                                         <  3:2067/140‑1 (MQ=255)
                                                                                                                    ccAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACAAGGGCGAGGTGAGTTCGGCCAACTCGCGGCTGTACGTGCAGCGAGCCATCCACGATGAGTTCATCGAGCGCCt                                >  5:34226/1‑149 (MQ=255)
                                                                                                                     cAAAGCTGGTTTTAGACAACTGCCAGGACCTTGACGTGGCTGCTGATAAGGACTCCTTATGCATTTTCTTCAACCAGGGCTATGTGTGTTCGTACAACTCGCGGCTGTACGTGCAGCGTGCAATCCACGATGAGTTCATCGAGCGCCTg                               <  8:143177/149‑1 (MQ=255)
                                                                                                                       aaCCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGTAAAAGTCCTCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCa                             <  1:294823/149‑1 (MQ=255)
                                                                                                                             gtgtTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCaa                       <  4:156577/149‑1 (MQ=255)
                                                                                                                                   gACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCg                 >  5:273405/1‑149 (MQ=255)
                                                                                                                                    aCAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGc                >  7:279760/1‑149 (MQ=255)
                                                                                                                                    aCAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCCCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGc                <  4:31577/149‑1 (MQ=255)
                                                                                                                                       gcagCCAGGTCCAGGTCCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAg             >  6:122163/4‑149 (MQ=255)
                                                                                                                                        cagCCAGGTCCAGGTCCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAg                            <  2:81912/131‑1 (MQ=255)
                                                                                                                                        cagCCAGGTCCAGGTCCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAg                            >  1:81912/3‑133 (MQ=255)
                                                                                                                                        cagCCAGGTCCAGGTCCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGt            <  5:122163/147‑1 (MQ=255)
                                                                                                                                             aGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGTGGCTg       <  3:102011/149‑1 (MQ=255)
                                                                                                                                                  ctggacctggCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGTGGCTGCCCgg  >  7:249720/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGTGGCTGCCCGG  >  NC_002947/2957935‑2958229

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