Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 181,275 A→T 5.4% H131L (CAC→CTC)  degP → serine endoprotease (protease Do), membrane‑associated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913181,2750AT 5.4% 357.8 / 7.1 111H131L (CAC→CTC) degPserine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base A (51/54);  new base T (2/4);  total (53/58)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 113 total aligned reads displayed.

CTGCCAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGCGGAT  >  NC_000913/181129‑181423
                                                                                                                                                  |                                                                                                                                                    
cTGCCAGGAAGGTTCTACGTTCCAGAGCTCGCCGTTCTGCCAGTGTGGCCAGGGCGGTAATGTTGGCGTCCATCAACAGAAATTCATGGCTCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCAc                                                                                                                                                     <  8:328602/148‑1 (MQ=255)
   ccAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttg                                                                                                                                                 >  8:82987/1‑149 (MQ=255)
    cAGGAAGGTTGTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgt                                                                                                                                                <  3:148219/149‑1 (MQ=255)
    cAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGTGGCCAGCTACAGAATTTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgt                                                                                                                                                <  2:3086/149‑1 (MQ=255)
    cAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgt                                                                                                                                                >  2:187724/1‑149 (MQ=255)
     aGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgtt                                                                                                                                               >  4:30391/1‑149 (MQ=255)
     aGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgtt                                                                                                                                               >  2:143569/1‑149 (MQ=255)
       gAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGa                                                                                                                                             <  6:5135/149‑1 (MQ=255)
       gAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGa                                                                                                                                             <  6:128844/149‑1 (MQ=255)
        aaGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGttgttggt                                                                                                                                            <  5:332669/149‑3 (MQ=255)
        aaGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGAt                                                                                                                                            <  2:62105/149‑1 (MQ=255)
            ttCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAAcg                                                                                                                                        <  4:346560/149‑1 (MQ=255)
              ctcCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAAcgc                                                                                                                                       <  1:340883/148‑1 (MQ=255)
                      cAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCAt                                                                                                                              <  1:143569/149‑1 (MQ=255)
                      cAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCAt                                                                                                                              <  8:123320/149‑1 (MQ=255)
                      cAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGGTCCGGCGTCATCATTGATGCCGATAACGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCAt                                                                                                                              <  8:129817/149‑1 (MQ=255)
                       agagCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCAt                                                                                                                              >  4:431652/1‑148 (MQ=255)
                       agagCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCAt                                                                                                                              >  1:252787/1‑148 (MQ=255)
                        gagCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCAtt                                                                                                                             >  1:167127/1‑148 (MQ=255)
                          gCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTaaa                                                                                                                          <  7:178846/149‑1 (MQ=255)
                              tcCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGtt                                                                                                                       >  6:261017/1‑148 (MQ=255)
                              tcCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGtt                                                                                                                       <  4:246916/148‑1 (MQ=255)
                              tcCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGtt                                                                                                                       >  3:246916/1‑148 (MQ=255)
                                       ccAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGa                                                                                                                <  7:420490/146‑1 (MQ=255)
                                        cAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGa                                                                                                            >  7:332724/1‑149 (MQ=255)
                                          gggTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATg                                                                                                          >  7:247384/1‑149 (MQ=255)
                                           ggTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATgg                                                                                                         <  5:261017/149‑1 (MQ=255)
                                             tGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGcc                                                                                                       <  7:82987/149‑1 (MQ=255)
                                                  aGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgtt                                                                                                                                                  >  7:149413/1‑101 (MQ=255)
                                                  aGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgtt                                                                                                                                                  <  8:149413/101‑1 (MQ=255)
                                                    ggCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCTTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGt                                                                                                 >  2:279023/1‑148 (MQ=255)
                                                       ggTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGa                                                                                                                <  1:344000/130‑1 (MQ=255)
                                                       ggTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGa                                                                                                                >  2:344000/1‑130 (MQ=255)
                                                         tAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGAcg                                                                                           >  7:323186/1‑149 (MQ=255)
                                                         tAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGAcg                                                                                           >  7:253069/1‑149 (MQ=255)
                                                           aTGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGttgttg                                                                                                                                              <  5:282148/96‑1 (MQ=255)
                                                           aTGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGttgttg                                                                                                                                              >  6:282148/1‑96 (MQ=255)
                                                                ggcggcCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAACGCCTATGTCGTCACCAACACCCACGATGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGATCGACGCGAAGAt                                                                                    >  2:212493/1‑149 (MQ=255)
                                                                   ggcCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGt                                                                                 >  7:86796/1‑149 (MQ=255)
                                                                   ggcCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGt                                                                                 <  3:431652/149‑1 (MQ=255)
                                                                    gcCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGtt                                                                                <  7:19641/149‑1 (MQ=255)
                                                                     ccAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTg                                                                               <  2:252787/149‑1 (MQ=255)
                                                                      cAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGttgttggtga                                                                                                                                          >  1:381544/1‑85 (MQ=255)
                                                                      cAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGGCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGttgttggtga                                                                                                                                          <  2:381544/89‑5 (MQ=255)
                                                                        gCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgt                                                                                                                                                <  3:365016/81‑1 (MQ=255)
                                                                         cAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGa                                                                                     >  5:95755/1‑139 (MQ=255)
                                                                         cAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGa                                                                                     <  6:95755/139‑1 (MQ=255)
                                                                          aaCAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCa                                                                            >  1:353040/1‑147 (MQ=255)
                                                                           aCAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCa                                                                            >  3:49941/1‑146 (MQ=255)
                                                                            cAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCa                                                                            >  7:75667/1‑145 (MQ=255)
                                                                             aGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGAt                                                                       <  3:30391/149‑1 (MQ=255)
                                                                             aGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGAt                                                                       <  8:332724/149‑1 (MQ=255)
                                                                               aaaTTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATcc                                                                     >  3:35071/1‑149 (MQ=255)
                                                                                 aTTCATGGCGCTGGGTTCCTGCGTCATCATTGATGCCGATAACCGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAAATGCGAGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCcgc                                                                   >  1:2086/1‑149 (MQ=255)
                                                                                  ttCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGtt                                                                                                >  4:104905/1‑119 (MQ=255)
                                                                                  ttCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGtt                                                                                                <  3:104905/119‑1 (MQ=255)
                                                                                     aTGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGctc                                                               <  1:187724/149‑1 (MQ=255)
                                                                                          gcTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGata                                                          >  2:49950/1‑149 (MQ=255)
                                                                                           cTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGatat                                                         <  8:75667/149‑1 (MQ=255)
                                                                                            tGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATc                                                        >  1:391974/1‑149 (MQ=255)
                                                                                                ttCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCgcgc                                                    <  1:49950/149‑1 (MQ=255)
                                                                                                ttCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCgcgc                                                    >  8:200951/1‑149 (MQ=255)
                                                                                                    ggCGTCATCATTGATGCCGATAAAGGCTATATCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAACGATCCGCGCTCTGATATCGCGCTGaa                                                >  3:350348/1‑148 (MQ=255)
                                                                                                     gCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATc                                               >  1:287135/1‑149 (MQ=255)
                                                                                                     gCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATc                                               >  3:183773/1‑149 (MQ=255)
                                                                                                      cGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATcc                                              <  2:167127/149‑1 (MQ=255)
                                                                                                           tcatTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATc                                         <  2:353040/149‑1 (MQ=255)
                                                                                                            catTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCaaa                                           >  6:1467/1‑146 (MQ=255)
                                                                                                            catTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCaaa                                           <  5:1467/146‑1 (MQ=255)
                                                                                                               tGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAg                                      <  8:253069/148‑1 (MQ=255)
                                                                                                                 aTGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGaa                                    <  6:42977/148‑1 (MQ=255)
                                                                                                                  tGCCGATAAAGGCTATGTCGTCACCAACAACCTCGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAAcc                                  <  3:101353/149‑1 (MQ=255)
                                                                                                                  tGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAAcc                                  <  2:391974/149‑1 (MQ=255)
                                                                                                                    ccGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAAccc                                 <  8:323186/148‑1 (MQ=255)
                                                                                                                     cGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGa                                                                                     >  5:68416/1‑95 (MQ=255)
                                                                                                                     cGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGa                                                                                     <  6:68416/95‑1 (MQ=255)
                                                                                                                     cGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCg                                >  4:21662/1‑148 (MQ=255)
                                                                                                                     cGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGa                               <  3:21662/149‑1 (MQ=255)
                                                                                                                           aGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGaaaaaac                         >  2:364737/1‑147 (MQ=255)
                                                                                                                            ggCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCt                        <  4:49941/149‑1 (MQ=255)
                                                                                                                              cTATGGCGGCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGa                      <  6:427536/149‑1 (MQ=255)
                                                                                                                                   tcgtcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGc                  >  2:285503/1‑148 (MQ=255)
                                                                                                                                   tcgtcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGc                  >  5:315738/1‑148 (MQ=255)
                                                                                                                                    cgtcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGCTCCAAATCCAGACACCGAAAAACCTCACCGCaa                >  2:465151/1‑149 (MQ=255)
                                                                                                                                       cACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCg                                >  3:314574/1‑130 (MQ=255)
                                                                                                                                       cACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCg                                <  4:314574/130‑1 (MQ=255)
                                                                                                                                        aCCAACATCCACGTTGTTGATTACGCGACGGTCATTATAGTTCTACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTGTGATATCGCGCTGATCCAAATCCAGAACCCGATAATCCTGACCGCAATTaa            <  2:2086/149‑1 (MQ=255)
                                                                                                                                        aCCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATcc                                        >  8:27168/1‑121 (MQ=255)
                                                                                                                                        aCCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATcc                                        <  7:27168/121‑1 (MQ=255)
                                                                                                                                         ccAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAAtt              >  5:163659/1‑146 (MQ=255)
                                                                                                                                         ccAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAAtt              >  6:135812/1‑146 (MQ=255)
                                                                                                                                             caacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGaa                                    <  4:231934/120‑1 (MQ=255)
                                                                                                                                              aacCACGTTGTTGATAACGCGACTGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGc      <  1:279023/149‑1 (MQ=255)
                                                                                                                                              aacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGc      <  6:315738/149‑1 (MQ=255)
                                                                                                                                              aacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGc      >  1:132389/1‑149 (MQ=255)
                                                                                                                                               acCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAAc                          >  1:205300/1‑128 (MQ=255)
                                                                                                                                               acCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAAc                          <  2:205300/128‑1 (MQ=255)
                                                                                                                                                 cACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGCGGa   <  7:200951/149‑1 (MQ=255)
                                                                                                                                                 cACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGCGGa   >  2:290033/1‑149 (MQ=255)
                                                                                                                                                  aCGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGCGGAt  >  7:105863/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTGCCAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGCGGAT  >  NC_000913/181129‑181423

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: