Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,001,039 G→A 5.7% A241T (GCC→ACC)  dmoA‑II → dimethyl‑sulfide monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,001,0390GA 5.7% 216.2 / 2.0 70A241T (GCC→ACC) dmoA‑IIdimethyl‑sulfide monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (34/32);  new base A (2/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GAACCCTCGCGGGTACACCCGATCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGTTGATCGTCGCCCCCA  >  NC_002947/4000892‑4001180
                                                                                                                                                   |                                                                                                                                             
gAACCCTCGCGGGTACACCCGATCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGc                                                                                                                                              <  3:20628/149‑1 (MQ=255)
 aaCCCTCGCGGGTACACCCGATCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGcc                                                                                                                                             >  4:243203/1‑149 (MQ=255)
 aaCCCTCGCGGGTACACCCGATCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGcc                                                                                                                                             >  5:243951/1‑149 (MQ=255)
    ccTCGCGGGTACACCCGATCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGaa                                                                                                                                          <  3:262407/149‑1 (MQ=255)
           ggTACACCCGATCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCGCCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCgtgt                                                                                                                                   <  7:70754/149‑1 (MQ=255)
           ggTACACCCGATCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATgtgtgt                                                                                                                                   <  2:27787/149‑1 (MQ=255)
                acCCGATCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTtcat                                                                                                                               >  2:149559/1‑148 (MQ=255)
                 cccGATCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTtcat                                                                                                                               >  1:132013/1‑147 (MQ=255)
                     aTCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCgg                                                                                                                         <  4:82450/149‑1 (MQ=255)
                         aCCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGcc                                                                                                                      >  1:214792/1‑148 (MQ=255)
                         aCCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGcc                                                                                                                      >  2:145566/1‑148 (MQ=255)
                           cATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAAc                                                                                                                   <  1:145566/149‑1 (MQ=255)
                           cATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAAc                                                                                                                   <  6:243951/149‑1 (MQ=255)
                                          ttCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGACCCGCGGCCCGCAGTTTGCCGCACGGCATGCCGACTGCGTGGTCATCAGCGGGCCAACGCCGACGGTGCTgcg                                                                                                    >  5:191669/1‑149 (MQ=255)
                                          ttCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGACACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTgcg                                                                                                    >  5:186958/1‑149 (MQ=255)
                                            ccAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTgcg                                                                                                    <  2:214792/147‑1 (MQ=255)
                                            ccAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGcgcc                                                                                                  >  7:136894/1‑149 (MQ=255)
                                             cAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGcgccg                                                                                                 >  5:182654/1‑149 (MQ=255)
                                                   cccGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGc                                                                                           >  7:6106/1‑149 (MQ=255)
                                                     cGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGcc                                                                                          >  4:202054/1‑148 (MQ=255)
                                                       gCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGAc                                                                                       <  1:299889/149‑1 (MQ=255)
                                                       gCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGAc                                                                                       >  8:222489/1‑149 (MQ=255)
                                                        cATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACg                                                                                      <  7:222489/149‑1 (MQ=255)
                                                           gCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGc                                                                                                                 <  4:68537/119‑1 (MQ=255)
                                                           gCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGc                                                                                                                 >  3:68537/1‑119 (MQ=255)
                                                           gCACCTGTGCCAGCCATCGCCACAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCa                                                                                   >  5:101132/1‑149 (MQ=255)
                                                             aCCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATc                                                                                 >  7:50969/1‑149 (MQ=255)
                                                               cTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCg                                                                               >  3:232894/1‑149 (MQ=255)
                                                               cTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCg                                                                               >  3:219687/1‑149 (MQ=255)
                                                                  tgCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAAc                                                                                                                   <  4:120333/110‑1 (MQ=255)
                                                                  tgCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAAc                                                                                                                   >  3:120333/1‑110 (MQ=255)
                                                                    ccagccaTCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAgg                                                                          >  8:30632/1‑149 (MQ=255)
                                                                    ccagccaTCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAg                                                                           <  4:163506/148‑1 (MQ=255)
                                                                         caTCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTgcgc                                                                                                   <  6:199473/119‑1 (MQ=255)
                                                                          aTCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGc                                                                     >  4:260957/1‑148 (MQ=255)
                                                                              ccgcAACGCACCCCGGTGTTGTTCCAGGACGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGc                                                                >  8:174022/1‑149 (MQ=255)
                                                                                    cGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCgg                                                          <  6:101132/149‑1 (MQ=255)
                                                                                     gCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATACCGAATGCGTGTTCATCAGCGGGCCAACGCCGACg                                                                                                            >  8:268983/1‑98 (MQ=255)
                                                                                     gCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATACCGAATGCGTGTTCATCAGCGGGCCAACGCCGACg                                                                                                            <  7:268983/98‑1 (MQ=255)
                                                                                         cccGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTg                                                     >  6:289746/1‑149 (MQ=255)
                                                                                          ccGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTgg                                                    <  3:243203/149‑1 (MQ=255)
                                                                                          ccGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTgg                                                    <  7:174022/149‑1 (MQ=255)
                                                                                           cGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCCGGCCTGCGAAGCGGCCGGGCGTGGc                                                   >  2:320281/1‑149 (MQ=255)
                                                                                            ggTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGcc                                                  <  6:182654/149‑1 (MQ=255)
                                                                                               gttgttCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGt                                                >  5:32311/1‑148 (MQ=255)
                                                                                                ttgttCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGa                                              <  8:50969/149‑1 (MQ=255)
                                                                                                      cAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGt                                        <  2:132013/149‑1 (MQ=255)
                                                                                                       aGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTg                                       <  8:6106/149‑1 (MQ=255)
                                                                                                       aGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTg                                       <  7:30632/149‑1 (MQ=255)
                                                                                                        ggCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGc                                      <  1:149559/149‑1 (MQ=255)
                                                                                                         gCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGc                                      >  4:41914/1‑148 (MQ=255)
                                                                                                         gCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGc                                      <  3:41914/148‑1 (MQ=255)
                                                                                                         gCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATACCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTg                                       >  8:72609/1‑147 (MQ=255)
                                                                                                         gCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATACCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTg                                       <  7:72609/147‑1 (MQ=255)
                                                                                                                cGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGATTGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTa                               <  5:289746/148‑1 (MQ=255)
                                                                                                                cGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTAc                              >  7:126233/1‑149 (MQ=255)
                                                                                                                    ggCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGcc                           >  3:114021/1‑148 (MQ=255)
                                                                                                                      cACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCCGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAg                        >  4:145842/1‑149 (MQ=255)
                                                                                                                             ccAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGt                 <  4:219687/149‑1 (MQ=255)
                                                                                                                               agcagTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGtt                >  2:223976/1‑148 (MQ=255)
                                                                                                                               agcagTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGt                 >  7:307579/1‑147 (MQ=255)
                                                                                                                               agcagTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGt                 <  8:307579/147‑1 (MQ=255)
                                                                                                                                  agTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGTTGAtc            <  1:223976/149‑1 (MQ=255)
                                                                                                                                    tttGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGt                                                >  1:88408/1‑111 (MQ=255)
                                                                                                                                    tttGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGt                                                <  2:88408/111‑1 (MQ=255)
                                                                                                                                       gCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGTTGAtcgtcg        <  8:182195/148‑1 (MQ=255)
                                                                                                                                         cGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGTTGATCGTCGccc     <  8:126233/149‑1 (MQ=255)
                                                                                                                                         cGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGTTGATCGTCGccc     <  8:136894/149‑1 (MQ=255)
                                                                                                                                          gCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGTTGATCGTCGcccc    >  5:347774/1‑149 (MQ=255)
                                                                                                                                            aCGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGTTGATCGTCGCCCCCa  >  1:257235/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GAACCCTCGCGGGTACACCCGATCAACCATGTGGGCGAGCACTTCCAGGTGCCCGGCATGCACCTGTGCCAGCCATCGCCGCAACGCACCCCGGTGTTGTTCCAGGCCGGTGCGTCGGCACGCGGCCAGCAGTTTGCCGCACGGCATGCCGAATGCGTGTTCATCAGCGGGCCAACGCCGACGGTGCTGCGCCGCTACGCCGACGGCATCCGCCAGGCCTGCGAAGCGGCCGGGCGTGGCCGTGACGAAGTGCTGATCTACGCCCAGGCGCTGTTGATCGTCGCCCCCA  >  NC_002947/4000892‑4001180

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