Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,447,029 A→T 7.7% K376M (AAG→ATG)  PP_1264 → fusaric acid resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,447,0290AT 7.7% 174.4 / 3.8 61K376M (AAG→ATG) PP_1264fusaric acid resistance protein
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base T (3/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGGGGGCTGCGGGCGGAGTACCTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGGCTGCTCAGCGATTGGCCCAGCGGCGCCATGATGAC  >  NC_002947/1446884‑1447176
                                                                                                                                                 |                                                                                                                                                   
gTGCGGGGGCTGCGGGCGGAGTACCTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGc                                                                                                                                                   <  8:260495/148‑1 (MQ=255)
         cTGCGGGCGGAGTACCTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGa                                                                                                                                         <  5:226962/149‑1 (MQ=255)
            cgggcggAGTACCTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCTCATGCTTGTGcgcggcca                                                                                                                                      <  2:300548/149‑5 (MQ=255)
             gggcggAGTACCTGGCAACGGACCCCAGCGTATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAg                                                                                                                                     >  1:31936/1‑149 (MQ=255)
              ggcggAGTACCTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAgcgc                                                                                                                                           <  1:289240/142‑1 (MQ=255)
              ggcggAGTACCTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAgcgc                                                                                                                                           >  2:289240/1‑142 (MQ=255)
                      aCCTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGaa                                                                                                                            >  6:272834/1‑149 (MQ=255)
                       ccTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGGGCAGTGGGACGAAc                                                                                                                           >  5:5180/1‑149 (MQ=255)
                        cTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAg                                                                                                                                              >  3:253408/1‑129 (MQ=255)
                        cTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAg                                                                                                                                              <  4:253408/129‑1 (MQ=255)
                                  aCCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGgcg                                                                                                                <  2:281939/149‑1 (MQ=255)
                                      cAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGa                                                                                                            >  6:295434/1‑149 (MQ=255)
                                      cAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGa                                                                                                            >  3:191483/1‑149 (MQ=255)
                                      cAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGa                                                                                                            >  2:21067/1‑149 (MQ=255)
                                      cAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGa                                                                                                            >  4:215562/1‑149 (MQ=255)
                                       aGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGAc                                                                                                           <  8:214291/149‑1 (MQ=255)
                                       aGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGAc                                                                                                           >  7:214291/1‑149 (MQ=255)
                                        gCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGGGCAGTGGGACGAACCCTACGTGGCGCAGAcc                                                                                                          <  6:5180/149‑1 (MQ=255)
                                                gACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAgctggct                                                                                                  >  8:178211/1‑149 (MQ=255)
                                                  ccTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGctggctgg                                                                                                >  5:248023/1‑149 (MQ=255)
                                                         gATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCt                                                                                         >  8:27903/1‑149 (MQ=255)
                                                           tttCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTgg                                                                                       <  7:27903/149‑1 (MQ=255)
                                                             tCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTggc                                                                                      <  6:248023/148‑1 (MQ=255)
                                                              ccATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTggcg                                                                                     <  7:178211/148‑1 (MQ=255)
                                                              ccATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTggcg                                                                                     <  2:162705/148‑1 (MQ=255)
                                                                   ccGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGgt                                                                               <  5:272834/149‑1 (MQ=255)
                                                                    cgcgTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGgtg                                                                              >  3:45459/1‑149 (MQ=255)
                                                                    cgcgTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGgtg                                                                              >  5:155632/1‑149 (MQ=255)
                                                                        ttCGATCTGCTGTACCGCTTCGTCGACGAGATGTTGTGCTACGCCGAATCCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTgcg                                                                          <  2:63134/149‑1 (MQ=255)
                                                                           gAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTgcgcgc                                                                       >  5:158283/1‑149 (MQ=255)
                                                                                gctgTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCgg                                                                  >  2:117236/1‑149 (MQ=255)
                                                                                gctgGACCGCTTCGTCGACGAGATGTTGCGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGTAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCgg                                                                  <  2:209231/149‑1 (MQ=255)
                                                                                    tACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCa                                                                <  1:117236/147‑1 (MQ=255)
                                                                                    tACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAg                                                               <  1:21067/148‑1 (MQ=255)
                                                                                                gacgaGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTg                                                             <  6:9927/138‑1 (MQ=255)
                                                                                                gacgaGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTg                                                             >  5:9927/1‑138 (MQ=255)
                                                                                                gacgaGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTgctgct                                                  <  5:295434/149‑1 (MQ=255)
                                                                                                   gagaTGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGCCACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTggg                                               >  5:218054/1‑149 (MQ=255)
                                                                                                    agaTGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGc                                              <  2:275285/149‑1 (MQ=255)
                                                                                                      aTGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAg                                            >  5:39521/1‑149 (MQ=255)
                                                                                                       tGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGc                                           <  6:155632/149‑1 (MQ=255)
                                                                                                        gTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGc                                           >  8:45434/1‑148 (MQ=255)
                                                                                                             gCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTgg                                     <  6:39521/149‑1 (MQ=255)
                                                                                                              cTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTAc                                        >  5:84578/1‑145 (MQ=255)
                                                                                                              cTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTAc                                        <  6:84578/145‑1 (MQ=255)
                                                                                                               tACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCgg                                                                  >  4:109735/1‑118 (MQ=255)
                                                                                                               tACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCgg                                                                  <  3:109735/118‑1 (MQ=255)
                                                                                                               tACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGgct                                   >  3:77041/1‑149 (MQ=255)
                                                                                                                aCGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGgctg                                  <  4:45459/149‑1 (MQ=255)
                                                                                                                                cTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTgg                                     <  1:287346/130‑1 (MQ=255)
                                                                                                                                cTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTgg                                     >  2:287346/1‑130 (MQ=255)
                                                                                                                                cTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGGCTGCTCAGCGATTGGc                     <  1:188070/146‑1 (MQ=255)
                                                                                                                                cTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGGCTGCTCAGCGATTGGc                     >  2:188070/1‑146 (MQ=255)
                                                                                                                                     tGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGTCAGCTACTGGCTGCTCAGCGATTGGCCCAgcggc             <  6:218054/149‑1 (MQ=255)
                                                                                                                                      gctcgcgCTCATGCTTGTGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAgctg                                                                                                     >  3:145607/4‑60 (MQ=25)
                                                                                                                                      gctcgcgCTCATGCTTGTGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAgctg                                                                                                     <  4:145607/57‑1 (MQ=25)
                                                                                                                                        tcgcgCTCATGCTTGTGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGGCTGCTCAGCGATTGGCCCAGCGgcgc           >  7:16036/2‑148 (MQ=255)
                                                                                                                                        tcgcgCTCATGCTTGTGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGAAGTGCTGTTGCTGCTGGGCAGCTACTGGCTGCTCAGCGATTGGCCCAGCGgcgc           >  3:176828/2‑148 (MQ=255)
                                                                                                                                        ccGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGGCTGCTCAGCGATTGGCCCAGCGGCGcc          <  6:158283/149‑1 (MQ=255)
                                                                                                                                          gcgcACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGGCTGCTCAGCGATTGGCCCAGCGGCGCCat        <  7:45434/149‑1 (MQ=255)
                                                                                                                                             cacaAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGGCTGCTCAGCGATTGGCCCAGCGGCGCCatgat     <  4:191483/149‑1 (MQ=255)
                                                                                                                                                aaGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGGCTGCTCAGCGATTGGCCCAGCGGCGCCATGATGAc  >  8:279294/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTGCGGGGGCTGCGGGCGGAGTACCTGGCAACGGACCCCAGCGAATCCGACCTGCTCGATTTCCATACCGCGTTCGAGCTGCTGTACCGCTTCGTCGACGAGATGTTGAGCTACGCCGAAACCCATGCCTCGCTGGCCGCGCACAAGCATGAGCGCGAGCAGTGGGACGAACCCTACGTGGCGCAGACCAGCTGGCTGGTATCGCTGGCGGCCGGTGTGCGCGCTTCGGCAGTGCTGTTGCTGCTGGGCAGCTACTGGCTGCTCAGCGATTGGCCCAGCGGCGCCATGATGAC  >  NC_002947/1446884‑1447176

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