Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,657,538 T→A 7.6% D319V (GAC→GTC)  PP_1450 ← TPS family activation/secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,657,5380TA 7.6% 163.2 / 3.3 61D319V (GAC→GTC) PP_1450TPS family activation/secretion protein
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base A (2/3);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

CGTCGGGCTCTGCTGGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCGGGCCAACGA  >  NC_002947/1657393‑1657686
                                                                                                                                                 |                                                                                                                                                    
cGTCGGGCTCTGCGGGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCgg                                                                                                                                                   <  1:105310/149‑1 (MQ=255)
 gTCGGGCTCTGCTGGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCgg                                                                                                                                                   >  6:374170/1‑148 (MQ=255)
       ctctGCTGGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGACGGCAGcagc                                                                                                                                            <  4:187188/149‑4 (MQ=255)
         ctgctgGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGACGGCAGcagcag                                                                                                                                          <  7:240010/149‑6 (MQ=255)
            ctgGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTc                                                                                                                                         <  8:106115/147‑1 (MQ=255)
                aCGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGGTg                                                                                                                                     >  8:299141/1‑147 (MQ=255)
                aCGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTg                                                                                                                                     >  7:22801/1‑147 (MQ=255)
                  gATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTagg                                                                                                                                 >  7:289348/1‑149 (MQ=255)
                    tAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTaggg                                                                                                                                >  7:245633/1‑148 (MQ=255)
                       ccGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTaggggagg                                                                                                                            >  2:189042/1‑149 (MQ=255)
                       ccGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTaggggagg                                                                                                                            >  6:171973/1‑149 (MQ=255)
                        cgcgCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTaggggaggg                                                                                                                           <  8:245633/149‑1 (MQ=255)
                        cgcgCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTaggggaggg                                                                                                                           >  1:84635/1‑149 (MQ=255)
                           gcCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACATTACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAgggggg                                                                                                                        >  1:128223/1‑149 (MQ=255)
                             cAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCt                                                                                                                      <  2:242164/149‑1 (MQ=255)
                             cAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCt                                                                                                                      <  5:374170/149‑1 (MQ=255)
                               gcgcTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCtt                                                                                                                     >  7:222784/1‑148 (MQ=255)
                               gagcTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTg                                                                                                                    <  2:84635/147‑1 (MQ=255)
                                  cTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGAta                                                                                                                 <  2:1914/149‑1 (MQ=255)
                                       gTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGAtagctta                                                                                                            <  8:289348/149‑1 (MQ=255)
                                            tCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGa                                                                                                                   >  8:222466/1‑137 (MQ=255)
                                            tCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGa                                                                                                                   <  7:222466/137‑1 (MQ=255)
                                            tCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTAGCTGTAGGGGAGGGGGGCTTGATAGCTTAgctgg                                                                                                       >  4:96398/1‑149 (MQ=255)
                                               gctgctTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAgctggctg                                                                                                    <  2:128223/149‑1 (MQ=255)
                                                ctgctTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCt                                                                                                                      >  6:316333/1‑130 (MQ=255)
                                                ctgctTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCt                                                                                                                      <  5:316333/130‑1 (MQ=255)
                                                   ctTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGAt                                                                                                                  >  7:331438/1‑131 (MQ=255)
                                                   ctTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGAt                                                                                                                  <  8:331438/131‑1 (MQ=255)
                                                   ctTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGgcgc                                                                                                <  8:16151/149‑1 (MQ=255)
                                                       cGGTCCAGACGGGGGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGa                                                                                             <  2:276000/148‑1 (MQ=255)
                                                       cGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGa                                                                                             <  8:22801/148‑1 (MQ=255)
                                                       cGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGa                                                                                             <  8:222784/148‑1 (MQ=255)
                                                                      ctggcgctgAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGgc                                                                             >  4:186691/1‑149 (MQ=255)
                                                                         gcgctgAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGATGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGcc                                                                          <  2:13310/149‑1 (MQ=255)
                                                                                    gCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGc                                                               >  8:320107/1‑149 (MQ=255)
                                                                                      cTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTg                                                             >  4:163740/1‑149 (MQ=255)
                                                                                         gCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTgta                                                           <  1:189042/148‑1 (MQ=255)
                                                                                          cTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTgtag                                                          <  3:163740/148‑1 (MQ=255)
                                                                                           tGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTgtagta                                                        >  3:234183/1‑149 (MQ=255)
                                                                                           tGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTgtagt                                                         >  3:239152/1‑148 (MQ=255)
                                                                                                gcgcCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCGTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAAc                                                    <  1:389109/148‑1 (MQ=255)
                                                                                                          aCCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTg                                         <  2:53558/149‑1 (MQ=255)
                                                                                                           cccGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTgg                                        >  7:132034/1‑149 (MQ=255)
                                                                                                            ccGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGc                                       >  8:308655/1‑149 (MQ=255)
                                                                                                                  gACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTgg                                  <  7:320107/148‑1 (MQ=255)
                                                                                                                     acTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCGGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCgg                              >  6:69449/1‑149 (MQ=255)
                                                                                                                       tCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTACAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGCCCGGGc                            >  2:324779/1‑149 (MQ=255)
                                                                                                                              cTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATc                       <  1:55063/147‑1 (MQ=255)
                                                                                                                              cTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATc                       >  2:55063/1‑147 (MQ=255)
                                                                                                                                gATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAg                    >  2:154010/1‑148 (MQ=255)
                                                                                                                                   aGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCAcc                >  2:262452/1‑149 (MQ=255)
                                                                                                                                    gctgctgctGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCAcc                <  8:132034/148‑1 (MQ=255)
                                                                                                                                    gacaaggctgcCGACGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCg               >  8:278000/7‑149 (MQ=255)
                                                                                                                                     ctgctgctgCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGt              <  3:96398/149‑1 (MQ=255)
                                                                                                                                     acaaggctgcCGACGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGc                            >  7:170521/6‑135 (MQ=255)
                                                                                                                                     acaaggctgcCGACGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGc                            <  8:170521/130‑1 (MQ=255)
                                                                                                                                          gctgcCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCggg         <  5:69449/149‑1 (MQ=255)
                                                                                                                                           ctgcCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCgg          <  1:154010/147‑1 (MQ=255)
                                                                                                                                             gcCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCgg          >  8:30311/1‑145 (MQ=255)
                                                                                                                                             gcCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCgg          <  7:30311/145‑1 (MQ=255)
                                                                                                                                              ccGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCGGGCCaa     >  7:335721/1‑149 (MQ=255)
                                                                                                                                                 tCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCGGGCCAACGa  >  7:4611/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGTCGGGCTCTGCTGGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCGGGCCAACGA  >  NC_002947/1657393‑1657686

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