Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,880,092 C→G 5.3% I234M (ATC→ATG PP_4289 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,880,0920CG 5.3% 245.5 / 3.2 75I234M (ATC→ATGPP_4289membrane protein
Reads supporting (aligned to +/- strand):  ref base C (33/38);  new base G (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TCGGTAACGTGTTCCGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGCAACAAGTACGCCTGGACCCTGCCCGTCGGC  >  NC_002947/4879944‑4880239
                                                                                                                                                    |                                                                                                                                                   
tCGGTAACGTGTTCCGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTtcatc                                                                                                                                                     <  4:181842/149‑1 (MQ=255)
tCGGTAACGTGTTCCGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTtcatc                                                                                                                                                     <  5:311357/149‑1 (MQ=255)
      aCGTGTTCCGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATc                                                                                                                                               <  5:307667/149‑1 (MQ=255)
        gtgtTCCGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAg                                                                                                                                             >  2:79894/1‑149 (MQ=255)
        gtgtTCCGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGATCATGATGAACAg                                                                                                                                             >  1:146400/1‑149 (MQ=255)
        gtgtTCCGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCGGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAg                                                                                                                                             >  2:74986/1‑149 (MQ=255)
          gtTCCGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTACCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGc                                                                                                                                            >  7:42966/1‑148 (MQ=255)
            tCCGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGATCATGATGAACAGCac                                                                                                                                          <  2:146400/148‑2 (MQ=255)
             ccGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAAc                                                                                                                                         >  6:254089/1‑148 (MQ=255)
                catcatcatGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCAtt                                                                                                                                     <  1:296397/149‑1 (MQ=255)
                    atcatGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTccc                                                                                                                                 <  6:48607/149‑1 (MQ=255)
                     tcatGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTccc                                                                                                                                 >  3:34890/1‑148 (MQ=255)
                       atGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGa                                                                                                                               <  4:190748/148‑1 (MQ=255)
                       atGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGATCATGATGAACAGCac                                                                                                                                          <  5:272374/137‑2 (MQ=255)
                       atGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGATCATGATGAACAGCac                                                                                                                                          >  6:272374/1‑136 (MQ=255)
                       atGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTACCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGa                                                                                                                               <  8:42966/148‑1 (MQ=255)
                        tGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTcc                                                                                                                                  >  4:330011/1‑144 (MQ=255)
                        tGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTcc                                                                                                                                  <  3:330011/144‑1 (MQ=255)
                            ggCCCATCGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCTGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCt                                                                                                                         <  5:175743/149‑1 (MQ=255)
                            ggCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCAcc                                                                                                                          >  8:320472/1‑148 (MQ=255)
                               ccAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCt                                                                                                                         <  4:289416/146‑1 (MQ=255)
                               ccAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCt                                                                                                                         >  3:289416/1‑146 (MQ=255)
                               ccAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGATGTTCATCATGATCAGCAACCAGTTCCCGAGCACCTACg                                                                                                                      >  2:240030/1‑149 (MQ=255)
                                   cGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTcc                                                                                                                                  >  5:352807/1‑133 (MQ=255)
                                   cGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTcc                                                                                                                                  <  6:352807/133‑1 (MQ=255)
                                        gCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAg                                                                                                              >  8:73595/1‑148 (MQ=255)
                                            gtggCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAAc                                                                                                                                         <  6:49067/117‑1 (MQ=255)
                                                 ggcgATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGc                                                                                                    >  5:269744/1‑149 (MQ=255)
                                                         agaACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCTCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTTGTACTGGCtaccgtag                                                                                            <  1:74986/149‑8 (MQ=255)
                                                         agaACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTTGTACTGGCtaccgtag                                                                                            <  1:79894/149‑8 (MQ=255)
                                                         agaACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTg                                                                                            <  8:239876/149‑1 (MQ=255)
                                                          gaACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTgg                                                                                           <  4:318460/149‑1 (MQ=255)
                                                            acaacCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTTCTGCGCTCTCGTCACAACAACTACTTCACCCTGCCTGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGc                                                                                          <  1:240030/148‑1 (MQ=255)
                                                                ccAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTa                                                                                     >  2:227052/1‑149 (MQ=255)
                                                                  aGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTg                                                                                            >  3:271698/1‑140 (MQ=255)
                                                                  aGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTg                                                                                            <  4:271698/140‑1 (MQ=255)
                                                                  aGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCTGTATc                                                                                   <  7:320472/149‑1 (MQ=255)
                                                                        ccGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAATCATTTCCCGAGCACCTACGGTAGCCAGTTGTACTGGCta                                                                                                  >  1:165481/1‑127 (MQ=255)
                                                                        ccGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAATCATTTCCCGAGCACCTACGGTAGCCAGTTGTACTGGCta                                                                                                  <  2:165481/128‑2 (MQ=255)
                                                                        ccGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAg                                                                               >  2:300712/1‑147 (MQ=255)
                                                                         cGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAg                                                                            >  7:187910/1‑149 (MQ=255)
                                                                             ccGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCcgc                                                                        <  7:73595/149‑1 (MQ=255)
                                                                                   cTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTc                                                                   <  8:187910/148‑1 (MQ=255)
                                                                                        ggCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATcc                                                             <  1:300712/149‑1 (MQ=255)
                                                                                           cAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCg                                                                                        <  1:289159/119‑1 (MQ=255)
                                                                                           cAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCg                                                                                        >  2:289159/1‑119 (MQ=255)
                                                                                           cAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGc                                                           >  2:209015/1‑148 (MQ=255)
                                                                                           cAAGGGCCTGATGAGCTCGAGTCACAACAAATACTTCACCCTGCCGGTGCTCTTCATCATGATCAGCAAGCAGTTCCCGAGCACCTACGGTAGCCAGTACAACTGACTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATACGcc                                                          >  2:239411/1‑149 (MQ=255)
                                                                                             aGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCac                                                        <  1:209015/149‑1 (MQ=255)
                                                                                              gggCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCa                                                                                <  7:82312/124‑1 (MQ=255)
                                                                                              gggCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCa                                                                                >  8:82312/1‑124 (MQ=255)
                                                                                              gggCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCa                                                                                <  1:145626/124‑1 (MQ=255)
                                                                                              gggCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCa                                                                                >  2:145626/1‑124 (MQ=255)
                                                                                                     ctgcGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTc                                                                   <  4:252320/130‑1 (MQ=255)
                                                                                                     ctgcGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTc                                                                   >  3:252320/1‑130 (MQ=255)
                                                                                                     ctgcGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGATAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTc                                                  <  4:34890/147‑1 (MQ=255)
                                                                                                                   caacaaCTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGa                                    >  5:18985/1‑147 (MQ=255)
                                                                                                                    aacaacTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTc                                                                   <  8:40240/115‑1 (MQ=255)
                                                                                                                    aacaacTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTc                                                                   >  7:40240/1‑115 (MQ=255)
                                                                                                                    aacaacTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACAccc                                           >  1:16859/1‑139 (MQ=255)
                                                                                                                    aacaacTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACAccc                                           <  2:16859/139‑1 (MQ=255)
                                                                                                                     acaacTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGc                                <  1:227052/149‑1 (MQ=255)
                                                                                                                        actactTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACtt                                                   >  1:128980/1‑127 (MQ=255)
                                                                                                                        actactTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACtt                                                   <  2:128980/127‑1 (MQ=255)
                                                                                                                         ctactTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGcaaca                            <  5:254089/149‑1 (MQ=255)
                                                                                                                            ctTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCg                                          <  7:154937/132‑1 (MQ=255)
                                                                                                                            ctTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCg                                          >  8:154937/1‑132 (MQ=255)
                                                                                                                            ctTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGCAACAAGt                         <  6:18985/149‑1 (MQ=255)
                                                                                                                              tCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGcaacaa                           <  7:164220/145‑1 (MQ=255)
                                                                                                                              tCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGcaacaa                           >  8:164220/1‑145 (MQ=255)
                                                                                                                                  ccTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGCAACAAGTACGCCt                   <  6:269744/149‑1 (MQ=255)
                                                                                                                                  ccTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGCAACAAGTACGCCt                   >  5:295794/1‑149 (MQ=255)
                                                                                                                                   cTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGCAACAAGTACGCCTg                  >  1:146950/1‑149 (MQ=255)
                                                                                                                                      ccGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGcaacaa                           <  2:58142/137‑1 (MQ=255)
                                                                                                                                      ccGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGcaacaa                           >  1:58142/1‑137 (MQ=255)
                                                                                                                                      caggTGCTGTTAATCATGATCAGCAACCATTTCCCGAGCACCTAAGGTAGCCAGTACAACTGGCTGATCCTGCCCGGTATCGCAGTAGCCGCTGTTCTGATCCGCAACTACTTCAACACCCGCCACGACAGCAACAAGTACGCCttcac               >  4:109767/3‑145 (MQ=255)
                                                                                                                                            cTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGCAACAAGTACGCCTGGACCCTGc          <  6:295794/148‑1 (MQ=255)
                                                                                                                                                  atcatGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCa                                                         >  7:273633/1‑95 (MQ=255)
                                                                                                                                                  atcatGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCa                                                         <  8:273633/95‑1 (MQ=255)
                                                                                                                                                   tcatGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGCAACAAGTACGCCTGGACCCTGCCCGTCGgc  >  1:114326/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGGTAACGTGTTCCGCATCATCATGCCGGCCCAACGCCAGCTGGTGGCGGCGATCGAGAACAACCAGACGCCCGACCCGGTCCTGCCGGCCAAGGGCCTGCTGCGCTCGCGTCACAACAACTACTTCACCCTGCCGGTGCTGTTCATCATGATCAGCAACCATTTCCCGAGCACCTACGGTAGCCAGTACAACTGGCTGATCCTGGCCGGTATCGCAGTAGCCGCGGTTCTGATCCGCCACTACTTCAACACCCGCCACGACAGCAACAAGTACGCCTGGACCCTGCCCGTCGGC  >  NC_002947/4879944‑4880239

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: