Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,958,085 A→T 8.5% D282V (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,0850AT 8.5% 163.2 / 5.4 59D282V (GAC→GTC) PP_2589aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base T (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 = 3.62e-01

GACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGTGGCTGCCCGGTA  >  NC_002947/2957946‑2958231
                                                                                                                                           |                                                                                                                                                  
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                                                <  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                                           >  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)
                                                                                                                                          gacctggCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGTGGCTGCCCGGTa  <  8:80886/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
GACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGGCAAGAGCCCAAACCTGGTGTTCGACAACTGCCAGGACCTGGACCTGGCTGCGGAAAAGGCCGCCTTCGGCATTTTCTTCAACCAGGGCGAGGTGTGTTCGGCCAACTCGCGGCTGTACGTGCAGCGTGCCATCCACGATGAGTTCATCGAGCGCCTGCAGGCCAAGGCTCGCCAGTGGCTGCCCGGTA  >  NC_002947/2957946‑2958231

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