Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,534,898 A→G 6.3% intergenic (+89/‑94) 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,8980AG 6.3% 289.2 / 6.1 96intergenic (+89/‑94)secA/argJprotein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (40/48);  new base G (4/2);  total (44/52)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

AGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAA  >  NC_002947/1534751‑1535045
                                                                                                                                                   |                                                                                                                                                     
aGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGa‑‑a                                                                                                                                                    >  4:94147/1‑149 (MQ=255)
  ccGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACg                                                                                                                                                  >  6:64071/1‑149 (MQ=255)
   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                                                                                                                             >  5:56288/1‑149 (MQ=255)
                       aagaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGggcg                                                                                                                             >  4:249753/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:276394/149‑1 (MQ=255)
                                   cATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcgc                                                                                                                 >  4:275894/1‑149 (MQ=255)
                                   cATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcgc                                                                                                                 >  4:270980/1‑149 (MQ=255)
                                       gCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATg                                                                                                             >  5:28961/1‑149 (MQ=255)
                                        ccACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATgg                                                                                                            <  6:56288/149‑1 (MQ=255)
                                         cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                                           >  4:7119/1‑149 (MQ=255)
                                         cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                                           <  3:322427/149‑1 (MQ=255)
                                         cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                                           <  3:7119/149‑1 (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                                                                                                                              <  6:94782/117‑1 (MQ=14)
                                                      aGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGagc                                                                                                                              >  5:94782/1‑117 (MQ=14)
                                                              gAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCCCGATGTGGTGCGAAGCGGCCCCTTTTTc                                                                                      >  5:219089/1‑149 (MQ=255)
                                                               aaCCTTACGTTGTACATGCCTGGGGTCGCTTTGCGGCCCTTTCGCATCACACGGCCGGTCCCACAGGTAGCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCttttt                                                                                       <  3:94147/147‑1 (MQ=255)
                                                               aaCCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTc                                                                                      <  7:295360/148‑1 (MQ=255)
                                                               aaCCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTc                                                                                      <  3:368668/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                                                                                                                                   >  7:250951/1‑83 (MQ=255)
                                                                                   tGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTtgtg                                                                                                                                   <  8:250951/83‑1 (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                                           <  5:185496/148‑1 (MQ=255)
                                                                                                          gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                                           <  2:5645/148‑1 (MQ=255)
                                                                                                          gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                                           <  3:87562/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                                          <  7:258767/130‑1 (MQ=255)
                                                                                                                             acaGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtct                                          >  8:258767/1‑130 (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                                   <  8:64071/133‑1 (MQ=255)
                                                                                                                                 gTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                                   >  7:64071/1‑133 (MQ=255)
                                                                                                                                  gatcgcgCCCGTTCATGG‑‑ACGAGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAccc                  >  7:243391/4‑149 (MQ=39)
                                                                                                                                   aCCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCTATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAccc                  <  4:118176/148‑1 (MQ=255)
                                                                                                                                               cATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtttt     <  2:57928/149‑1 (MQ=255)
                                                                                                                                               cATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtttt     <  6:363619/149‑1 (MQ=255)
                                                                                                                                               cATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtttt     <  1:308843/149‑1 (MQ=255)
                                                                                                                                                   a‑‑aCGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGaa  <  2:155973/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAA  >  NC_002947/1534751‑1535045

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: