Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,644,062 A→T 5.6% N310I (AAC→ATC)  PP_4108 → 4‑aminobutyrate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,644,0620AT 5.6% 216.5 / 2.9 71N310I (AAC→ATC) PP_41084‑aminobutyrate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base T (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

ATCGCTGGCGGTATGCCGCTGGGTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAATGCCGATGGCAGCCC  >  NC_002947/4643914‑4644209
                                                                                                                                                    |                                                                                                                                                   
aTCGCTGGCGGTATGCCGCTGGGTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGaaaa                                                                                                                                                     >  3:343844/1‑149 (MQ=255)
aTCGCTGGCGGTATGCCGCTGGGTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGaaaa                                                                                                                                                     <  7:377290/149‑1 (MQ=255)
     tGGCGGTATGCCGCTGGGTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTg                                                                                                                                                >  6:243082/1‑149 (MQ=255)
     tGGCGGTATGCCGCTGGGTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTg                                                                                                                                                >  7:1722/1‑149 (MQ=255)
       gCGGTATGCCGCTGGGTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTg                                                                                                                                                <  5:38627/147‑1 (MQ=255)
       gCGGTATGCCGCTGGGTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGc                                                                                                                                               <  8:33806/148‑1 (MQ=255)
                    gggTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTggg                                                                                                                                       <  8:256141/143‑1 (MQ=255)
                    gggTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTggg                                                                                                                                       >  7:256141/1‑143 (MQ=255)
                      gTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATAACCGACGAAAACCTTGCAACCTGGGGTGAACgg                                                                                                                               <  5:405616/149‑1 (MQ=255)
                       tGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGc                                                                                                                              <  3:352650/149‑1 (MQ=255)
                            tGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGaa                                                                                                                         >  6:410768/1‑149 (MQ=255)
                               tCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAg                                                                                                                      <  2:99198/149‑1 (MQ=255)
                                   gcgcAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCa                                                                                                                  <  5:243082/149‑1 (MQ=255)
                                    cgcAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCAt                                                                                                                 <  5:367087/149‑1 (MQ=255)
                                     gcAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATc                                                                                                                >  7:295615/1‑149 (MQ=255)
                                      cAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCg                                                                                                               >  3:303833/1‑149 (MQ=255)
                                       aaGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGt                                                                                                              >  3:140557/1‑149 (MQ=255)
                                         ggAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCa                                                                                                            <  8:1722/149‑1 (MQ=255)
                                           aaTTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGGAAGCGGGGGTGAACGGCCGGAACAGGCCATCGTCAGc                                                                                                          >  5:399496/1‑149 (MQ=255)
                                            aTTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAATCCTTGCAACCTg                                                                                                                                         <  3:300874/117‑1 (MQ=255)
                                            aTTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAATCCTTGCAACCTg                                                                                                                                         >  4:300874/1‑117 (MQ=255)
                                            aTTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGcc                                                                                                         <  4:140557/149‑1 (MQ=255)
                                               gATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGc                                                                                                       <  2:135541/148‑1 (MQ=255)
                                               gATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGc                                                                                                       <  8:39604/148‑1 (MQ=255)
                                                 tGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCGCCTATTCGGGCAACCCGATATCCAGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGCCCGACGAAAACCTTGCAACCTCGGGTGAACGGCACGAACAGGCCATCGTCAGCCGCTAt                                                                                                    >  3:236765/1‑149 (MQ=255)
                                                  ggcggcGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGcc                                                                                                         >  3:70827/1‑143 (MQ=255)
                                                  ggcggcGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGcc                                                                                                         <  4:70827/143‑1 (MQ=255)
                                                            cccAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTgg                                                                                           >  7:145158/1‑147 (MQ=255)
                                                            cccAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTgg                                                                                           >  7:262150/1‑147 (MQ=255)
                                                            cccAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGa                                                                                          <  4:209674/148‑1 (MQ=255)
                                                              cAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAATCCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAgg                                                                                       >  3:22279/1‑149 (MQ=255)
                                                              cAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAgg                                                                                       >  6:108355/1‑149 (MQ=255)
                                                                  ggcggcCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTc                                                                                   <  4:343844/149‑1 (MQ=255)
                                                                    cggcCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCgg                                                                                 <  4:247590/149‑1 (MQ=255)
                                                                        cTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCt                                                                             >  4:329614/1‑149 (MQ=255)
                                                                         tacaaggCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCa                                                                                                                             >  6:396551/6‑100 (MQ=255)
                                                                         tacaaggCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCa                                                                                                                             <  5:396551/95‑1 (MQ=255)
                                                                           ggTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGa                                                                           <  4:303833/148‑1 (MQ=255)
                                                                            gTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGc                                                                         >  3:257877/1‑149 (MQ=255)
                                                                               gCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCg                                                                      <  3:162929/149‑1 (MQ=255)
                                                                                 aCCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCgg                                                                                 <  3:351235/136‑1 (MQ=255)
                                                                                 aCCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCgg                                                                                 >  4:351235/1‑136 (MQ=255)
                                                                                         gggCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAATCCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGgcc                                                            <  4:22279/149‑1 (MQ=255)
                                                                                         gggCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGgcc                                                            >  2:380519/1‑149 (MQ=255)
                                                                                          ggCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTgg                                                                                           >  3:182377/1‑117 (MQ=255)
                                                                                          ggCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTgg                                                                                           <  4:182377/117‑1 (MQ=255)
                                                                                           gCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGgccgc                                                          <  4:257877/149‑1 (MQ=255)
                                                                                                ccGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCt                                                                             >  6:253942/1‑125 (MQ=255)
                                                                                                ccGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCt                                                                             <  5:253942/125‑1 (MQ=255)
                                                                                                 cGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGAcc                                                    >  3:289000/1‑149 (MQ=255)
                                                                                                 cGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGAcc                                                    <  1:201585/149‑1 (MQ=255)
                                                                                                   atatCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCgg                                                  >  1:343271/1‑149 (MQ=255)
                                                                                                      tCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGt                                               >  3:289389/1‑149 (MQ=255)
                                                                                                        cTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTg                                              <  5:410768/148‑1 (MQ=255)
                                                                                                         tgtgCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTggg                                            >  6:145308/1‑149 (MQ=255)
                                                                                                             cggcggCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGcc                                        <  4:406527/149‑1 (MQ=255)
                                                                                                               gcggcCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCAt                                      <  6:116296/149‑1 (MQ=255)
                                                                                                                 ggcCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGc                                    <  3:329614/149‑1 (MQ=255)
                                                                                                                   ccTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGa                                  <  4:289389/149‑1 (MQ=255)
                                                                                                                    cTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAg                                 <  8:302709/149‑1 (MQ=255)
                                                                                                                    cTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAg                                 <  1:380519/149‑1 (MQ=255)
                                                                                                                        gCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGa                                  <  3:349787/144‑1 (MQ=255)
                                                                                                                        gCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGa                                  >  4:349787/1‑144 (MQ=255)
                                                                                                                          cAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATc                            >  7:239907/1‑148 (MQ=255)
                                                                                                                                   gcAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGcc                   >  4:232870/1‑148 (MQ=255)
                                                                                                                                     aGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAAt                >  5:101051/1‑149 (MQ=255)
                                                                                                                                     aGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAAt                <  3:232870/149‑1 (MQ=255)
                                                                                                                                           ccGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAATGCCGAt          <  2:343271/149‑1 (MQ=255)
                                                                                                                                           ccGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAATGCCGAt          >  4:165151/1‑149 (MQ=255)
                                                                                                                                           ccGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAATGCCGAt          <  4:289000/149‑1 (MQ=255)
                                                                                                                                                gaAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAATGCCGATGGCAg     >  7:387890/1‑149 (MQ=255)
                                                                                                                                                gaAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAATGCCGATGGCAg     >  7:82461/1‑149 (MQ=255)
                                                                                                                                                gaAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAATGCCGATGGCAg     >  7:221307/1‑149 (MQ=255)
                                                                                                                                                   aaCCTTGCAACCTGGGGTGAACGGCAGGTACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAATGCCGATGGCAGccc  >  1:228178/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCGCTGGCGGTATGCCGCTGGGTGCAGTGGTCGGGCGCAAGGAATTGATGGCGGCGCTGCCCAAGGGCGGCCTGGGTGGCACCTATTCGGGCAACCCGATATCCTGTGCGGCGGCCTTGGCCAGCCTGGCGCAGATGACCGACGAAAACCTTGCAACCTGGGGTGAACGGCAGGAACAGGCCATCGTCAGCCGCTATGAGCGCTGGAAGGCTTCGGGCCTGAGCCCGTACATCGGCCGCCTGACCGGCGTGGGCGCCATGCGAGGGATCGAGTTCGCCAATGCCGATGGCAGCCC  >  NC_002947/4643914‑4644209

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: