Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,657,538 T→A 7.5% 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.5% 186.0 / 3.7 66D319V (GAC→GTC) PP_1450TPS family activation/secretion protein
Reads supporting (aligned to +/- strand):  ref base T (35/26);  new base A (2/3);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

GCGTCGGGCTCTGCTGGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCGGGCCAACGA  >  NC_002947/1657392‑1657686
                                                                                                                                                   |                                                                                                                                                       
gCGTCGGGCTCTGCTGGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCg                                                                                                                                                       >  3:208900/1‑149 (MQ=255)
        ctctGCTGGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCtt                                                                                                                                               >  5:117056/1‑149 (MQ=255)
           tgctgGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGACGGCAGcagca                                                                                                                                              <  6:182013/147‑5 (MQ=255)
             ctgGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGt                                                                                                                                             <  2:60025/146‑1 (MQ=255)
             ctgGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGACGGCAGcagcag                                                                                                                                             <  1:236089/146‑6 (MQ=255)
                gACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGt                                                                                                                                       <  4:208900/149‑1 (MQ=255)
                   gATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTagg                                                                                                                                    >  3:325456/1‑149 (MQ=255)
                      aaCCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTaggg                                                                                                                                   >  8:14883/1‑147 (MQ=255)
                        ccGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTaggggagg                                                                                                                               <  1:363170/149‑1 (MQ=255)
                         cgcgCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTaggggaggg                                                                                                                              >  4:327976/1‑149 (MQ=255)
                              cAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCt                                                                                                                         <  1:181476/149‑1 (MQ=255)
                               aGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCtt                                                                                                                        <  5:183733/149‑1 (MQ=255)
                               aGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCtt                                                                                                                        <  2:334653/149‑1 (MQ=255)
                                 cgcTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑Ga                                                                                                                   >  7:300034/1‑149 (MQ=255)
                                  gcTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GAt                                                                                                                  <  7:119710/149‑1 (MQ=255)
                                      cGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GAtagct                                                                                                              <  2:100758/149‑1 (MQ=255)
                                      cGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGAAGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GAtagct                                                                                                              >  1:289911/1‑149 (MQ=255)
                                       ggTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GAtagctt                                                                                                             <  8:202920/149‑1 (MQ=255)
                                             tCAGCTGCTTGCGGTCCAGACGGGCGCTGTCGCTGAGCATTCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAgctgg                                                                                                       <  8:113722/149‑1 (MQ=255)
                                             tCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAgctgg                                                                                                       >  5:226008/1‑149 (MQ=255)
                                                 ctgctTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑g                                                                                                                    >  3:98350/1‑132 (MQ=255)
                                                 ctgctTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑g                                                                                                                    <  4:98350/132‑1 (MQ=255)
                                                        cGGTCCAGACGGGCGCTGGCGCTGAGCATTCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGACTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAg                                                                                            <  3:15877/149‑1 (MQ=255)
                                                        cGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAg                                                                                            <  7:14883/149‑1 (MQ=255)
                                                        cGGTCCAGACGGGCGCTGGCGCTGAGCAGGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGTTAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGa                                                                                             <  2:289911/148‑1 (MQ=255)
                                                          gTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAggg                                                                                          <  8:300034/149‑1 (MQ=255)
                                                           tCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAggg                                                                                          >  8:348011/1‑148 (MQ=255)
                                                              aGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCa                                                                                        >  8:19253/1‑147 (MQ=255)
                                                               gACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAg                                                                                     <  6:117056/149‑1 (MQ=255)
                                                                      gctggcgctgAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAgg                                                                              >  3:175020/1‑149 (MQ=255)
                                                                          gcgctgAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTCGGACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTagggga                                                                                                                                 >  1:259726/1‑98 (MQ=255)
                                                                          gcgctgAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTagggga                                                                                                                                 <  2:259726/97‑1 (MQ=255)
                                                                             ctgAGCATGCCTTGATGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCCTCGGCAGCCTTGTCGCTGTAGGGGATCGGGGCTT‑‑‑GATAGCTTAGCTGGCAGGCGATGAGGGCAAAGGTCCAGGAGCCg                                                                         >  8:280054/1‑147 (MQ=255)
                                                                              tgAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAg                                                                               <  3:299219/140‑1 (MQ=255)
                                                                              tgAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAg                                                                               >  4:299219/1‑140 (MQ=255)
                                                                               gAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAgg                                                                     >  3:193924/1‑149 (MQ=255)
                                                                                     gCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGc                                                               >  4:94029/1‑149 (MQ=255)
                                                                                           cTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCACAGGTCCAGGCGCCGGAGGGGAGGCTgtag                                                          >  8:46711/1‑148 (MQ=255)
                                                                                                     cacaATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGatat                                               >  1:943/1‑149 (MQ=255)
                                                                                                          tACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GAtag                                                                                                                >  3:101308/1‑79 (MQ=255)
                                                                                                          tACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GAtag                                                                                                                <  4:101308/79‑1 (MQ=255)
                                                                                                            cccGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTg                                         <  3:94029/148‑1 (MQ=255)
                                                                                                             ccGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGc                                       <  6:226008/149‑1 (MQ=255)
                                                                                                              cGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCt                                      <  8:18289/149‑1 (MQ=255)
                                                                                                               gCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTTGAGGATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTgg                                        >  1:116077/1‑149 (MQ=255)
                                                                                                                 tC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTgtagta                                                        <  7:310123/128‑1 (MQ=255)
                                                                                                                 tC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTgtagta                                                        >  8:310123/1‑128 (MQ=255)
                                                                                                                       acTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCgg                              >  7:70709/1‑149 (MQ=255)
                                                                                                                          cAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCg                           >  7:196079/1‑149 (MQ=255)
                                                                                                                            ggCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTa                         >  5:201835/1‑149 (MQ=255)
                                                                                                                            ggCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTa                         >  5:196753/1‑149 (MQ=255)
                                                                                                                               ccTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCg                      <  7:19253/149‑1 (MQ=255)
                                                                                                                                cTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGACGGGCGTATCGa                     >  8:146640/1‑149 (MQ=255)
                                                                                                                                 tGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAg                    >  1:13978/1‑149 (MQ=255)
                                                                                                                                  gATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAg                    >  7:281727/1‑148 (MQ=255)
                                                                                                                                  gATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAg                    >  5:103797/1‑148 (MQ=255)
                                                                                                                                   ttAGCTGCGGCTGCCGTCGGCAGGCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCa                  <  7:280054/148‑1 (MQ=255)
                                                                                                                                   aTAGCTGCTGCTGCCGTCGGCAGCCTTCTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCa                  <  7:348011/149‑1 (MQ=255)
                                                                                                                                       ctgctgctgCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGt              >  7:340831/1‑149 (MQ=255)
                                                                                                                                       ctgctgctgACGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGGTT‑‑‑GATAGCTTAGCTGGCTTGAGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTATGTCCGGGCGTATCGATCACCGt              >  4:202813/1‑149 (MQ=255)
                                                                                                                                       acaaggctgcCGACGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAAc                                                    <  7:66533/106‑1 (MQ=255)
                                                                                                                                       acaaggctgcCGACGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAAc                                                    >  8:66533/6‑111 (MQ=255)
                                                                                                                                       acaaggctgcCGACGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGt              >  2:222246/6‑149 (MQ=255)
                                                                                                                                                ccGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCGGGCCaa     >  4:331798/1‑149 (MQ=255)
                                                                                                                                                ccGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCGGGCCaa     >  3:245927/1‑149 (MQ=255)
                                                                                                                                                  gTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCGGGCCAAc    >  7:4514/1‑148 (MQ=255)
                                                                                                                                                   tCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCGGGCCAACGa  <  8:70709/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
GCGTCGGGCTCTGCTGGACGATAACCGCGCCAGCGCTGCGGTTGATCAGCTGCTTGCGGTCCAGACGGGCGCTGGCGCTGAGCATGCCTTGCTGGTTGCGCCACAATACCCGCTC‑GACACTCAGGCCCTGATAGCTGCTGCTGCCGTCGGCAGCCTTGTCGCTGTAGGGGAGGGGGGCTT‑‑‑GATAGCTTAGCTGGCTGGCGCTGAGGGCAAAGGTCCAGGCGCCGTAGGGGATGCTGTAGTACAACGATATCCCTTGGCTTTGGCCGGGCGTATCGAGCACCGTGCGGGCCAACGA  >  NC_002947/1657392‑1657686

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