Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,262,496 T→A 7.0% T223T (ACT→ACA PP_4638 → methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,262,4960TA 7.0% 209.6 / 5.4 71T223T (ACT→ACAPP_4638methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (31/35);  new base A (3/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

CAGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAGGGGTTC  >  NC_002947/5262356‑5262641
                                                                                                                                            |                                                                                                                                                  
cAGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGcc                                                                                                                                            <  6:232830/149‑1 (MQ=255)
 aGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGATCACGTGGCCGGCGTGACTGTCACGcca                                                                                                                                           <  1:12471/149‑1 (MQ=255)
               gccgccAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACAt                                                                                                                             >  3:156796/1‑149 (MQ=255)
                 cgccAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGc                                                                                                                           >  4:123101/1‑149 (MQ=255)
                  gccAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCt                                                                                                                          >  8:348990/1‑149 (MQ=255)
                   ccAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTc                                                                                                                         >  3:191532/1‑149 (MQ=255)
                     aaGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCAc                                                                                                                                               <  1:387574/125‑1 (MQ=255)
                     aaGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCAc                                                                                                                                               >  2:387574/1‑125 (MQ=255)
                     aaGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGAc                                                                                                                               >  1:406241/1‑141 (MQ=255)
                     aaGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGAc                                                                                                                               <  2:406241/141‑1 (MQ=255)
                         aTCAGGAAGCCATACACTGGATGCGCCGTGAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCg                                                                                                                  <  3:342947/149‑1 (MQ=255)
                         aTCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCg                                                                                                                  <  7:348990/149‑1 (MQ=255)
                          tCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCTGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCgg                                                                                                                 <  2:386501/149‑1 (MQ=255)
                              gAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGAc                                                                                                             >  6:163698/1‑149 (MQ=255)
                                    aTACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGgcg                                                                                                         >  4:153166/1‑147 (MQ=255)
                                    aTACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCa                                                                                                       <  2:140665/149‑1 (MQ=255)
                                        acTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCAGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGAGGAGACGGCGCAGTc                                                                                                    >  1:139144/1‑148 (MQ=255)
                                         cTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCa                                                                                                  <  4:156796/149‑1 (MQ=255)
                                          tGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCt                                                                                                                          <  2:166632/125‑1 (MQ=255)
                                          tGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCt                                                                                                                          >  1:166632/1‑125 (MQ=255)
                                          tGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTcc                                                                                                   >  4:93798/1‑147 (MQ=255)
                                            gATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAg                                                                                                <  8:209823/148‑1 (MQ=255)
                                            gATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCATCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGGGACATGCTCGAAGCGGAGACGGCGCAGTCCAAg                                                                                                <  4:49003/148‑1 (MQ=255)
                                                cgcCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCt                                                                                            <  3:93798/148‑1 (MQ=255)
                                                 gcCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCACACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCttt                                                                                          >  8:165251/1‑149 (MQ=255)
                                                       aGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAg                                                                                    >  5:402819/1‑149 (MQ=255)
                                                       aGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATACATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAg                                                                                    >  8:366886/1‑149 (MQ=255)
                                                          ccTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTg                                                                                         <  7:58144/141‑1 (MQ=255)
                                                          ccTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTg                                                                                         >  8:58144/1‑141 (MQ=255)
                                                           ttCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACAGTCACGccgg                                                                                                                                          >  1:93551/2‑90 (MQ=255)
                                                           ttCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACAGTCACGccgg                                                                                                                                          <  2:93551/91‑3 (MQ=255)
                                                              aCCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCg                                                                                                                  >  7:317289/1‑112 (MQ=255)
                                                              aCCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCg                                                                                                                  <  8:317289/112‑1 (MQ=255)
                                                              aCCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCa                                                                             <  2:298240/149‑1 (MQ=255)
                                                                   aGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGtt                                                                        <  4:191532/149‑1 (MQ=255)
                                                                    gCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCTCGCAGTGtt                                                                        >  3:178647/1‑148 (MQ=255)
                                                                           cTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCt                                                                <  6:402819/149‑1 (MQ=255)
                                                                              gcctgcctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTg                                                               <  5:21408/147‑1 (MQ=255)
                                                                                ctgcctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCt                                                           >  2:111018/1‑149 (MQ=255)
                                                                                   cctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTg                                                                                         >  8:268188/1‑116 (MQ=255)
                                                                                   cctCATGAGCCTCACCCACGGGCGCGCGACGATGCGCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTg                                                                                         <  7:268188/116‑1 (MQ=255)
                                                                                     tCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTg                                                               >  4:232923/1‑140 (MQ=255)
                                                                                     tCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTg                                                               <  3:232923/140‑1 (MQ=255)
                                                                                     tCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTg                                                         >  4:241558/1‑146 (MQ=255)
                                                                                         gAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTa                                                  >  5:190123/1‑149 (MQ=255)
                                                                                           gCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATc                                                >  1:296001/1‑149 (MQ=255)
                                                                                              tCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGtt                                                                        <  8:267456/122‑1 (MQ=255)
                                                                                              tCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGtt                                                                        >  7:267456/1‑122 (MQ=255)
                                                                                                aCCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGt                                             <  2:296001/147‑1 (MQ=255)
                                                                                                aCCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTa                                            <  6:190123/148‑1 (MQ=255)
                                                                                                aCCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTAc                                           <  3:153166/149‑1 (MQ=255)
                                                                                                   cACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACaaa                                        <  1:380696/149‑1 (MQ=255)
                                                                                                     cGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGc                                      <  2:139144/149‑1 (MQ=255)
                                                                                                         cgcgcgACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGAt                                  <  5:163698/149‑1 (MQ=255)
                                                                                                          gcgcgACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATg                                 <  1:111018/149‑1 (MQ=255)
                                                                                                                    gCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTg                                                                       >  4:57519/1‑101 (MQ=255)
                                                                                                                    gCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTg                                                                       <  3:57519/101‑1 (MQ=255)
                                                                                                                        gACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGgt                     <  7:165251/147‑1 (MQ=255)
                                                                                                                        gACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGgt                     <  7:366886/147‑1 (MQ=255)
                                                                                                                         aCCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAg                                                                                                <  7:30181/71‑1 (MQ=255)
                                                                                                                         aCCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAg                                                                                                >  8:30181/1‑71 (MQ=255)
                                                                                                                            aCGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCAc               >  1:323939/1‑149 (MQ=255)
                                                                                                                            aCGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCAc               >  8:244876/1‑149 (MQ=255)
                                                                                                                                 gatggCGTGACAGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTg                                    >  4:403867/4‑123 (MQ=255)
                                                                                                                                 gatggCGTGACAGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTg                                    <  3:403867/120‑1 (MQ=255)
                                                                                                                                 gatggCGTGACAGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGt           >  3:154387/4‑148 (MQ=255)
                                                                                                                                 gCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGt           <  3:241558/148‑1 (MQ=255)
                                                                                                                                   cGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGAAGATTATCGGTACAAAGCTGATGGGCTACTCCGGTGTGCACCTGTCAg        >  3:293251/1‑149 (MQ=255)
                                                                                                                                      cGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAgggg     >  3:327874/1‑149 (MQ=255)
                                                                                                                                         gACTGTCACGCCATCCATGCGCGACATGCTCGATAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAGGGGtt   >  1:55459/1‑149 (MQ=255)
                                                                                                                                          aCTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCTCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAGGGGTTc  <  4:178647/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CAGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAGGGGTTC  >  NC_002947/5262356‑5262641

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: