Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 834,159 T→C 5.2% W279R (TGG→CGG)  PP_0717 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947834,1590TC 5.2% 229.7 / 2.2 76W279R (TGG→CGG) PP_0717membrane protein
Reads supporting (aligned to +/- strand):  ref base T (31/41);  new base C (2/2);  total (33/43)
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.55e-01

CTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTATGCCAGCGCCATGTTGATCTGGAT  >  NC_002947/834012‑834307
                                                                                                                                                   |                                                                                                                                                     
cttgcCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTg                                                                                                                                                      <  8:114101/149‑1 (MQ=255)
   gcCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTgg                                                                                                                                                     <  2:351521/147‑1 (MQ=255)
    ccGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGt                                                                                                                                                   <  3:131855/148‑1 (MQ=255)
    ccGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGt                                                                                                                                                   <  5:3673/148‑1 (MQ=255)
     cGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGt                                                                                                                                                   >  7:9275/1‑147 (MQ=255)
        ggTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGtt                                                                                                                                              <  8:330745/149‑1 (MQ=255)
           gtTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTgg                                                                                                                                                     <  4:80068/139‑1 (MQ=255)
           gtTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTgg                                                                                                                                                     >  3:80068/1‑139 (MQ=255)
           gtTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGcc                                                                                                                                           <  3:291399/149‑1 (MQ=255)
            ttGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCt                                                                                                                                         >  6:58525/1‑149 (MQ=255)
            ttGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCt                                                                                                                                         >  6:52732/1‑149 (MQ=255)
            ttGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCt                                                                                                                                         >  8:208272/1‑149 (MQ=255)
            ttGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCt                                                                                                                                         <  6:298100/149‑1 (MQ=255)
              ggCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTga                                                                                                                                       >  5:49319/1‑149 (MQ=255)
                 gATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTgacga                                                                                                                                    >  4:364035/1‑149 (MQ=255)
                       gcTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTAGCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGt                                                                                                                              <  1:3288/149‑1 (MQ=255)
                       gcTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGt                                                                                                                              <  7:146228/149‑1 (MQ=255)
                        cTGTGGGTATCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGtt                                                                                                                             <  3:205759/149‑1 (MQ=255)
                        cTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGtt                                                                                                                             <  3:76104/149‑1 (MQ=255)
                        cTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGtt                                                                                                                             <  3:210639/149‑1 (MQ=255)
                             ggTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTGTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATc                                                                                                                         <  8:90292/148‑1 (MQ=255)
                               tAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTg                                                                                                                      >  5:246556/1‑149 (MQ=255)
                                 gcgcTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGc                                                                                                                     <  5:196286/148‑1 (MQ=255)
                                 gcgcTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCt                                                                                                                    <  6:49319/149‑1 (MQ=255)
                                   gcTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTgg                                                                                                                  <  3:364035/149‑1 (MQ=255)
                                    cTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTggg                                                                                                                 <  6:246556/149‑1 (MQ=255)
                                       aGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTgggg                                                                                                                <  3:368451/147‑1 (MQ=255)
                                             ttcgttcTCCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCg                                                                                                          <  5:58525/147‑1 (MQ=255)
                                             ttcgttcACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCg                                                                                                          <  5:52732/147‑1 (MQ=255)
                                                 ttcACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAg                                                                                                    <  5:104543/149‑1 (MQ=255)
                                                     ccTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGGAGGCTCGGCGCCTGTACGAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACa                                                                                                <  8:167915/149‑1 (MQ=255)
                                                     ccTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCactgtc                                                                                                <  7:208272/149‑6 (MQ=255)
                                                        gtgGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTg                                                                                             >  4:66500/1‑149 (MQ=255)
                                                           gATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGt                                                                                                                                                   <  1:125119/93‑1 (MQ=255)
                                                           gATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGt                                                                                                                                                   >  2:125119/1‑93 (MQ=255)
                                                              gTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCATGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTg                                                                                       >  6:195993/1‑149 (MQ=255)
                                                                     ttGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGa                                                                                >  5:186118/1‑149 (MQ=255)
                                                                              ctgctgTTCGCCCTTTGGGAGGCTCGGCGGCTGTACCAGGCGGGGCCGAGCCGGCAGATCCCGGCGTACTGGGTCAGTTGCCTGTCGACCATGGTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACaa                                                                       <  4:56260/149‑1 (MQ=255)
                                                                               tgctgTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAg                                                                      >  4:246215/1‑149 (MQ=255)
                                                                                 ctgTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACCGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGAtgt                                                                             >  6:20785/1‑140 (MQ=255)
                                                                                 ctgTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACCGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGAtgt                                                                             <  5:20785/140‑1 (MQ=255)
                                                                                 ctgTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACCGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGAtgt                                                                             >  6:26094/1‑140 (MQ=255)
                                                                                 ctgTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACCGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGAtgt                                                                             <  5:26094/140‑1 (MQ=255)
                                                                                     acGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGTGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGAtgt                                                                             <  4:58971/135‑1 (MQ=255)
                                                                                     acGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGTGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGAtgt                                                                             >  3:58971/2‑136 (MQ=255)
                                                                                        cccTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCg                                                             >  7:126331/1‑149 (MQ=255)
                                                                                        cccTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCg                                                             <  8:126331/149‑1 (MQ=255)
                                                                                         ccTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGAtgt                                                                             >  3:246427/1‑132 (MQ=255)
                                                                                         ccTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGAtgt                                                                             <  4:246427/132‑1 (MQ=255)
                                                                                          cTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGt                                                           <  8:9275/149‑1 (MQ=255)
                                                                                               gCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGt                                                       <  3:246215/148‑1 (MQ=255)
                                                                                                 aGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATg                                                    >  6:114784/1‑149 (MQ=255)
                                                                                                     tCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTg                                                                                             <  8:86533/104‑1 (MQ=255)
                                                                                                     tCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGTAGGATGTGTACAAGGCATTTGc                                                              >  6:213524/1‑135 (MQ=255)
                                                                                                     tCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGTAGGATGTGTACAAGGCATTTGc                                                              <  5:213524/135‑1 (MQ=255)
                                                                                                     tCGGCGCCTCTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTg                                                                                             >  7:86533/1‑104 (MQ=255)
                                                                                                        gcgcCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGt                                             <  1:45637/149‑1 (MQ=255)
                                                                                                               tACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCGTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTgg                                       >  1:86680/1‑149 (MQ=255)
                                                                                                               tACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGa                                                                                                  <  7:123962/89‑1 (MQ=255)
                                                                                                               tACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGa                                                                                                  >  8:123962/1‑89 (MQ=255)
                                                                                                               tACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACaa                                                                       <  2:38686/116‑1 (MQ=255)
                                                                                                               tACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACaa                                                                       >  1:38686/1‑116 (MQ=255)
                                                                                                               tACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTgg                                       <  1:157219/148‑1 (MQ=255)
                                                                                                                  cAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGcggt                                   >  4:89602/1‑149 (MQ=255)
                                                                                                                       ggTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGt                                                                                        <  4:34230/91‑1 (MQ=255)
                                                                                                                       ggTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGt                                                                                        >  3:34230/1‑91 (MQ=255)
                                                                                                                           ccGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTa                          <  8:39173/149‑1 (MQ=255)
                                                                                                                             gAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTgggg                                                                                                                <  6:108751/61‑1 (MQ=255)
                                                                                                                             gAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTgggg                                                                                                                >  5:108751/1‑61 (MQ=255)
                                                                                                                             gAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTTTATg                        >  6:261025/1‑149 (MQ=255)
                                                                                                                             gAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTATg                        >  5:251292/1‑149 (MQ=255)
                                                                                                                             gAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTATg                        >  5:245816/1‑149 (MQ=255)
                                                                                                                               gCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTCGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTATGcc                      >  3:298859/1‑149 (MQ=255)
                                                                                                                                     agagCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTATGCCAgcgc                 >  8:20371/1‑148 (MQ=255)
                                                                                                                                     agagCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTATGCCAgcgc                 >  4:366950/1‑148 (MQ=255)
                                                                                                                                                  cTGGGTCAGTTGCCTGACGACCATGTTGATCTCGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTATGCCAGCGCCATGTTGATCTGGa   <  4:298859/149‑1 (MQ=255)
                                                                                                                                                   tGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTATGCCAGCGCCATGTTGATCTg     >  8:15026/1‑146 (MQ=255)
                                                                                                                                                   tGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTATGCCAGCGCCATGTTGATCTGGAt  >  7:171793/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTGTACCAGGCGGTGCCGAGCCGGCAGAGCCCGGCGTACTGGGTCAGTTGCCTGACGACCATGTTGATCTTGCTGGGGCAGTCGGTGCAGGACAGTGCAGGTGGGCAGGATGTGTACAAGGCATTTGCGGTGCGTATGGCGCTGTTTCTGGCGGTGTCGGTGTATGCCAGCGCCATGTTGATCTGGAT  >  NC_002947/834012‑834307

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