Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,958,079 A→T 9.7% D280V (GAC→GTC)  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,0790AT 9.7% 159.5 / 5.4 59D280V (GAC→GTC) PP_2589aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base T (2/4);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGTGGCTG  >  NC_002947/2957931‑2958224
                                                                                                                                                    |                                                                                                                                                 
cTGCACCCGGAAGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTCGTGTTCGACAACTGCCAgga                                                                                                                                                   >  4:1314/1‑149 (MQ=255)
    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                                                                               >  3:156577/1‑148 (MQ=255)
                                                                     tCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCa                                                                               >  6:368655/1‑148 (MQ=255)
                                                                               ccAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTg                                                                    >  2:234685/1‑149 (MQ=255)
                                                                                     tCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTg                                                              >  7:29412/1‑149 (MQ=255)
                                                                                     tCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTg                                                              >  6:208238/1‑149 (MQ=255)
                                                                                     tCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTg                                                              >  6:203211/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                                                                                    <  8:133148/112‑1 (MQ=255)
                                                                                                    tGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTc                                                                                    >  7:133148/1‑112 (MQ=255)
                                                                                                      gAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGTAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGa                                             <  5:368655/149‑1 (MQ=255)
                                                                                                           cggcggCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAg                                         <  7:338934/148‑1 (MQ=255)
                                                                                                           cggcggCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAg                                         <  6:53417/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                                    <  3:2067/140‑1 (MQ=255)
                                                                                                                        ccAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCAt                                    >  4:2067/1‑140 (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                       >  1:81912/3‑133 (MQ=255)
                                                                                                                                            cagCCAGGTCCAGGTCCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAg                       <  2:81912/131‑1 (MQ=255)
                                                                                                                                            cagCCAGGTCCAGGTCCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGt       <  5:122163/147‑1 (MQ=255)
                                                                                                                                                 aGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGTGGCTg  <  3:102011/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGTGGCTG  >  NC_002947/2957931‑2958224

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: ‑