Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,327,519 T→C 5.9% S163P (TCG→CCG)  PP_3797 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,327,5190TC 5.9% 201.2 / 2.6 68S163P (TCG→CCG) PP_3797membrane protein
Reads supporting (aligned to +/- strand):  ref base T (35/29);  new base C (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

ATTCAGGGCAGGACCCCAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTGGTTCATTGCGGTCATT  >  NC_002947/4327384‑4327662
                                                                                                                                         |                                                                                                                                                
aTTCAGGGCAGGACCCCAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCgg                                                                                                                                     >  4:213087/1‑149 (MQ=255)
 ttCAGGGCAGGACCCCAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGTTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGc                                                                                                                                    >  8:109252/1‑149 (MQ=255)
    aGGGCAGGACCCCAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCctcgctc                                                                                                                                             <  2:5849/137‑1 (MQ=255)
    aGGGCAGGACCCCAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGCCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCctcgctc                                                                                                                                             >  1:5849/1‑137 (MQ=255)
      ggCAGGACCCCAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCtggt                                                                                                                               <  8:207004/149‑1 (MQ=255)
          ggACCCCAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGtt                                                                                                                           <  7:303270/149‑1 (MQ=255)
           gACCCCAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGtt                                                                                                                           >  4:200995/1‑148 (MQ=255)
              cccAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTAc                                                                                                                        >  6:340208/1‑148 (MQ=255)
                cAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACGAAGTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTAc                                                                                                                        >  1:120423/1‑148 (MQ=255)
                cAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTAc                                                                                                                        >  5:104267/1‑146 (MQ=255)
                 aGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACaaga                                                                                                                    <  3:213087/149‑1 (MQ=255)
                              gAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTTTGGAAAg                                                                                                        <  5:120434/148‑1 (MQ=255)
                              gAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTTTGGAAAg                                                                                                        >  6:120434/1‑148 (MQ=255)
                              gAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAgg                                                                                                       <  1:232756/149‑1 (MQ=255)
                              gAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAgg                                                                                                       <  4:367490/149‑1 (MQ=255)
                              gAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAg                                                                                                        <  4:255018/148‑1 (MQ=255)
                                 gCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCtt                                                                                                    >  2:30518/1‑149 (MQ=255)
                                   cTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCttt‑t                                                                                                 >  7:33191/1‑149 (MQ=255)
                                    ttCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCttt‑tt                                                                                                >  1:83154/1‑149 (MQ=255)
                                    ttCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCttt‑t                                                                                                 >  4:131463/1‑148 (MQ=255)
                                         ccGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAAc                                                                                           <  5:340208/149‑1 (MQ=255)
                                         ccGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAAc                                                                                           <  2:120423/149‑1 (MQ=255)
                                         ccGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAAc                                                                                           >  5:113055/1‑149 (MQ=255)
                                           gcgcCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACT‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCTTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAAccc                                                                                         <  3:105738/149‑1 (MQ=255)
                                           gcgcCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAAccc                                                                                         >  3:286131/1‑149 (MQ=255)
                                             gcCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGt                                                                                       <  4:286131/149‑1 (MQ=255)
                                               ccTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTAGAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGc                                                                                     >  2:295562/1‑149 (MQ=255)
                                                   cACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACa                                                                                                                       <  7:74332/112‑1 (MQ=255)
                                                   cACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACa                                                                                                                       >  8:74332/1‑112 (MQ=255)
                                                     cGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCg                                                                               <  3:200995/149‑1 (MQ=255)
                                                         tggtggctggtgCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTc                                                                                               >  3:269496/1‑129 (MQ=255)
                                                         tggtggctggtgCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTc                                                                                               <  4:269496/129‑1 (MQ=255)
                                                             ggctggtgCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGt                                                                       >  2:92510/1‑149 (MQ=255)
                                                                tggtgCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTTCTTCAAACCCGTGCGTACCGGCCATGGTTc                                                                     <  1:30518/149‑1 (MQ=255)
                                                                  gtgCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCTGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTcgc                                                                   >  1:156874/1‑148 (MQ=255)
                                                                  gtgCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTcg                                                                    <  8:386204/147‑1 (MQ=255)
                                                                        ttCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAAc                                                                                           <  4:209422/118‑1 (MQ=255)
                                                                        ttCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAAc                                                                                           >  3:209422/1‑118 (MQ=255)
                                                                           aCCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTcgc                                                                   <  6:225390/139‑1 (MQ=255)
                                                                           aCCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTcgc                                                                   >  5:225390/1‑139 (MQ=255)
                                                                           aCCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTg                                                         >  8:246854/1‑149 (MQ=255)
                                                                           aCCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTg                                                         >  8:90631/1‑149 (MQ=255)
                                                                           aCCAACG‑‑CTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCCTACCGCCCATGCTTCGCGAATTTTCTg                                                         >  2:26973/1‑149 (MQ=255)
                                                                             cAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTggg                                                       <  7:90631/149‑1 (MQ=255)
                                                                                    gTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTAc                                                                                                                        <  4:131014/80‑1 (MQ=255)
                                                                                    gTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTAc                                                                                                                        >  3:131014/1‑80 (MQ=255)
                                                                                          gCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCAc                                            <  4:176726/149‑1 (MQ=255)
                                                                                           cTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGCCCTGCAcc                                           >  8:110952/1‑149 (MQ=255)
                                                                                           cTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGAccc                                                >  1:247265/1‑143 (MQ=255)
                                                                                           cTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGAccc                                                <  2:247265/144‑2 (MQ=255)
                                                                                               ggCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCt                                       <  1:92510/149‑1 (MQ=255)
                                                                                                 cTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCCCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCt                                       >  5:306241/1‑147 (MQ=255)
                                                                                                 cTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCCCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCt                                       <  6:311757/147‑1 (MQ=255)
                                                                                                 cTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCCCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCt                                       <  6:306241/147‑1 (MQ=255)
                                                                                                 cTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCCCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCt                                       >  5:311757/1‑147 (MQ=255)
                                                                                                 atGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGATTT‑ATCAAACCAGTGGGTACACGCAAAGCTTCGAGAATTTTAGGGGGTGACCGGAACCGTCTg                                      >  8:184770/2‑148 (MQ=255)
                                                                                                    gTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTcgc                                                                   <  3:27663/116‑1 (MQ=255)
                                                                                                    gTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTcgc                                                                   >  4:27663/1‑116 (MQ=255)
                                                                                                       ccTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCggg                                <  2:83154/148‑1 (MQ=255)
                                                                                                       ccTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGGTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCgggg                               >  5:201435/1‑149 (MQ=255)
                                                                                                                    ctCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTg                  >  7:132291/1‑149 (MQ=255)
                                                                                                                    ctCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTg                  >  6:282030/1‑149 (MQ=255)
                                                                                                                       ggCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTTCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTgg                 <  3:120160/147‑1 (MQ=255)
                                                                                                                       ggCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTGGtt               <  8:33191/149‑1 (MQ=255)
                                                                                                                       ggCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTGGt                <  7:246854/148‑1 (MQ=255)
                                                                                                                        gCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTGGtt               >  3:111197/1‑148 (MQ=255)
                                                                                                                                aTGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTGGTTCATTGCGGt      <  1:214530/149‑1 (MQ=255)
                                                                                                                                    tGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTGGTTCATTGCGGTCAtt  >  8:164042/1‑149 (MQ=255)
                                                                                                                                    tGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTGGTTCATTGCGGTCAtt  >  6:49137/1‑149 (MQ=255)
                                                                                                                                    tGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTGGTTCATTGCGGTCAtt  >  6:356817/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                
ATTCAGGGCAGGACCCCAGACTTCAACTTCGAAGCCTTCACCCGCGCCCTGCACGGCTGGTGGCTGGTGCCGTTCACCAACG‑‑GTTTCAGCTGGGGCTGGTACCTGGTATCGCTGCTCGGCCTGCCGATGCTGGCCTCGCTCGTCACCGGCCTGGTGGTTTACAAGAAGTTCTGGAAAGGCTTT‑TTCAAACCCGTGCGTACCGGCCATGGTTCGCGAATTTTCTGGGGTGACCTGCACCGTCTGTCCGGGGTGTGGTCAATCTGGTTCATTGCGGTCATT  >  NC_002947/4327384‑4327662

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