Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,849,921 T→A 6.4% L209Q (CTG→CAG)  ilvD → dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,849,9210TA 6.4% 183.8 / 2.5 62L209Q (CTG→CAG) ilvDdihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base A (2/2);  total (33/29)
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.80e-01

GACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTCCTGCCACGCAGCATCG  >  NC_002947/5849776‑5850067
                                                                                                                                                 |                                                                                                                                                  
gACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCTCCTGCCCGACCTGGGTTTCGTGCTCCGGCATGTTCACCGCCAACTCGTTGAACTGCCTGTCCGTTGCACTgg                                                                                                                                                  <  4:135047/148‑1 (MQ=255)
gACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTggg                                                                                                                                                 <  4:382237/149‑1 (MQ=255)
gACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTgg                                                                                                                                                  <  3:172900/148‑1 (MQ=255)
gACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTg                                                                                                                                                   <  1:241157/147‑1 (MQ=255)
        cGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGccc                                                                                                                                         >  1:98163/1‑149 (MQ=255)
           tGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGc                                                                                                                                      >  1:136499/1‑149 (MQ=255)
              cATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCgg                                                                                                                                   <  3:269436/149‑1 (MQ=255)
              cATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCgg                                                                                                                                   <  2:373253/149‑1 (MQ=255)
               aTGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCggg                                                                                                                                  >  5:320551/1‑149 (MQ=255)
               aTGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCggg                                                                                                                                  <  4:369440/149‑1 (MQ=255)
               aTGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCggg                                                                                                                                  >  1:189707/1‑149 (MQ=255)
                tGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCggg                                                                                                                                  >  8:100031/1‑148 (MQ=255)
                     aTCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACgg                                                                                                                            <  7:131445/149‑1 (MQ=255)
                          ggCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGa                                                                                                                       <  3:253019/149‑1 (MQ=255)
                            ccGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGAcc                                                                                                                     <  6:320551/149‑1 (MQ=255)
                                   cACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGcc                                                                                                               <  5:353331/148‑1 (MQ=255)
                                                     aGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGag                                                                                             >  2:340710/1‑148 (MQ=255)
                                                           cGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGt                                                                                      <  2:98163/149‑1 (MQ=255)
                                                           cGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGt                                                                                      >  3:107424/1‑149 (MQ=255)
                                                             aaTACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTc                                                                                    <  2:136499/149‑1 (MQ=255)
                                                                 cGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTc                                                                                 <  7:100031/148‑1 (MQ=255)
                                                                  gAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCAc                                                                               >  6:6358/1‑149 (MQ=255)
                                                                   agcgcagcgcCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCAcg                                                                              >  2:22990/1‑149 (MQ=255)
                                                                    gcgcagcgcCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCAGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCAcgg                                                                             >  6:30389/1‑149 (MQ=255)
                                                                    gcgcagcgcCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCAGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCAcgg                                                                             >  6:24967/1‑149 (MQ=255)
                                                                     cgcagcgcCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCAcggc                                                                            >  5:25128/1‑149 (MQ=255)
                                                                           gcCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGgccggccg                                                                      <  2:189707/149‑1 (MQ=255)
                                                                             cTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTa                                                                    <  6:25128/149‑1 (MQ=255)
                                                                                cccGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCa                                                                 <  4:107424/149‑1 (MQ=255)
                                                                                  cGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCACCTCGATGAACTGCCTGACCGAAGCCCTGCGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACAGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCAtc                                                               >  6:293876/1‑149 (MQ=255)
                                                                                  cGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCAtc                                                               >  5:42618/1‑149 (MQ=255)
                                                                                    aCCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCAtcgt                                                             >  4:162377/1‑149 (MQ=255)
                                                                                      cTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCAtcgtc                                                            >  8:271844/1‑148 (MQ=255)
                                                                                      cTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCAtcgtc                                                            >  6:316368/1‑148 (MQ=255)
                                                                                      cTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCAtcgtc                                                            <  7:271844/148‑1 (MQ=255)
                                                                                           gTTCGTGCTCCGGCATGTTCACCGCCAAGTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCtg                                                      >  1:221691/1‑149 (MQ=255)
                                                                                            ttCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCtgt                                                     <  1:22990/149‑1 (MQ=255)
                                                                                            ttCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCtgt                                                     <  6:42618/149‑1 (MQ=255)
                                                                                                       gCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGctac                                          >  2:119143/1‑149 (MQ=255)
                                                                                                       gCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGctac                                          >  8:320617/1‑149 (MQ=255)
                                                                                                            ttCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGACCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACgg                                     >  7:11288/1‑149 (MQ=255)
                                                                                                               aCCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCAtcgt                                                             <  7:313325/122‑1 (MQ=255)
                                                                                                               aCCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCAtcgt                                                             >  8:313325/1‑122 (MQ=255)
                                                                                                                   ccAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAAt                              <  1:340710/149‑1 (MQ=255)
                                                                                                                     ttggcggtGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACACACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGa                            <  3:238832/142‑1 (MQ=255)
                                                                                                                      aCTCGATGAACTGCCTGACCGAAGCCCAGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGAc                           <  5:30389/149‑1 (MQ=255)
                                                                                                                      aCTCGATGAACTGCCTGACCGAAGCCCAGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGAc                           <  5:24967/149‑1 (MQ=255)
                                                                                                                        tCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGa                         >  1:153108/1‑149 (MQ=255)
                                                                                                                             gAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGACTCgg                    >  5:296731/1‑149 (MQ=255)
                                                                                                                              aaCTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGt                   >  7:87961/1‑149 (MQ=255)
                                                                                                                               aCTTCCTTACCGAAGCCCTTGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGACACCCACTCCGACCTCTAGCAGCTGTTCCTCACGTCCGGCCGTACCATCGTCGAGCTCTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTc                  <  7:150441/149‑1 (MQ=255)
                                                                                                                               aCTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTc                  >  8:17924/1‑149 (MQ=255)
                                                                                                                               aCTGCCTGACCGAAGCCCTGGACCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTc                  <  8:11288/149‑1 (MQ=255)
                                                                                                                                  gCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGCCCCTGGCCCCCCACTCCGACCGCGAGCAGCTGGTCCTCAAGGCCGGCCGTAACATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGCATGACGAGTCGGTCCTg               >  2:328792/1‑149 (MQ=255)
                                                                                                                                   ccTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTCCTGc              >  7:196755/1‑149 (MQ=255)
                                                                                                                                   ccTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTCCTGc              >  8:294457/1‑149 (MQ=255)
                                                                                                                                        ccGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTCCTGCCACgc         <  7:320617/149‑1 (MQ=255)
                                                                                                                                         cGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTCCTGCCACg          >  2:58555/1‑147 (MQ=255)
                                                                                                                                           aaGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTTGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTCCTGCCACgcag       >  7:258494/1‑148 (MQ=255)
                                                                                                                                              cccTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCCCGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTCCTGCCACGCAGCATc   >  4:202204/1‑149 (MQ=255)
                                                                                                                                               ccTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTCCTGCCACGCAGCATCg  <  5:6358/149‑1 (MQ=255)
                                                                                                                                               ccTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTCCTGCCACGCAGCATCg  >  1:322091/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGAGCGCAGCGCCTGCCCGACCTGTGGTTCGTGCTCCGGCATGTTCACCGCCAACTCGATGAACTGCCTGACCGAAGCCCTGGGCCTGGCCCTGCCGGGCAACGGTTCGACCCTGGCCACCCACTCCGACCGCGAGCAGCTGTTCCTCACGGCCGGCCGTACCATCGTCGAGCTGTGCAAGCGCTACTACGGCGAGAATGACGAGTCGGTCCTGCCACGCAGCATCG  >  NC_002947/5849776‑5850067

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