Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,761,303 A→G 9.2% T137A (ACC→GCC)  PP_1567 → capsid protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,761,3030AG 9.2% 212.9 / 9.7 76T137A (ACC→GCC) PP_1567capsid protein
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base G (3/4);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

GCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTCCAGA  >  NC_002947/1761154‑1761450
                                                                                                                                                     |                                                                                                                                                   
tcTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACcaccggca                                                                                                                                                     <  1:23269/149‑1 (MQ=255)
         ccAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGACGGAAGGCCGTCACCAACATTACTCGCCCTACCACCGGCACCGGTGGCg                                                                                                                                            >  2:106883/1‑150 (MQ=255)
          cAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGTGCCGGTggtag                                                                                                                                           <  2:1210560/150‑4 (MQ=255)
            aTGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCg                                                                                                                                         <  1:1204267/150‑1 (MQ=255)
             tGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGTGCCGGTggtaga                                                                                                                                          >  1:1396113/1‑144 (MQ=255)
             tGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGTGCCGGTggtaga                                                                                                                                          <  2:1396113/148‑5 (MQ=255)
                tCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGt                                                                                                                                                >  1:118396/1‑139 (MQ=255)
                tCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGt                                                                                                                                                <  2:118396/139‑1 (MQ=255)
                  cAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAg                                                                                                                                   >  2:1275785/1‑150 (MQ=255)
                     ttGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCa                                                                                                                                >  2:615016/1‑150 (MQ=255)
                     ttGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGACATCAATATTACCAGCTCTACCACCGGCACCCGTGGCGTCGGCGTAGCCa                                                                                                                                <  1:1159910/150‑1 (MQ=255)
                         cTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTa                                                                                                                                    >  1:283259/1‑142 (MQ=255)
                         cTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTa                                                                                                                                    <  2:283259/142‑1 (MQ=255)
                          ttCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAg                                                                                                                           >  2:904669/1‑150 (MQ=255)
                          ttCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAg                                                                                                                           >  2:1986204/1‑150 (MQ=255)
                                  ggCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCcg                                                                                                                   >  2:1681909/1‑150 (MQ=255)
                                  ggCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCcg                                                                                                                   >  2:698694/1‑150 (MQ=255)
                                  ggCTGGCGAAGAACTGACAGAGGGTGACGGCTTCAAGTCGCTGCAAACTCGCGGCCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCcg                                                                                                                   >  2:1918244/1‑150 (MQ=255)
                                     tGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGt                                                                                                                >  1:1311774/1‑150 (MQ=255)
                                     tGGCGAAGAACGGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGt                                                                                                                <  2:1577158/150‑1 (MQ=255)
                                       gCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTAGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGCTGGCGTCAGCGTAGCCATCCAGCCGACCCGCGTTc                                                                                                              >  1:308417/1‑150 (MQ=255)
                                         gaagaaCTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTccc                                                                                                            >  1:178037/1‑150 (MQ=255)
                                         gaagaaCTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTccc                                                                                                            <  2:1311774/150‑1 (MQ=255)
                                               cTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATc                                                                                                      >  1:2066374/1‑150 (MQ=255)
                                               cTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATc                                                                                                      <  1:1072481/150‑1 (MQ=255)
                                                  aCCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACcaccggcaccg                                                                                                                                                  <  2:2365489/103‑1 (MQ=255)
                                                  aCCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACcaccggcaccg                                                                                                                                                  >  1:2365489/1‑103 (MQ=255)
                                                  aCCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCCTCGTg                                                                                                   >  2:994192/1‑150 (MQ=255)
                                                  aCCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTg                                                                                                   >  1:2304350/1‑150 (MQ=255)
                                                  aCCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTg                                                                                                   <  1:2188544/150‑1 (MQ=255)
                                                       gtgtGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTgg                                                                                              <  1:698694/150‑1 (MQ=255)
                                                         gtgACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGc                                                                                            >  2:471683/1‑150 (MQ=255)
                                                          tgACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGcc                                                                                           >  1:1172995/1‑150 (MQ=255)
                                                          tgACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGATCGAGTCAGCAGCGGAGCCGTCACCAACAGTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGcc                                                                                           >  1:1500970/1‑150 (MQ=255)
                                                            acgacTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCgg                                                                                         >  2:2290817/1‑150 (MQ=255)
                                                                 ttCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAgcgc                                                                                    <  1:1153962/150‑1 (MQ=255)
                                                                 ttCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAgcgc                                                                                    <  1:904669/150‑1 (MQ=255)
                                                                   cAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGa                                                                                  <  2:2066374/150‑1 (MQ=255)
                                                                   cAAGTCGCTGCAAAATTGCGGGCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGa                                                                                  <  2:1640828/150‑1 (MQ=255)
                                                                      gTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGtt                                                                               <  2:648162/150‑1 (MQ=255)
                                                                      gTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGtt                                                                               <  1:129294/150‑1 (MQ=255)
                                                                       tCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGTGCCGGTggta                                                                                                                                            <  1:134734/88‑3 (MQ=255)
                                                                       tCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGTGCCGGTggta                                                                                                                                            >  2:134734/1‑86 (MQ=255)
                                                                              aaaCTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTc                                                                       <  1:2038034/150‑1 (MQ=255)
                                                                                            ggCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATc                                                                                                                              <  1:2108959/81‑1 (MQ=255)
                                                                                            ggCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATc                                                                                                                              >  2:2108959/1‑81 (MQ=255)
                                                                                            ggCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATg                                                         >  2:1641020/1‑150 (MQ=255)
                                                                                             gCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGc                                                        <  1:994192/150‑1 (MQ=255)
                                                                                                  gTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGc                                                   >  1:595363/1‑150 (MQ=255)
                                                                                                   ttCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGTGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGcc                                                  >  2:568514/1‑150 (MQ=255)
                                                                                                   ttCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCAt                                                                                                                               >  1:1580480/1‑73 (MQ=255)
                                                                                                   ttCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCAt                                                                                                                               <  2:1580480/73‑1 (MQ=255)
                                                                                                   ttCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGa                                                                                        >  1:1679547/1‑112 (MQ=255)
                                                                                                   ttCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGa                                                                                        <  2:1679547/112‑1 (MQ=255)
                                                                                                   ttCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGcc                                                  >  2:853250/1‑150 (MQ=255)
                                                                                                        gTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACt                                             <  2:456816/150‑1 (MQ=255)
                                                                                                         tAGGCCGGAAGGCCGTCACAAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTg                                            <  2:2199529/150‑1 (MQ=255)
                                                                                                           ggCTGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGc                                          <  1:470590/150‑1 (MQ=255)
                                                                                                           ggCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGc                                          >  2:218289/1‑150 (MQ=255)
                                                                                                           ggCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTAGCGACCTGATCATGCCAGGCCGTACTGGc                                          >  2:1161518/1‑150 (MQ=255)
                                                                                                            gCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCGAACCCGCGTTCCCGGCATAGTAGTTGGTCCGGAGCGCGAGTTCTCGATTGGCGACCTGATCATGCCAGGCCGTACTGGCt                                         >  2:1709475/1‑150 (MQ=255)
                                                                                                            gCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCt                                         >  1:2196744/1‑150 (MQ=255)
                                                                                                            gCCGGAAGGCCGTCACCAACATTACAAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCt                                         >  2:1411531/1‑150 (MQ=255)
                                                                                                                 aaGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAAc                                    <  1:106883/150‑1 (MQ=255)
                                                                                                                     ccGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAg                                <  1:1681909/150‑1 (MQ=255)
                                                                                                                     ccGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAg                                >  2:1133558/1‑150 (MQ=255)
                                                                                                                       gTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTc                              >  1:2199706/1‑150 (MQ=255)
                                                                                                                        tCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGc                                                                                                                          <  2:1695258/57‑1 (MQ=255)
                                                                                                                        tCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGc                                                                                                                          >  1:1695258/1‑57 (MQ=255)
                                                                                                                            cAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGtt                         >  1:511618/1‑150 (MQ=255)
                                                                                                                             aaCATTACCAGCTCTACCACAGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTc                        <  2:2199706/150‑1 (MQ=255)
                                                                                                                                 ttACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGc                    >  2:2110846/1‑150 (MQ=255)
                                                                                                                                   aCCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAg                  <  1:1709475/150‑1 (MQ=255)
                                                                                                                                   aCCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAg                  <  1:1986204/150‑1 (MQ=255)
                                                                                                                                    ccAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAg                                                                                 <  2:327306/86‑1 (MQ=255)
                                                                                                                                    ccAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAg                                                                                 >  1:327306/1‑86 (MQ=255)
                                                                                                                                         tctACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGAcc            <  1:853250/150‑1 (MQ=255)
                                                                                                                                         gacgccaccGGTGCCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCAt                                                                         <  2:2140863/85‑1 (MQ=255)
                                                                                                                                         gacgccaccGGTGCCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCAt                                                                         >  1:2140863/5‑89 (MQ=255)
                                                                                                                                              caccggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGaa                                     <  1:952148/120‑1 (MQ=255)
                                                                                                                                              caccggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGaa                                     >  2:952148/1‑120 (MQ=255)
                                                                                                                                               accggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTc      >  2:831494/1‑150 (MQ=255)
                                                                                                                                                 cggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTCGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTCCa    <  1:2290817/150‑1 (MQ=255)
                                                                                                                                                  ggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTCCAg   <  2:1500970/150‑1 (MQ=255)
                                                                                                                                                   gcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTCCAGa  >  2:561128/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTCCAGA  >  NC_002947/1761154‑1761450

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑