Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,798,289 A→C 6.7% L525R (CTG→CGG)  PP_5077 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,798,2890AC 6.7% 278.0 / 8.4 90L525R (CTG→CGG) PP_5077hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/40);  new base C (2/4);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

TGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTT  >  NC_002947/5798143‑5798437
                                                                                                                                                  |                                                                                                                                                     
tGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGc                                                                                                                                                     >  3:81282/1‑149 (MQ=255)
      cAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCTGGCCACGGCCCGctgata                                                                                                                                                <  3:61018/148‑5 (MQ=255)
      cAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttgt                                                                                                                                               >  6:22678/1‑149 (MQ=255)
       aaCTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttgtt                                                                                                                                              <  3:288414/149‑1 (MQ=255)
       aaCTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttgtt                                                                                                                                              <  4:81282/149‑1 (MQ=255)
              ggCTTGTGGACGGCAGTCCGGGACAGAAAAGCCTGCCGCTCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTTTTGGACGCt                                                                                                                                       <  5:279071/149‑1 (MQ=255)
                cTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGc                                                                                                                                     <  5:22678/149‑1 (MQ=255)
                    tgGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGAc                                                                                                                                 >  1:133590/1‑149 (MQ=255)
                      gCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACg                                                                                                                                <  8:105872/148‑1 (MQ=255)
                       ccGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGt                                                                                                                              <  2:133590/149‑1 (MQ=255)
                        cGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCTGGCCACGGCCCGctga                                                                                                                                                  <  8:63162/128‑3 (MQ=255)
                        cGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCTGGCCACGGCCCGctga                                                                                                                                                  >  5:104106/1‑126 (MQ=255)
                        cGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCTGGCCACGGCCCGctga                                                                                                                                                  <  6:104106/128‑3 (MQ=255)
                        cGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCTGGCCACGGCCCGctga                                                                                                                                                  >  7:63162/1‑126 (MQ=255)
                            aGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGATGGACGCTGCCGACGGTACGTg                                                                                                                         >  2:112242/1‑149 (MQ=255)
                                      aGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGt                                                                                                                 >  7:321836/1‑147 (MQ=255)
                                      aGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGt                                                                                                                 <  8:321836/147‑1 (MQ=255)
                                        aaaaGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAc                                                                                                             >  2:202414/1‑149 (MQ=255)
                                        aaaaGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAc                                                                                                             >  6:14150/1‑149 (MQ=255)
                                              cTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCt                                                                                                       >  7:320042/1‑149 (MQ=255)
                                              cTGCCGCGCTCGCGCGCATCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCt                                                                                                       <  6:224120/149‑1 (MQ=255)
                                                gCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTgg                                                                                                     >  3:13273/1‑149 (MQ=255)
                                                 ccGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCGAGAATTCGCTTAGCTACTGGGCGGCGGGGTGCTGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGa                                                                                                    >  1:42801/1‑149 (MQ=255)
                                                  cgcgCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGAc                                                                                                   <  2:42603/149‑1 (MQ=255)
                                                   gcgcTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGAcc                                                                                                  >  8:268886/1‑149 (MQ=255)
                                                   gcgcTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGAcc                                                                                                  >  7:245320/1‑149 (MQ=255)
                                                       tCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCg                                                                                              >  6:169264/1‑149 (MQ=255)
                                                         gggCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCg                                                                                              <  2:309403/147‑1 (MQ=255)
                                                          ggCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTg                                                                                           >  4:142237/1‑149 (MQ=255)
                                                           gcgcAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTa                                                                                                                             >  2:135166/1‑114 (MQ=255)
                                                           gcgcAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTa                                                                                                                             <  1:135166/114‑1 (MQ=255)
                                                               agcagGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAg                                                                                      <  6:21046/149‑1 (MQ=255)
                                                               agcagGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAg                                                                                      <  1:202414/149‑1 (MQ=255)
                                                                        gTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGa                                                                               <  8:245320/147‑1 (MQ=255)
                                                                         tCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTg                                                                            <  5:14150/149‑1 (MQ=255)
                                                                           gAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAc                                                                                                             >  6:128393/1‑114 (MQ=255)
                                                                           gAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAc                                                                                                             <  5:128393/114‑1 (MQ=255)
                                                                           gAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGcc                                                                          <  2:25482/149‑1 (MQ=255)
                                                                           gAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGcc                                                                          >  8:117943/1‑149 (MQ=255)
                                                                             aaGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGcc                                                                          <  4:228524/147‑1 (MQ=255)
                                                                                 tAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAAcc                                                                    >  4:315651/1‑149 (MQ=255)
                                                                                 tAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAAcc                                                                    >  2:238861/1‑149 (MQ=255)
                                                                                   gactccGCTTAGCTACTGGGCGGCGGGGGGATGCCGCCCGGAAGGTCTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGGGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGc                                                                  <  2:306161/144‑1 (MQ=255)
                                                                                   gAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCg                                                                   <  4:13273/148‑1 (MQ=255)
                                                                                      ttCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGc                                                               <  2:42801/149‑1 (MQ=255)
                                                                                      ttCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGc                                                               <  5:169264/149‑1 (MQ=255)
                                                                                         gCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTc                                                            >  6:134488/1‑149 (MQ=255)
                                                                                           ttAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgc                                                          >  7:153849/1‑149 (MQ=255)
                                                                                               cTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGtt                                                      >  8:229390/1‑149 (MQ=255)
                                                                                               cTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGtt                                                      >  2:315754/1‑149 (MQ=255)
                                                                                               cTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGtt                                                      >  5:48824/1‑149 (MQ=255)
                                                                                                tACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCttta                                                                               <  5:99948/123‑3 (MQ=255)
                                                                                                tACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGa                                                                               >  6:99948/1‑123 (MQ=255)
                                                                                                  cTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTg                                                     <  6:48824/147‑1 (MQ=255)
                                                                                                  cTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGc                                                   >  4:192484/1‑149 (MQ=255)
                                                                                                   tGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGc                                                                  >  2:11965/1‑133 (MQ=255)
                                                                                                   tGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGc                                                                  <  1:11965/133‑1 (MQ=255)
                                                                                                     ggcggcGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTg                                                     >  8:132558/1‑144 (MQ=255)
                                                                                                     ggcggcGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTg                                                     <  7:132558/144‑1 (MQ=255)
                                                                                                     ggcggcGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑a                                                >  2:334192/1‑148 (MQ=255)
                                                                                                       cggcggGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAAc                                             >  2:231873/1‑149 (MQ=255)
                                                                                                             ggTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGc                                                               >  4:168708/1‑126 (MQ=255)
                                                                                                             ggTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGc                                                               <  3:168708/126‑1 (MQ=255)
                                                                                                               tGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTaa                                     <  5:134488/149‑1 (MQ=255)
                                                                                                                    ccgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCac                                >  8:120520/1‑149 (MQ=255)
                                                                                                                    ccgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCac                                >  6:50179/1‑149 (MQ=255)
                                                                                                                    ccgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCa                                 <  7:229390/148‑1 (MQ=255)
                                                                                                                      gccCGGAATGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTTCCGTCGGTACGTGGCCACGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCacga                              <  1:112242/149‑1 (MQ=255)
                                                                                                                      gccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTaa                                     >  5:170302/1‑142 (MQ=255)
                                                                                                                      gccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTaa                                     <  6:170302/142‑1 (MQ=255)
                                                                                                                      gccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCacga                              <  3:192484/149‑1 (MQ=255)
                                                                                                                           gAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTAc                         <  7:268886/149‑1 (MQ=255)
                                                                                                                             aGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTa                          >  7:295239/1‑146 (MQ=255)
                                                                                                                              ggTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACaga                      <  7:117943/149‑1 (MQ=255)
                                                                                                                              ggTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACaga                      <  1:231873/149‑1 (MQ=255)
                                                                                                                                 aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCg                                        >  5:113399/1‑128 (MQ=255)
                                                                                                                                 aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCg                                        <  6:113399/128‑1 (MQ=255)
                                                                                                                                 aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGc                   >  8:108803/1‑149 (MQ=255)
                                                                                                                                       gggCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTg                >  2:66215/1‑146 (MQ=255)
                                                                                                                                       gggCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTg                <  1:66215/146‑1 (MQ=255)
                                                                                                                                        ggCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCt            >  2:326079/1‑149 (MQ=255)
                                                                                                                                         gCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTg           <  7:108803/149‑1 (MQ=255)
                                                                                                                                          ccGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTg           >  6:287931/1‑148 (MQ=255)
                                                                                                                                           cGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCa         <  1:238861/149‑1 (MQ=255)
                                                                                                                                            acggCCCGCTCTTGTTGTACGCTGCAGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGTCAATTTCTTGATTGCAGCAACCGCTGCGTCGCGTTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAg        <  2:202374/147‑1 (MQ=255)
                                                                                                                                               gCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGtt     >  3:166745/1‑149 (MQ=255)
                                                                                                                                                ccAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTa                          <  5:336507/127‑1 (MQ=255)
                                                                                                                                                ccAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTa                          >  6:336507/1‑127 (MQ=255)
                                                                                                                                                 cAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGtttt   >  2:174777/1‑149 (MQ=255)
                                                                                                                                                 cAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGtttt   >  4:287870/1‑149 (MQ=255)
                                                                                                                                                  aGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGttt    <  1:71638/148‑1 (MQ=255)
                                                                                                                                                  aGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGtttta  >  2:71638/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCG‑AAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTT  >  NC_002947/5798143‑5798437

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