Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,534,902 G→A 6.3% intergenic (+93/‑90) secA → / → argJ protein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,534,9020GA 6.3% 281.9 / 3.6 95intergenic (+93/‑90)secA/argJprotein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (40/48);  new base A (4/2);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

CGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGG  >  NC_002947/1534754‑1535050
                                                                                                                                                      |                                                                                                                                                    
cGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACgg                                                                                                                                                      <  5:373306/149‑1 (MQ=255)
 gTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACggg                                                                                                                                                     <  5:64071/149‑1 (MQ=255)
            gTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTtg                                                                                                                                          >  3:118176/1‑149 (MQ=255)
                gggCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTgg                                                                                                                                       <  3:93942/148‑1 (MQ=255)
                 ggCAAGAAGTTCAAGCATTGCCTCGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGGCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTgggg                                                                                                                                     <  8:355661/149‑1 (MQ=255)
                  gCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTggggg                                                                                                                                    >  7:13133/1‑149 (MQ=255)
                    aagaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGggcg                                                                                                                                  >  4:249753/1‑149 (MQ=255)
                    aagaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGggcg                                                                                                                                  >  5:56288/1‑149 (MQ=255)
                     agaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGggcgg                                                                                                                                 >  5:284679/1‑149 (MQ=255)
                     agaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGggcgg                                                                                                                                 >  5:345770/1‑149 (MQ=255)
                     agaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGagcgg                                                                                                                                 <  4:285183/149‑1 (MQ=21)
                        acttCAAGCATTGCCACGGTCAGATCGAGTGATCTTAACCTTACGTTGTACAAGCCTGGGGTCTCTTTGCGGCCCTTTCGCAACACACGGCCGCTCCCACAGGGACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCt                                                                                                                              <  5:216831/147‑1 (MQ=255)
                         gTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCtt                                                                                                                             <  8:211371/149‑1 (MQ=255)
                          ttCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTtg                                                                                                                            <  8:13133/149‑1 (MQ=255)
                          ttCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTtg                                                                                                                            >  2:308843/1‑149 (MQ=255)
                            cAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTtgtg                                                                                                                          >  5:287629/1‑149 (MQ=255)
                            cAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTtgt                                                                                                                           <  6:345770/148‑1 (MQ=255)
                             aaGCATTGCCACGGTCGGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTgtgt                                                                                                                         <  4:261588/149‑1 (MQ=255)
                              aGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                                        >  7:30468/1‑149 (MQ=255)
                               gCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcg                                                                                                                       <  6:287629/149‑1 (MQ=255)
                               gCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcg                                                                                                                       >  5:263703/1‑149 (MQ=255)
                               gCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcg                                                                                                                       >  1:57928/1‑149 (MQ=255)
                                cATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcgc                                                                                                                      >  4:270980/1‑149 (MQ=255)
                                cATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcgc                                                                                                                      >  4:275894/1‑149 (MQ=255)
                                cATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcgc                                                                                                                      <  4:276394/149‑1 (MQ=255)
                                    gCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATg                                                                                                                  >  5:28961/1‑149 (MQ=255)
                                     ccACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATgg                                                                                                                 <  6:56288/149‑1 (MQ=255)
                                      cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                                                <  3:322427/149‑1 (MQ=255)
                                      cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                                                <  3:7119/149‑1 (MQ=255)
                                      cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                                                >  4:7119/1‑149 (MQ=255)
                                       aCGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATgggg                                                                                                               >  1:237045/1‑149 (MQ=255)
                                         ggTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTg                                                                                                             >  4:87562/1‑149 (MQ=255)
                                               aTCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGc                                                                                                       <  1:2591/149‑1 (MQ=255)
                                               aTCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCAATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGc                                                                                                       >  5:152057/1‑149 (MQ=255)
                                                  gAGTGATCTGAACCTTACGTTGTCCAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCCACACCAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGC‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCgg                                                                                                     <  8:260623/148‑1 (MQ=255)
                                                   aGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGagc                                                                                                                                   >  5:94782/1‑117 (MQ=14)
                                                   aGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGagc                                                                                                                                   <  6:94782/117‑1 (MQ=14)
                                                           gAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCCCGATGTGGTGCGAAGCGGCCCCTTTTTc                                                                                           >  5:219089/1‑149 (MQ=255)
                                                            aaCCTTACGTTGTACATGCCTGGGGTCGCTTTGCGGCCCTTTCGCATCACACGGCCGGTCCCACAGGTAGCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCttttt                                                                                            <  3:94147/147‑1 (MQ=255)
                                                            aaCCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTc                                                                                           <  3:368668/148‑1 (MQ=255)
                                                            aaCCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTc                                                                                           <  7:295360/148‑1 (MQ=255)
                                                                  aCGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAcatc                                                                                     <  6:28961/148‑1 (MQ=255)
                                                                         aCAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttt                                                                             >  4:250641/1‑149 (MQ=255)
                                                                           aaGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttttt                                                                           <  3:275894/149‑1 (MQ=255)
                                                                           aaGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttttt                                                                           <  3:270980/149‑1 (MQ=255)
                                                                             gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTc                                                                          >  6:166860/1‑148 (MQ=255)
                                                                             gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTc                                                                          >  4:6543/1‑148 (MQ=255)
                                                                             gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCt                                                                         <  2:237045/149‑1 (MQ=255)
                                                                               cTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                       <  3:249753/149‑1 (MQ=255)
                                                                                tGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTtgtg                                                                                                                                        <  8:250951/83‑1 (MQ=255)
                                                                                tGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTtgtg                                                                                                                                        >  7:250951/1‑83 (MQ=255)
                                                                                   ggTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAgc                                                                   >  1:5645/1‑149 (MQ=255)
                                                                                    gTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCATTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAgcg                                                                  <  5:143941/149‑1 (MQ=255)
                                                                                       gCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGctc                                                               <  6:5490/149‑1 (MQ=255)
                                                                                            gCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAAt                                                           <  3:250641/148‑1 (MQ=255)
                                                                                                   tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATTA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGt                                                     <  3:160090/147‑1 (MQ=255)
                                                                                                   tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttg                                                   >  1:155973/1‑149 (MQ=255)
                                                                                                   tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGt                                                     >  4:160090/1‑147 (MQ=255)
                                                                                                     tCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttg                                                   <  5:166860/147‑1 (MQ=255)
                                                                                                     tCGCAACACAAGGCCGCTCCCACACGTACCGCGCCCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTATCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttg                                                   <  4:343232/147‑1 (MQ=255)
                                                                                                       gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                                                <  3:87562/148‑1 (MQ=255)
                                                                                                       gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                                                <  5:185496/148‑1 (MQ=255)
                                                                                                       gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                                                <  2:5645/148‑1 (MQ=255)
                                                                                                        aaaCAAACGGCCGCTCCCACCGGTACCGCGCTCGTCATTTACCACGCCATCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                                                <  1:316238/148‑1 (MQ=255)
                                                                                                         aaCACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttg                                             >  6:301811/1‑149 (MQ=255)
                                                                                                           cacaAGGCCGCTCCAACAGGTTCCGCGCTCGTCGATGT‑‑ACGGGCGCGATCTTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttg                                             <  6:219089/147‑1 (MQ=255)
                                                                                                            acaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                                                <  7:310024/79‑1 (MQ=255)
                                                                                                            acaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                                                >  8:310024/1‑79 (MQ=255)
                                                                                                             caAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTTGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTcc                                         >  3:153127/1‑149 (MQ=18)
                                                                                                              aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                                        >  5:363619/1‑149 (MQ=255)
                                                                                                               aGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCtt                                       <  3:6543/149‑1 (MQ=255)
                                                                                                                ggCCGCTCCCACAGGTACCGCGCTCGGCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCttt                                      <  8:296424/149‑1 (MQ=255)
                                                                                                                 gCCGCTCCCACAGGTACCGCGCGCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTg                                     <  8:30468/149‑1 (MQ=255)
                                                                                                                  ccGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGc                                    >  2:369941/1‑149 (MQ=255)
                                                                                                                     cccccACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCa                                 >  1:355360/3‑149 (MQ=21)
                                                                                                                        ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAc                                                                                        >  6:42453/1‑91 (MQ=255)
                                                                                                                        ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAc                                                                                        <  5:42453/91‑1 (MQ=255)
                                                                                                                        ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAccc                              >  3:162065/1‑149 (MQ=255)
                                                                                                                        ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAcc                               <  6:263703/148‑1 (MQ=255)
                                                                                                                        ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAc                                <  6:236407/147‑1 (MQ=255)
                                                                                                                        ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttggtc                                          >  6:51857/1‑136 (MQ=255)
                                                                                                                        ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttggtc                                          <  5:51857/136‑1 (MQ=255)
                                                                                                                          acaGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtct                                               >  8:258767/1‑130 (MQ=255)
                                                                                                                          acaGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtct                                               <  7:258767/130‑1 (MQ=255)
                                                                                                                          acaGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTg                            >  2:225290/1‑149 (MQ=255)
                                                                                                                          acaGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTg                            >  4:198298/1‑148 (MQ=255)
                                                                                                                            aGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCa                          >  5:200450/1‑149 (MQ=255)
                                                                                                                            aGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCa                          >  1:66487/1‑149 (MQ=255)
                                                                                                                              gTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                                        >  7:64071/1‑133 (MQ=255)
                                                                                                                              gTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                                        <  8:64071/133‑1 (MQ=255)
                                                                                                                               gatcgcgCCCGTTCATGG‑‑ACGAGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAccc                       >  7:243391/4‑149 (MQ=39)
                                                                                                                                aCCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCTATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAccc                       <  4:118176/148‑1 (MQ=255)
                                                                                                                                            cATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtttt          <  6:363619/149‑1 (MQ=255)
                                                                                                                                            cATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtttt          <  2:57928/149‑1 (MQ=255)
                                                                                                                                            cATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtttt          <  1:308843/149‑1 (MQ=255)
                                                                                                                                                a‑‑aCGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGaa       <  2:155973/148‑1 (MQ=255)
                                                                                                                                                   aCGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACt     >  7:48466/1‑149 (MQ=255)
                                                                                                                                                     gggCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCg   <  3:109868/149‑1 (MQ=255)
                                                                                                                                                     gggCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCg   <  8:48466/149‑1 (MQ=255)
                                                                                                                                                      ggCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCgg  <  6:284679/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGG  >  NC_002947/1534754‑1535050

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: